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Abstract 

This paper reviews quantum spin squeezing, which characterizes the sensitivity of a state with respect 
to SU(2) rotations, and is significant for both entanglement detection and high-precision metrology. We 
first present various definitions of spin squeezing parameters, explain their origin and properties for typical 
states, and then discuss spin-squeezed states produced with nonlinear twisting Hamiltonians. Afterwards, we 
explain pairwise correlations and entanglement in spin-squeezed states, as well as the relations between spin 
squeezing and quantum Fisher information, where the latter plays a central role in quantum metrology. We 
also review the applications of spin squeezing for detecting quantum chaos and quantum phase transitions, 
as well as the influence of dccoherence on spin squeezing. Finally, we review several experimental realizations 
of spin squeezing, as well as their corresponding theoretical backgrounds, including: producing spin-squeezed 
states via particle collisions in Bosc-Einstcin condensates, transferring photon squeezing to atomic ensembles, 
and generating spin squeezing via quantum non-demolition measurements. 
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1. Introduction 

In the past two decades, spin squeezing [1, 2, 3] has attracted considerable attention, both theoretically 
and experimentally. The notion of spin squeezing has arisen mainly from two considerations: The study 
of particle correlations and entanglement [4, 5, 6], and the improvement of measurement precision in ex- 
periments [2, 3, 7, 8]. However, the definition of spin-squeezing parameter is not unique, and we will here 
review several definitions that have been used in the past. The most widely studied squeezing parameters 
were proposed by Kitagawa and Ueda in Ref. [1], denoted by £|, and by Wineland et al. in Ref. [2, 3], 
denoted by The spin squeezing parameter £| [1] was inspired by the well-known photon squeezing; while 
the parameter £|j [2, 3] was introduced naturally in the standard Ramsey spectroscopy experiment, where 
the squeezing parameter is the ratio of the phase resolution when using a correlated state versus that when 
using a coherent spin state (CSS). 

One application of spin squeezing is to detect quantum entanglement [6, 9, 10], which plays a key role in 
both the foundations of quantum physics and quantum-information processing [11, 12]. Parameter £|, found 
to be related to negative pairwise correlations [13] and concurrence [14, 15], is able to characterize pairwisc 
entanglement for a class of many-body spin-1/2 states. Reference [4] proved that a many-body spin-1/2 
state is entangled if it is spin squeezed corresponding to ^ < 1. Indeed, spin-squeezing parameters are 
multipartite entanglement witnesses [6]. Another important reason for choosing spin-squeezing parameters 
as measures of multipartite correlations is that spin squeezing is relatively easy to be generated and measured 
experimentally [16, 17, 18, 19, 20, 21, 22, 23, 24, 25]. Spin-squeezing parameters only involve the first and 
second moments of the collective angular momentum operators, and in many practical cases, e.g., in Bose- 
Einstein condensations (BECs), particles cannot be addressed individually, and only the collective operators 
can be measured. To detect entanglement for realistic experiments, a class of spin squeezing inequalities 
were proposed [26, 27, 28], based the first, second, and even the third moments of collective spin operators. 
The use of spin-squeezing parameters as entanglement detectors has already been discussed in a recent 
review [6]. Besides, spin squeezing could also be useful for quantum computation (see, e.g., the reviews in 
[29, 30, 31]), as well as for quantum simulations (see, e.g., the reviews in [32]). 

Another application of spin squeezing, especially in experiments, is to improve the precision of measure- 
ments, e.g., in the Ramsey spectroscopy [2, 3, 33, 34, 35, 36, 37, 38, 8, 39], and in making more precise atom 
clocks [3, 5, 7, 23, 40, 41, 42, 43, 44] and gravitational-wave interferometers [45, 46, 47]. Therefore, many 
efforts have been devoted to the generation of squeezing in atomic systems. Basically, these works can be 
sorted by two categories. 

(i) Generating spin squeezing in atomic ensembles via atom-photon interactions [48]. Within this cate- 
gory, a very natural idea is to transfer squeezing from light to atoms [3, 16, 21, 49, 50, 51, 52, 53, 54, 55, 
56, 57, 58, 59, 60, 61, 62], and many proposals focused on transferring squeezing via electromagnetically 
induced transparency [63, 64, 65, 66, 67, 68, 69, 70, 71, 72]. When considering light-atom interactions, 
the detuning between the light and the atoms are very important. In the large detuning regime, an ef- 
fective Hamiltonian is obtained, consisting of a dispersive interaction term and a nonlinear interaction 
term, and the magnitude of the two terms can be adjusted, (a) The dispersive interaction between the 
light and atoms, results in a Faraday rotation of the polarization of the light [73] . Then by performing a 
quantum nondemolition (QND) measurement of the output light, the atomic ensembles can be squeezed 
conditioned to the measurement results [74, 75, 76, 77, 78, 79, 80, 73, 81, 82, 38, 83, 84, 85, 86, 40, 87, 
88, 89, 90, 42, 91, 92, 93, 94, 95, 23, 24, 96, 97, 25, 98, 99, 100, 101], and this method was realized in 
experiments [74, 63, 21, 23, 24, 102, 103, 104]. (b) The nonlinear term, which describes interactions be- 
tween atoms, is a one-axis twisting Hamiltonian [1, 105, 106], and can be used to generate spin squeezing 
[107, 108, 109, 110, 111, 22]. There are also some other proposals, including placing atoms in a high-Q 
cavity such that the atoms interact with a single field mode (not squeezed) repeatedly [112], illuminating 
bichromatic light on the atoms in a bad cavity [113] without requiring strong atom-cavity coupling. The 
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intrinsic spin squeezing in a large atomic radiating system was studied in Ref. [114], where spin-squeezed 
states are generated due to photon-exchanging induced strong interatomic correlations. Spin squeezing can 
also be produced via squeezing exchange between the motional and internal degrees of freedom[115]. In 
addition, squeezing can also be transferred from the atomic ensembles to photons [116, 117, 118]. 

(ii) Generating spin squeezing in the BEC via atomic collisions. In the last decade, generating spin- 
squeezed states in a BEC attracted many interests. The nonlinear atom-atom collisions in a two-component 
BEC can be described by a one-axis twisting Hamiltonian [4, 17, 119, 120, 121, 122, 123, 124, 125, 126, 127, 
128, 129, 130, 131, 132, 133, 18, 134, 135, 136], which can generate spin-squeezed states. Besides, two-axis 
twisting Hamiltonian can be realized via Raman processes [137, 138], or via two-component condensates in 
a double- well potential [139]. Bragg scattering induced spin squeezing was studied in Refs. [140, 141]. Spin 
squeezing can also be created in a spinor-1 BEC [142, 143], via spin-exchange interactions [144, 145], or free 
dynamical evolution [146, 147]. As studied in Ref. [148], the ground state of a F = 1 dipolar spinor BEC with 
transverse external magnetic field is spin squeezed. The squeezing of the ground state of spin-1 condensates 
in an optical lattice has been studied in [149, 150], where spin squeezing occurs in the Mott insulator phase. 
The particle collisions induced dccohcrcnce in the generation of spin squeezing was studied in [151], and 
spin-squeezed states generated in BEC are robust to dissipation [152], and particle losses [126, 153, 154]. 
Photodissociation induced spin squeezing was studied in Refs. [155, 156]. The spin-squeezed states in BEC 
have been used to perform sub-shot-noise measurements [46, 157], detect weak forces [158]. Experimental 
realizations of spin-squeezed states in BEC were reported in Refs. [17, 159, 160, 161, 18, 19, 20]. Recently, 
spin squeezing in BEC was re-examined by considering the effects of indistinguishable of particles [162, 163], 
where spin squeezing is neither necessary to achieve sub-shot noise accuracies in quantum metrology, nor 
in general an entanglement witness. Below, we only consider distinguishable cases. Besides the above 
two categories, other proposals such as using Ising Hamiltonian to generate squeezed states was studied in 
Ref. [41, 164, 165, 166]. 

This review is organized as follows. First, in Sec. 2, we present a full list of definitions of spin squeezing, 
which were proposed for different tasks, and give the relations between the spin squeezing and the bosonic 
squeezing. We make a comparison between bosonic squeezing and spin squeezing. We demonstrate their 
relations under the large-iV limit and small excitations. In this case spin squeezing will reduce to bosonic 
squeezing. 

In Sec. 3, we review the generation of spin-squeezed states via the one-axis twisting and two-axis twisting 
Hamiltonians, which have both fixed parity and particle-exchange symmetry. The proposal of using nonlinear 
Hamiltonians to generate spin-squeezed states originates from quantum optics, where the nonlinear Kerr 
effect is employed to prepare squeezed light. The dynamics controlled by the one- axis twisting Hamiltonian 
can be solve analytically. This type of Hamiltonian is widely studied in the regime of spin squeezing, 
and could be implemented in BEC [4, 17, 18, 19, 20], and large detuning atom-field interaction models 
[167, 107, 22], etc. The two-axis twisting has advantages in generating spin squeezing as compared to the 
one-axis twisting, however, it is not easy to implemented in experiments, and analytical results cannot be 
obtained for arbitrary system size. 

Since the nonlinear Hamiltonians involve two-body interactions, one can infer that spin-squeezed states 
may be pairwise correlated or entangled. This is reviewed in Sec. 4. We study the relations between spin 
squeezing and negative pairwise correlation. Then we review the relations between spin squeezing and 
entanglement, and this is one of the most important branches in the study of spin squeezing. The spin- 
squeezing parameter £| is closely related to pairwise entanglement. Moreover, the squeezing parameter 
proposed in [28, 15] has been found to be qualitatively equivalent to the concurrence for systems with 
exchange symmetry and parity. The parameter £^ can provide a criterion for multipartite entanglement. 
Thus, spin squeezing parameters can also be used to detect entanglement beyond pairwise entanglement. 
Spin-squeezing parameters can also be used to detect entanglement in a system of spin-j particles. In the end 
of this section, we also discuss the concept of two-mode spin squeezing, of which the bipartite entanglement 
is a valuable quantum information resource. 

In Sec. 5, we discuss spin squeezing and quantum Fisher information (QFI) in quantum metrology. The 
spin squeezing parameter £^ arises in the Ramsey interferometer, and the QFI lies at the heart of parameter 
estimation theory. Both quantities characterize the entanglement enhanced quantum metrology. Their 



4 



relations are reviewed and compared in the Ramsey interferometry. In Sec. 6, we discuss the applications of 
spin squeezing in quantum phase transitions (QPTs) and quantum chaos. In this section, we first employ 
the Lipkin-Meshkov-Glick model, which has a second-order QPT, to show that spin squeezing can be used 
to identify the critical point. Then, we review spin squeezing in quantum chaos. Two typical chaotic models 
are described: the quantum kicked-top model and the Dicke model. Some quantities, like entropy and 
concurrence, have been studied in these two models, and spin-squeezing parameters can also reflect their 
chaotic behaviors. The lifetime of a spin-squeezed state is much shorter when the system is chaotic than 
when the system is regular. 

In Sec. 7, we consider the effects of decoherence. Firstly, in Sec. 7.1, we discuss how decoherence affects 
the generation of spin squeezing. Then in Sec. 7.2 and 7.3, we discuss the robustness and lifetime of spin 
squeezing in the presence of decoherence. Finally in Sec. 7.4, we discuss how decoherence affects the ability 
of spin-squeezed states to the improvement of phase resolution in Ramsey processes. 

Finally, in Sec. 8, we review the generation of spin squeezing in experiments, including BEC, light- 
matter interaction, QND measurements. The closely related theoretical backgrounds are provided as well. 
For each physical system, we give a brief yet comprehensive summary of a series experimental results, as 
well as detailed discussions of several notable experimental improvements. We first discuss experiments on 
BEC, which involves Mott insulator transitions and the one-axis twisting. Then we review the transfer of 
squeezing, from light to atomic ensemble, via absorption of squeezed vacuum. Atomic ensembles can be 
viewed as quantum memories, and light is a very promising information carrier. Spin-squeezed states have 
been prepared in an optically thick Cesium atomic ensemble by using a squeezed vacuum. We also discuss 
a QND-measurement scheme, which generates conditional spin-squeezed states; and by utilizing feedback 
techniques, one can produce unconditional spin-squeezed states. Finally, we summarize some useful results 
and derivations in the appendices. 



2. Definitions of spin squeezing 

The definition of spin squeezing is not unique, and it depends on the context where squeezing is consid- 
ered. The most popular definitions are proposed by Kitagawa and Ueda [1], in analogy to photon squeezing, 
and by Wineland et al. [2, 3], in Ramsey experiments. The latter one is directly associated with quantum 
metrology. Besides these two widely studied definitions, there are some other definitions of spin squeezing, 
which were introduced for certain considerations. Before we discuss spin squeezing, we first give a brief 
review on bosonic squeezing, widely studied in quantum optics. These will give some initial indications on 
the possible definitions and the generation of spin squeezing. 

2.1. Coherent state and bosonic squeezing 

Let us first review some basic concepts of bosonic squeezing, which is widely studied in quantum optics 
[168, 169, 170, 171], and has applications in practical precision measurements [172, 173]. Squeezed states 
have also been studied in other contexts, including squeezed phonons in condensed matter physics [174, 175, 
176, 177] and squeezed states of microwave radiation in a superconducting resonant circuit [178]. 

Consider the bosonic creation and annihilation operators a) and a, with the corresponding commutation 
relation 

[o,ot]=l, (1) 

satisfying a\n) = \fn\n — 1) and a)\n) = \fn + l\n + 1), where \n) is the number state, i.e., a^a\n) = n\n). 
Below we consider two dimensionless operators X and P given by 

X = a + a\ P=°—^, (2) 
i 

which are the position and momentum amplitudes, respectively The commutator 

[X,P]=2i (3) 
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results in the Hciscnberg uncertainty relation 

MAP > 1, (4) 



where AA = y (A 2 ) — (A) 2 is the standard deviation, and (A^4) 2 is the variance. 

2.1.1. Coherent state 

The coherent state [179, 180, 181], usually called the minimum uncertainty state, satisfies the following 
simple condition 

AX = AP = 1. (5) 

The coherent state and its applications have been reviewed in [182, 183]. From the above criterion (5), the 
simplest coherent state is the vacuum state |0). A coherent state is defined as 

a\a) = a\a), (6) 

from which the first moments of X and P for the coherent state are directly obtained as Eq. (5). 

The uncertainty relation (4) cannot be violated. However, we can choose a state to make AX (or AP) 
smaller than 1, and this type of state is called coherent squeezed state. Generally, squeezing occurs when 
variance is less than 1 in any direction of the X-P plane. We introduce an operator in the X-P plane as 

Xg = e i0^a Xe -i9a^a 

= ae- w + aV e , (7) 

with X = Xq and P = X^/ 2 being special cases. We see that the generalized operator Xg can be obtained 
by rotating operator X in the phase space. The so-called principal quadrature squeezing [184, 185] is 
characterized by a parameter 

C * = ^) (AXe)2 ' (8) 

which is the minimum value of (AXg) 2 with respect to 0, and £^ < 1 indicates bosonic principal squeezing. 
The minimization in the above definition can be easily performed (See Appendix A) 

Cl = l + 2((aM-|(a)| 2 )-2|(a 2 )-(a) 2 |, (9) 

which only contains expectation values (a), (a 2 ), and (a* a). 

2.1.2. Generation of bosonic squeezed state 

A bosonic squeezed state can be generated by applying the following nonlinear Hamiltonian 

H = i (go) 2 - g*a 2 ) (10) 

on a coherent states as 

\a,r,)=S(r,)\a), (11) 

with squeezing operator S: 

S (77) = exp (-iHt) = exp {^l*a 2 - ^ flt2 ) ( 12 ) 

and complex number n = r exp (iff) , where r = —2 \g\ t. 

In Sec. 3.2, we will see that the two-axis twisting Hamiltonian that generates spin-squeezed states is 
inspired from this squeezing operator S. The nonlinear operator S mixes the X and P components and 
thus performs not only a rotation, but also squeezing. As shown in Ref. [170], the standard deviations for 
the rotated operators are derived as 

AX e/2 = e~ r , AP 0/2 = e r , (13) 
6 



The variance in Xg/2 is reduced with a factor e~ r , and the squeezing is determined by r, which is usually 
called the squeezing parameter. Computer animations showing various types of squeezed states, including 
squeezed wave packets are available online at [186]. 

Bosonic squeezed states can also be generated from the Kerr interaction Hamiltonian [187] 

H = n{a i a) 2 . (14) 

Choosing an initial state be the coherent state, we have the state at time t, 

\ip(t))KeTr = exp [~iKt{a)a) 2 ] \a). (15) 



2.2. Coherent spin state 

Now we turn to spin systems. Hereafter, we will mainly consider the system consisting of N spin-1/2 
particles, since in many practical cases wc deal with two-level systems or qubits. A two-level atom interacting 
with the radiation field can be treated as a spin-1/2 particle in a magnetic filed. The angular momentum 
operators for this ensemble of spin-1/2 particles are given by 

1 N 

J a =2*52 aion a = x iU> z (16) 
1=1 

where o~i a is the Pauli matrix for the Z-th particle. Before we discuss spin squeezing, we introduce the CSS. 
The CSS is defined as a direct product of single spin states 



N 



coe-lO^+e^sin-ll)! 



(17) 



where |0); and |1)/ are the eigenstates of o~\ z with eigenvalues 1 and —1, respectively. The above definition 
(17) gives an intuitive geometric description as shown in Fig. 1, i.e., all the spins point in the same direction. 
The CSS can also be written as 

|„) = |fl J 0) = (l + |t ? | a )-J J2 { 2 1 ) l V +m iJ», 1£C, (18) 

in terms of the Dicke states \j, m) , which are the eigenstates of J z with eigenvalue m, and 

Q 

T) = — tan — exp(— icfr). (19) 



Note that, when 6 = n, r\ is divergent, and only the coefficient before \j,j) is nonzero, thus \9 = tt, <fi) = 
The above form is expressed as the projection of a coherent state onto Dicke states \j, m), and m\r)}\ 
are related to the projection noise [188], obeying a binomial distribution, shown in Fig. 1. 
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(a) Coherent spin state 
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(b) Spin squeezed state 
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Figure 1: (Color online) Bloch sphere representations for a coherent spin state \CSS) shown in (a), and a spin-squeezed state 
| SSS) in (b). The radius of the Bloch sphere is \{J)\. The fluctuations of the spin components are represented by the circular 
multi-color disks in (a) and the multi-color elliptical disks in (b). The phase variance A(j>, determined by spin fluctuation disks 
and spin lengths, is the resolution of the spin state with respect to rotations, which characterizes the frequency precision in 
Ramsey spectroscopy. For a CSS, the projection noise is characterized by a binomial distribution, as shown in (a), while for a 
SSS, the projection noise is not binomial, and may be a sub-binomial distribution (b). 



Below, we shall derive Eq. (18). By using the identity 



exp[z(4(T + + f]U-)\ = cos i/£?7 + i — -^—{t i a + + ?ycr_), 



one can write the CSS in the following form 



where 



N N 

\8,<t>) = ®R0,4>Mi - (g)exp(Ca ;+ - C*<r,_) |0) ; , 
i=i i=i 



(20) 



(21) 



(22) 



It can be further written as 



R{0, 4>)\j,j) = exp(CJ+ -CJ-) 



(23) 



where \ j,j) = |0)/ is the eigenstate of J z with eigenvalue j — N/2, which represents all the spins polarize 

1=1 

to the z-direction. Here 

t/-J- J x i iJy 



are the ladder operators and R(6, 4>) is the rotation operator, which can be also given by 

R(9, (j)) — exp (—iOJft) — exp [i9 (J x sin <j> — J y cos <j))] , 

where n = (— sine/), cos 0, 0). 

Now we write Eq. (23) in a more explicit form by using the formula [189, 183] 

exp(C</+ - CJ-) = exp(?7J+)exp In (l + |?7| 2 ^ J z exp(-r?*J_), 

where r/ is defined as Eq. (19). The above relations can be verified as below. The rotation operator R{ 
can be expressed in a product form as 



(24) 
(25) 

(26) 



N 

i=i 

and Ri (8, <fi) can be evaluated readily as 

Ri(e,4>) = 

similarly, we find 



lexp(Ccr i+ -C*CTi-), 



i=i 



cos |C| -e^* sin |C| 
sin Id cos |d 



(27) 



(28) 



exp (t]<t + ) exp 



lWl + M a* 



exp (— rfcrJ) 



-1/2 



V" 



-v* (i + H 



-1/2 



i + M 




(29) 



Then, taking Eqs. (22) and (19) into the above two matrices, the relation (26) is obtained. 

2.3. Spin-squeezing parameters based on the Heisenberg uncertainty relation 

The definition of spin squeezing is not unique. When talking about spin squeezing, we should specify a 
certain spin-squeezing parameter. The uncertainty relation for angular momentum operators results from 
the commutation relation 

[J a , Jfj] = iEa^Jy, (30) 

where a, /3, 7 denote the components in any three orthogonal directions, and e Q ^ 7 is the Levi-Civita symbol. 
The uncertainty relation is 

(31) 



(AJ Q ) 2 (AJ^) 2 >|(J 7 )| 2 /4. 



In analogy to bosonic squeezing, spin squeezing could be defined when one of the fluctuations in the left-hand 
side of Eq. (31) satisfies [45, 190] (AJ a ) 2 < |(J 7 )| /2, and consequently, a squeezing parameter £jf is given 
by 



2 2(AJ a ) 



a 7^7 S (x,y,z) , 



(32) 
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and if C|f < 1? the state is squeezed. The subscript H refers to the Heisenberg uncertainty relations. More 
generally, one can define 

_ 2(AJ Tll ) 2 

iH ' ' (33) 

where n\ and n*2 are two orthogonal unit vectors. 

Below, we calculate t; 2 H for a CSS |0, <j>). According to the definition (32), the value of t; 2 H depends on 
the choice of the directions n\ and H2, and also on the parameter 9 and <fi of the CSS. From Eq. (B.6), one 
can find that 

e H = 1 (34) 

\n ■ n 2 \ 

More specifically, let no be the z-dircction and n\ the x-direction.The squeezing parameter then becomes 

sin 9 



cos< 



(35) 



where we have used Eq. (B.8). It is found that £jj may be less than 1 for a CSS, and this is not expected, 
since a CSS should not be spin squeezed. 

Now we discuss a generalized spin-squeezing parameter, introduced in Refs. [191, 192], for which spin 
squeezing could occur in two orthogonal directions simultaneously. The subscript H' indicates that this 
parameter is also related to the Heisenberg uncertainty relation. We first discuss spin squeezing in the 
x-axis, considering two orthogonal angular momentum operators in the y-z plane 

Jg = J y cos 9 + J z sin 6, 
Je+yr/2 = ~ Jy s i n + Jz cos 9, (36) 

and the commutation relations [J x , Jg] = iJg+n/2: [Jx, Je+Tr/2\ — ~iJe- The state is squeezed in the a:-axis 
if 

(AJ x f < 1 \(Je+n/2)\ , and/or (AJ X ) 2 < i \(Jg)\ . (37) 



Since the maximum value of (Jg) is \J (J z ) 2 + (Jy) 2 , we can say that the state is J^-squeezed if (AJ X ) 2 < 

1/22 2 l / 2 2 

I V {Jz) + {Jy) i otherwise, the state is J y -squeezed if (AJ y ) < 5 y (J z ) + (Jx) ■ The above two inequal- 
ities could hold simultaneously. 

More generally, we could write the above squeezing parameter as follows 

a = 2(AJ gl ) 2 

For the CSS \9, <j>), from Eq. (B.6), one can find that = y/T— {no ■ ni) 2 - More specifically, let no be the 
^-direction, we have 

&, = |sin0|. (39) 

The above equation indicates that the CSS can be squeezed, implying that is not a desirable definition 
for spin squeezing. 

2.4. Squeezing parameter ^| given by Kitagawa and Ueda 

Unlike the bosonic systems where the variance is equal in any direction for a bosonic coherent state, 
for a CSS the variance of spin operators depends on n, and there exists a prior direction: the mean-spin 
direction (MSD) 

s » -m-mv 
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where the second inequality results from the exchange symmetry. 

Below, we use ri± to denote the direction perpendicular to the MSD. For a CSS, we have (AJft ± ) 2 = j/2 
(See Appendix B), thus a state is spin squeezed if the variance of Jn ± is less than j/2. Compared with the 
principal squeezing (8), we arrive at the spin-squeezing parameter [1] 

_ min(AJjJ _ 4 min ( A Jjj 
Zs ~ jj- 2 - ^ , (41) 

where j — N/2, and n± refers to an axis perpendicular to the MSD and the minimization is over all directions 
n±. It is desirable that the spin-squeezing parameter £| is equal to 1 for the CSS. 

To calculate the parameter £|, the first step is to compute the MSD determined by the expectation 
values ( J Q ), with a G {x, y, z}. The MSD Hq can be written in spherical coordinates as 

no = (sin 9 cos <f>, sin 9 sin <fi, cos 9) , (42) 

where 9 and <fi are polar and azimuthal angles, respectively. The angles 9 and 4> are given by [193] 



9 = arccos 



\J\ 



arccos 
2tt — arccos 



(43) 



where | J] = J K ) 2 + (J y ) 2 + ( J z ) 2 is the magnitude of the mean spin. With respect to no, the other two 
orthogonal bases are given as 

Hi = (— sin^, cos0, 0), (44) 

712 = (cos 6* cos <j>, cos # sin 0, — sin 9) . (45) 

The above expressions are valid for 9 ^ 0,7r. For 9 = 0,7r, the mean spin is along the ±z direction, and the 
possible choices of <fi can be or ir. 

The second step now is to find the minimal variance of Jn ± = J • n±_ ■ The direction ri± can be represented 

as 

rij_ = Hi T = Hi cos ip + «2 sin <p, (46) 

where O is a 2 x 2 orthogonal matrix that performs rotations in the normal plane. The variance (AJ^) 2 
can be written as 

(AJ Hx ) 2 = (4 J = H ± THl, (47) 

where the symmetric matrix 

F=L ( J p T , Coy(J f J n2 )\ 
\Cov{J ni ,Jfi 2 ) (Ji 2 ) J 



in which 



1 



Cov(JfJ 1 , J,t 2 ) — — ([Jitn Jn 2 } + ) — fe) (Jk 2 ) 



( [>A"Ji 7 ^n 2 ] + ) J (49) 



is the covariance between and J^ 2 , and [X, Y] , = XY + YX is the anti-commutator. In the above 
equation, {Jn-s) = (Jn 2 ) = since Hi and n*2 are perpendicular to the MSD. The variance can be written as 

(AJnJ 2 =HiO T YOnl, (50) 
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and the matrix O can be chosen such that 



T TO = diag{A_, A + }, (51) 

where the eigenvalues 



(52) 



(Jl + 4 2 ) ± ^ (( Jl - J| 2 ) j + 4Cov ( J fll , J fl2 ) 2 
and min(AJ^j i ) z = A_ (See Appendix C), thus the squeezing parameters becomes 

(53) 



(4 + -4) ± J2 2 )) 2 + 4Cov(J Sl ,J S2 ) 2 



The optimal squeezing angle in Eq. (46) is given as 



I aim is 



3s ( 7 : A n2 ) if S < 0, 
7T- |arccos [ VA2+B2 ) if B > 0, 

where we define 

A =<-4-4.), B^CovfeJ^). (55) 

We know that the parameter £| = 1 for the uncorrelated pure CSS |6>, </>) in Eq. (17). Thus, if there 
are certain quantum correlations among the elementary spins, we may have £| < 1, i.e., the fluctuation in 
one direction is reduced, as shown in Fig. 1. Therefore, the squeezing parameter £| has natural connections 
with quantum correlations (entanglement). Indeed, it has been found that £| has a very close relation with 
quantities such as negative correlations [13] and concurrence [15], and we will discuss them in Sec. 4. 

2.5. Squeezing parameter £^ given by Wineland et al. 

Now, we discuss the spin-squeezing parameter proposed by Wineland et al. [2, 3] in the study of Ramsey 
spectroscopy. The squeezing parameter £^ is f ne ratio of the fluctuations between a general state and a 
CSS in the determination of the resonance frequency in Ramsey spectroscopy. The CSS here acts as a 
noise reference state. In contrast with £|, which is the analog of bosonic squeezing, the parameter £^ is 
substantially connected to the improvement of the sensitivity of angular-momentum states to rotations, and 
thus is attractive for experiments. 

Here we use a simple graphical way to describe this type of parameter, as shown in Fig. 1, while the 
mathematical form and essential physics are the same as for Ramsey spectroscopy and the Mach-Zender 
interferometer employed in Ref. [3]. More details of the Ramsey process will be discussed in Sec. 5.3. 

Consider now a spin state \ip). Without loss of generality, we assume the MSD to be along the z direction, 
and thus (J x ) — (J y ) — 0. This spin state is shown in Fig. 1, the circular or elliptical disks represent the 
variance AJ^, which is also called the projection noise. This state can be represented by a cone ending in 
a Bloch sphere. 

Let us now rotate the state around the x-axis. Then, in the Heiscnbcrg picture, we have 

J° ut = exp(i(f)J x )Jy exp(—i(j)J x ) = cos(j)J y — sin(j)J z . (56) 
From the above equation, we can immediately obtain 

(J° ut ) = -sin0(J z ), (57) 
(AJ° ut ) 2 = cos 2 0(AJ y ) 2 + sin 2 0(AJ 2 ) 2 - \ sin(2^)([J„ J z }+). (58) 
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According to the error propagation formula, Ax = Af{x)/ \ d(f(x)) /dx\, the phase sensitivity A<p can be 
calculated as 

AJ° ut _ AJ° ut 

<t> ~ \d(j^)/d^\ ~ icos^)r (59) 

For the rotation angle 0^0, (AJ° ut ) 2 ~ (AJ y ) 2 and cos</> ~ 1. Thus, the above equation reduces to 

a *=ra> (60) 

which allows us to know the phase sensitivity from the expectations and variances of the collective spin 
operators. 

From the above procedure, for the more general case that the MSD is not along the z direction, we can 
obtain the phase sensitivity as 

A^= AJ ^. (61) 
(J) 

For the CSS, we can obtain the phase sensitivity from Eqs. (B.8) and (B.9) as 

(A^)css = (62) 



which is the so-called standard quantum limit (SQL) or shot-noise limit. This is the limit of precision in 
atomic interferometry experiments when we use uncorrelated atoms. 

The squeezing parameter proposed by Wineland et al. is defined as [2, 3] 



(A0) 2 N(AJ R± f 



2 



(A0) 



CSS 



(J) 



2 



(63) 



This is the ratio of the phase sensitivity of a general state versus the CSS. Here, we choose the direction n± 
where AJ R± is minimized. This definition is related to the £| by Kitagawa and Ueda via 

Since j = y > \{J)\, we have £| < £ R . Even though these two parameters are similar (when j — \{J)\, 
£r = £l)> their physical meanings are different. When £ R < 1, the state is spin squeezed, and its phase 
sensitivity to rotation is improved over the shot-noise limit. According to Eq. (61), the phase sensitivity can 
be written as 

A0=-^L. (65) 



If d R < 1, A(f> < (A(j)) css beats the shot-noise limit. 

The lower bound of the phase sensitivity is given by the Heisenberg uncertainty relation (A J R± ) 2 (A ) 2 > 
l\(Jfi)\ 2 , from which we have £| 4(A J n , ± ) 2 /N > 1. Using the relations N 2 /A = j 2 > (J|) > {AJ n f and the 
fact that the largest eigenvalue of jH is j 2 , we obtain 

e R > (66) 



By using Eq. (65), we further have 

1 

iV' 
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A0>(A0) HL = -, (67) 



where (A0) HL is the Heisenberg limit (HL) [194]. 
Reference [2] proposed another definition, 

&« = -^?s = 2{A,hJ \ (68) 
\(J)\ \(J)\ 

We may choose another set of three orthogonal directions {n^, n^, n }- Then the uncertainty relation can 
be written as 

(AJ r ,J 2 (AJ sl ) 2 >|(J„- )| 2 /4. (69) 

Thus, the parameter £|j„ is naturally defined in this coordinate system. We will show that £jj„ = 1 for 
the CSS, indicating that this parameter is appropriate for characterizing spin squeezing. For a CSS, the 
variance of Jfj is zero, i.e., no fluctuations along the MSD (See Appendix B), (AJ^ ) 2 = 0. The variances 
of the angular momenta in the plane perpendicular to no behave similar to bosonic operators, for which we 
obtain (See Appendix B) 

(AJ Q ) 2 = (AJ^) 2 = j/2, (70) 

where the subscripts a and /3 denote two orthogonal axes perpendicular to the mean-spin direction Hq. So, 
we have = 1. 

The three squeezing definitions £g, £fp/, and £^ discussed above are not equivalent, however, they satisfy 
the following relation 

e s < e B » < &, 

since j > \(J)\ always holds. Therefore, if a state is spin squeezed according to criterion £|j < 1, it is 
definitely squeezed according to parameter t^,, and £fj, and thus £fj < 1 is the most stringent condition of 
squeezing among these three parameters in Eq. (71). 

We also note that, the projection noise can be schematically visualized in Fig. 1, where we show the 
projection noise distributions for a CSS and a SSS. For CSS, take \j,j) for example, and the MSD is along 
the z direction. If we measure J y , the expectation value (J y ) = 0, while the variance (AJ y ) 2 = j/2. The 
probability of the outcome states \ j, m) y obeys the binomial distribution 



P(m)= {j,j\j,m) y = {j,m exp{-i-J x ) ~ 2~^ \ N — m J 1 ^ 

where \j,m) y — exp(i^J x )\j,m) is the eigenstate of J y with eigenvalue m. The phase sensitivity A</> is 
determined by the width of the binomial distribution, and a spin-squeezed state may thus have a sub- 
binomial distribution [75]. 

2.6. Other spin- squeezing parameters 

In this section, we review some other definitions of spin squeezing, and discuss their applications. Ref- 
erence [4] proposed a spin-squeezing parameter which is a criterion for multipartite entanglement. 
The subscript B! indicates the close relationship between the spin squeezing parameters and £fj. The 
parameter £fj, is defined as 

(Jn 2 ) 2 + (Jn 3 ) 2 

For spin-1/2 many-body systems, it has been proved [4] that, if < 1 the state is entangled. However, 
an entangled state may not be necessarily squeezed. Based on this parameter, spin squeezing is directly 
connected to multipartite entanglement. Inspired by this finding, many works studied the relations between 
spin squeezing and entanglement (sec, e.g., [6]). The choices of the directions rii are arbitrary, and if n^ is 
chosen to be the MSD, while n\ is chosen to minimize the variance AJ^, then reduces to 
Below, we show that for the CSS, £fj, = 1. Substituting Eq. (B.6) into Eq. (73), we obtain 

& = - f";! 1 ' 2 , = i. (74) 

(no-n 2 ) +{n -n 3 ) 
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In summary, can be viewed as a generalization of and also provides a useful criterion for many-body 
entanglement. 

In the study of Dicke states \j, m) in BEC, Raghavan et al. [121] proposed a new kind of squeezing 
parameter 

e„ = N(AJa) \ , (75) 

where the subscript D indicates that this parameter can detect entanglement in Dicke states. If £f> < 1, 
the state is squeezed along ft. Actually, in Ref. [121], the authors only consider the variance (AJ z y alo ne 
the 2-axis. For a CSS \9,<p), £,% = 1, which can be directly proved by substituting Eq. (B.6) into Eq. (75). 
For all entangled symmetric Dicke states \j,m), we will show that this parameter can detect entanglement 
in Dicke states in Sec. 3. 

As discussed previously, an important application of the spin squeezing parameter is to detect entan- 
glement. This kind of entanglement criterion is based on collective-spin inequalities, and it is attractive 
for experiments, because in practice we cannot always address individual particles, while the expectation 
values and variances for collective spin operators are easier to measure, such as in population spectroscopies. 
Therefore, Tth et al. [28] generalized the spin squeezing definitions and gave a set of spin inequalities. We 
find that one of the spin inequalities is suitable to be rewritten as a new type of spin-squeezing parameter, 
and this inequality reads 

(A-l)(AJ Sl ) 2 >(j| 2 ) + (j| 3 )-iV/2, (76) 

which holds for any separable states, and the violation of this inequality indicates entanglement. This 
inequality can be further written in the following form 

N {AJ Hl f > (j 2 ) - (J Hl ) 2 - N/2, (77) 

and one can define a spin squeezing parameter related to entanglement 

A - , " (AJ "' 2 , . (78) 



whose form is similar to £f>. In fact, for symmetric states with only Dicke states populated, {J 2 

N/2(N/2 + 1), and then £|, reduces to £fj. Thus, for the CSS, = £fj = 1. In Table 1, we give a summary 
of the spin squeezing parameters discussed above, and we also show their values for the CSS state \9,<j>) in 
Eq. (17). 
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Table 1: This table shows the definitions for different spin-squeezing parameters. Except for £ff and £, 2 H i, the other spin- 
squeezing parameters all equal to 1 for the CSS. 



Squeezing parameters 


Definitions 


Coherent spin state 


References 


P 2 


2(AJ fll ) 2 
\(Jfh)\ 


1 - (n • ni) 2 
\(ft ■ n 2 )\ 


[45] 


P 2 
SH' 


2(AJ Sl ) 2 


y/l - (n • n^ 2 


[191] 


P 2 


2(A Jfi ± ) 2 nin 
1^)1 


1 


[2] 


P 2 


4(A Jfi ± ) 2 nin 
N 


1 


[1] 


& 


N 2 a 
4(J) 2 S 


1 


[3] 




A^(AJ fil ) 2 
(^ 2 > 2 + fe) 2 


1 


[4] 


e D 


iV(AJ s ) 2 
iV 2 /4 - (J a ) 2 


1 


[121] 


a 


NiAJftJ 2 
(j 2 )-N/2-(J Kl ) 2 


1 


[28] 



2. 7. Rotationally invariant extensions of squeezing parameters 

As discussed above, for a given state, the denominators of parameters £fj/, £f>, and £§; depend on the 
choice of the directions. Unlike parameters £i and for which the denominators are constants. Thus it 
is difficult to determine the minima of £fj, and Inspired by Toll's discussions [28, 195, 196], we 
will give slightly different definitions for P 2 ^, and £f,, and the new definitions provide us a simple way to 
determine whether a state is spin squeezed. 

Spin squeezing with respect to £f;, , £|j , and £f; are equivalent to the following three inequalities 

iV(AJ Sl ) 2 < (J S2 ) 2 + (J- 3 ) 2 , 
AT 2 

N (AJ r t) 2 < — (Jft) 2 , 

N (AJ a ) 2 < (J)' - y - (Jnf , (79) 

The right-hand sides of the above three inequalities can be written in a rotation-invariant form by adding 
(Jni) 2 to both sides of the first inequality, and (Jft) 2 to the second and third ones. Then we obtain 

(7V-l)(AJ„-) 2 + (jl)<(J) 2 , 

N 2 

(N-l)(AJ n ) 2 + (4)<— , 

(N - 1) (AJ„-) 2 + (4) < ( J ) 2 - ^ , (80) 
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which are equivalent to the original inequalities, and we replace ni with n in the first inequality. Note that 
the right-hand side of the inequality is invariant under rotations. Thus, to detect entanglement, we shall 
find the minimum value of the left-hand side by rotating the direction n. So, the spin squeezing parameters 
can be defined as 

p = min a [(N -1) (A Jnf + (4)} 

& = w> ™ n[( - N - 1} (AJs)2 + ( (82) 

,2 



mm 



m n [(N-l)(AJ H y + (4)] 



e E = - , : 2 ' — (83) 



( / y - n/2 



These parameters are rotationally invariant. 

The minimization procedure works as below. First the term (N — 1) (AJ^) 2 + ( J|) can be written as [28] 

(N - 1) (AJ n ) 2 + (jf ) = nTn T , (84) 
where the superscript T denotes the transpose, and T is a 3 x 3 matrix, which is defined as 

r=(JV-l)7 + C, (85) 

where the covariance matrix 7 is given as 

7W = C w -<Jfc)<Jj) for k,le{x,y,z} = {1,2,3}, (86) 

with a correlation matrix 

C w = i(J, J k + J k Ji). (87) 

The minimum value of nTn T is the minimum eigenvalue of Y (See Appendix B), and thus the spin-squeezing 
parameter based on these inequality can be defined as [197] 

> {fi) _ N/2 ^ 

If £|, < 1, the state is spin squeezed and entangled. In the case when (f 2 big) = j (j + 1) and when A m i n is 
obtained in the 7?j_-direction, we find £|, = and thus can be regarded as a generalization of For a 
CSS, A min =j 2 and || = 1. 

It is interesting to note that the left-hand sides of the above inequalities are equal to nYn T , and thus the 
minimum values are just A m ; n . Therefore, the new definitions of squeezing parameters become proportional 
to A m i n . Although these new parameters are not quantitatively equal to their original ones, they are 
qualitatively equivalent to the original ones in the sense that they all can detect whether a state is spin 
squeezed or not. 

In the above discussion, spin squeezing was defined for many-qubit states belonging to the maximum 
multiplicity subspace of the collective angular momentum operator J. This means that these states exhibit 
particle exchange symmetry. The concept of spin squeezing is therefore restricted to symmetric many-body 
systems that are accessible to collective operations alone. There is no a priori reason for this restriction 
in real experiments, thus Ref. [198] explored the possibility of extending the concept of spin squeezing to 
multi-qubit systems, where individual qubits are controllable in the sense that they are accessible to local 
operations. This requires a criterion of spin squeezing that exhibits invariance under local unitary operations 
on the qubits. 

Now we introduce a local unitary invariant spin squeezing criterion for N qubits. We denote the mean- 
spin direction for a local qubit i by 

a» = i7§\T (89) 
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Associating a mutually-orthogonal set {hi±, hi±_i, ftjo} of unit vectors with each qubit, we may define the 
collective operators 



N , JV , N 



J± = ^^Vi ■ n i± , J x < = ' ™ iJ -' ^0=0^0^- n i0 , (90) 

i—1 i—1 i—1 

which satisfy the usual angular momentum commutation relations as Eq. (30), and this leads to the uncer- 
tainty relation 

(AJ ± )(AJ ± ,)>±(Jo). (91) 

Analogous to the above spin definitions, like £| and we may define the corresponding squeezing param- 
eters [198] as 

^ S= N ' *« = {Jo f ■ (92) 

One important advantage of these squeezing parameters [198] is that they arc locally unitary invariant due 
to the definition of J in Eq. (90). Immediately, one can verify that for all TV-body separable states, the 
above two parameters in Eq. (92) equal to one, since the single spin-1/2 particle is never spin squeezed. 
This is in contrast with the original squeezing parameters £| and £fg, which are equal to 1 only for CSS. 

2.8. Spin-squeezing parameters for states with parity 

In preceding discussions, we reviewed several different spin-squeezing parameters presented in the litera- 
ture. In the following, we mainly focus on two typical spin-squeezing parameters, £| and A, which have wide 
applications in detecting entanglement, and in quantum metrology, etc. Besides these two parameters, we 
also consider parameter £|, in some parts of the paper, since £jj has close relations with £| and can be viewed 
as a generalization of Furthermore, ^ is significant in discussing spin squeezing and entanglement. 

We shall see that most typical spin-squeezed states are of a fixed parity. In this subsection we will study 
spin-squeezing parameters £|, £fj. and £g for states with even (odd) parity, and give relations among them. 
For general states, it is hard to find explicit relations of these three parameters. At first we explain the 
parity of a spin state and restrict ourselves to states with exchange symmetry and with only Dicke states 
being populated. For these states, we may define the parity operator 

p = {-\y h+J , (93) 

with eigenvalues p — 1 and —1 corresponding to even and odd parity, respectively. Thus the state \j, —j), 
which means each spin pointing to the — z direction, has even parity, and for a Dicke state \j, — j +n), if n 
is even (odd), the state has even (odd) parity. In general, if a spin state is spanned only by Dicke states 
of even (odd) parity, this state is of even (odd) parity. States with parity are very common, e.g., the state 
generated by the one-axis twisting model, which we will discuss in Sec. 3. 

2.8.1. Spin- squeezing parameters £| and ^ 

Here, we consider spin-squeezing parameters £| and For states with parity, we have the following 
relations, 

(J a ) = (J a J z ) = (J z J a ) =0, a = x,y, (94) 
which mean the MSD is along the z direction. The above equation leads to the following zero covariance 

Cov(j 2 , J H± ) = i ([J„, J fl J + ) - (J z ) (J a± ) = 0, (95) 

which implies that there are no spin correlations between the longitudinal (z) and transverse directions (x-y 
plane). The MSD is along the z-axis, then the general expression for the spin-squeezing parameter £| in 
Eq. (53) reduces to [15] 

e s = ^{{Jl + Jl)-\{J 2 -)\). (96) 
18 



Due to the exchange symmetry, we find 



(J a ) 


TV 
= 2 








(97) 


{Jl) 


N 

~~ 7 


N(N-l) 
f 4 


(O"l a £72o) 7 " = X,y,Z, 


(98) 




3N 

_ ~7~ 


N(N-1 
+ 4 


~ (01 ■ <7 2 ) , 




(99) 


(AJ Q ) 2 


JV 

~ 7 


1 + (JV-1) 


(o-i Q cr 2Q ) - TV (CTi a ) 2 








= JV 
4 


1 + NC aa - 


!ci Q cr 2a )] , 




(100) 






= (c r i a C72a) - 


- (fla) (0"2a) 




(101) 



where 

is the correlation function along the a-direction. Furthermore, we obtain 

(J2> =JV(JV-l)<<7i_o-a->, 

r2 , T 2\ ^ , N (N — 1) 



AT 

2" ' 4 
JV N(N-1) 
2~ + 2 



(102) 



(103) 
(104) 



By substituting Eqs. (102) and (103) into Eq. (96), we obtain 

^1 = 1 - 2(N- 1) (|<<7l_ff 2 _)| - (ff 1+ <7 2 _)) , 

where we have used (cxi + <7 2 _) = (o"i_<T2+) ; which results from the exchange symmetry. From Eq. (97) and 
the relation between the parameters £fj and we find 

N 



& 



2|J 2 



2 ,2 

2 _ SS 



{(Tlx 



\2 ' 



(105) 



The above expressions for the two spin-squeezing parameters establish explicit relations between spin squeez- 
ing and local two-spin correlations, and we will further discuss these relations in Sec. 4. 



Table 2: Comparison between different squeezing parameters for states with parity (extended from Ref. [197]). In the third 
column, simplified expressions are displayed for squeezing parameters for states with parity. The squeezing parameters are also 
expressed in terms of local expectations (fourth column). 



Parameters 


Definitions 


States with parity 


In terms of local expectations 


f 2 


4(AJfj x ) min 
N 


^{{jl + J 2 y)-\{J-)\) 


1 - 2(JV- 1) (K«7l_0-3-}l - (*U-*2-» 


f 2 


— s.<? 

4(J) 2 


e s 


e s 


4(J Z ) 2 /N 2 


(*1 Z ) 2 


f 2 


•^min 


min {Cs^ 2 } 


min{e s , 1 + (N- l)C zz } 


(J) 2 


i(J z ) 2 /N 2 


(Viz) 2 




AT2 


nun{f| >? 2 } 


min{£ 2 , 1 + (JV - l)C zz } 


a 


^min 


min {^I^ 2 } 


mixing, 1 + (JV-1)C,J 


(J 2 ) - JV/2 


4(J 2 )/JV 2 -2/N 


(1 - N- 1 )^ ■ 8 2 ) + JV" 1 
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2.8.2. Spin- squeezing parameter £ 2 



E 

Here, we derive explicit expression of £fj for states with parity. By using Eq. (95), the correlation matrix 
(87) is simplified to the following form 

'{Jl) C xy \ 
C xy {J 2 y ) , (106) 
(Jl)) 

where C xy = ( [J x , J y \ + )/2. From the correlation matrix C and the definition of covariance matrix 7 given 
by Eq. (86), one finds 

<N(Jl) NC xy \ 

NC xy N(J*) . (107) 

N {AJ Z ) 2 + (J z ) 2 J 

This matrix has a block-diagonal form and the eigenvalues of the 2x2 block are obtained as 

A± = f (<JJ + J, 2 >±|(^>|). (108) 
Therefore, the smallest eigenvalue A m ; n of T is obtained as 

A min =min{A_,iV(AJ z ) 2 + (J 2 ) 2 }. (109) 

It is interesting to see that the eigenvalue A_ is simply related to the spin squeezing parameter £| in Eq. (96) 
via 

Thus, from the definition of f| given by Eq. (83), we finally find 

~ minjfo,^ 2 ) 

el- /V J . (in) 

AUA /N 2 ~2/N 



where 



_2_ 4 



r = _(iV(AJ z ) 2 + (J z ) 2 ). (112) 

The meaning of <; 2 will be clear by substituting Eqs. (97) and (100) into the above equation. Then, we 
obtain 

e = N (AJ 2 ) 2 + (J 2 Z ) = 1 + (N - 1) C zz . (113) 

Parameter ^ 2 is just a linear function of the correlation function C zz given in Eq. (101). From Eq. (99) and 
the above expression, we write parameter £fj in terms of expectations of local operators as 

tE = (l-7V-i )(oW2)+7V -i- ( 114 ) 

In the special case when only spin j = N/2 is populated, i.e., (^J 2 ^ — N/2(N/2 + 1), we have (<xi • 02) = 1, 
and therefore, the above equation reduces to 

11 =min{^ ) c 2 }. (115) 

Thus, the relations among the spin-squeezing parameters £|, £fj, and are clear for states with parity, and 
in this case, the spin-squeezing parameters are determined by pairwise correlations in the z-axis and the x-y 
plane. The above results are summarized in Table 2 extended from Ref. [199]. 
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2.8.3. Dicke States 

Below, we first consider Dicke states and simple superpositions of two Dicke states. The Dicke states 
\j,m) axe entangled state except for m = ±j. Since we mainly consider an ensemble of spin-1/2 particles, 
the state \j, to) can be written as 



\j, to) 





u 


— to)! 


u 


+ 


m)!(2j)! 






— to)! 


u 


+ 


m)!(2j)! 



T3+ m \ 




j+m 

|1\®iV 



(116) 



Since the Dicke states \j, to) are eigenstates of J z with eigenvalue to, the MSD is along the z-axis. A single 
Dicke state has either even or odd parity. Thus, we can use formulas shown in Table 2, and we only need to 
calculate £| and q 2 . 

From Eq. (96), the squeezing parameter £| can be written as 

e s = i+j--({JD + (H7) 



j 

which is determined by two expectation values: (J 2 ) and (J_). For the Dicke states, we obtain 

e! = i + £^>i, (us) 

since j > to. The equal sign holds for to = ±j, when the Dicke state becomes a CSS. The quantity <; 2 is 
obtained as 



e = r 



N (AJ Z ) 2 + ( J z f = ^ < 1. (119) 



Thus, from Table 2, one finds 



a=fi = min(^ 2 )^ 2 <l, (120) 
indicating that the state is squeezed and entangled for m ^ ±j. The parameters £^ and are given by 



e R = ( ■' ) e s > i, (.121) 



711 

ll> = -4r2 = 1' ( 122 ) 

m z /j z 

thus spin squeezing according to the three parameters £|, and ^ 2 R , cannot reflect the underlying entan- 
glement in the Dicke states. 

2.9. Relations between spin squeezing and bosonic squeezing 

Above, we reviewed some basic concepts about bosonic and spin squeezing, and now we will demonstrate 
the relationship between these two mathematically distinct, yet intuitively connected squeezing [200, 201]. 
It has been shown in Ref. [201] that the spin-squeezing parameter £| reduces to the bosonic squeezing in the 
limit of large number of atoms and small excitations. For this purpose, we consider the principal quadrature 
squeezing defined as (8). The definition of Cfj provides an atomic squeezing counterpart to bosonic squeezing, 
and is similar to the definition of both of them searching for minimum squeezing. 

It is well-known that the Heisenberg-Weyl algebra describing the bosonic mode can be obtained by 
contraction from the SU(2) algebra describing the ensemble of atoms [202] . To see this, we define b = J_ / 
and fr! = J+/^/2j. From the commutation relation (1), we have 

[N^] = b\ W,b] = -b, [&,&t] = i__ (123) 

J 
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where Af = J z + j is the 'number operator', and its eigenvalues vary from to N, counting the number of 
excited atoms. In the limit of j — > oo and small (Af), the operators Af, b, and 6' satisfy the commutation 
relations of the Heisenberg-Weyl algebra. Note that, when we take this limit, the average number of excited 
atoms (Af) should be much less than the total number of atoms N. 
We can also use the usual Holstein-Primakoff transformation [203]: 

J+ = a^2j - at a , J_ =y/2j- at a a, J z =a i a-j. (124) 
In the limit of j — > oo, we have 

' a\ ^^a, -4^1, (125) 



\ Xl \ V J 

by expanding the square root and neglecting terms of 0(1/ j). We see that the bosonic system and the 
atomic spin system are connected by the large- j limit from an algebraic point of view. 

To display this connection, we consider even (odd) states. These states refer to those being a super- 
position of even (odd) Fock states for bosonic systems, and those being a superposition of Dicke states 
\n)j = \j,—j + n) with the even (odd) excitations for the atomic systems. The Dicke states \n)j satisfy 
Af\n)j = n\n)j. Specifically, even and odd bosonic coherent states have been realized experimentally in 
various physical systems. The even (odd) states serve as examples for demonstrating connections between 
bosonic and atomic squeezing. For even (odd) states, (a) = 0; thus, from Eq. (8), we obtain 

Cl = l+2<o t o> -2|(a 2 )|. (126) 

Obviously, one necessary condition for squeezing is that |(a 2 )| 7^ 0. 

For even (odd) atomic states, the squeezing parameter £| has already been given in Eq. (96), and we 
rewrite it as 

£! = l + 2(A0-^-^y^. (127) 

Using Eq. (125), in the limit of j — > 00, we find that Eq. (127) reduces to Eq. (126) for even (odd) states. 
Also, we may find that the squeezing parameter £|j also reduce to in this limit. This result displays a 
direct connection between bosonic squeezing and atomic squeezing. From an experimental point of view, 
the number of atoms is typically large enough, so the observed atomic squeezing is expected to approximate 
the bosonic quadrature squeezing. As a remark, Eqs. (126) and (127) obtained for even and odd states are 
also applicable to arbitrary states, which is discussed in Ref. [201]. 



3. Generation of spin squeezing with nonlinear twisting Hamiltonians 

In this section, we discuss generating spin-squeezed states with the one-axis twisting and two-axis twisting 
Hamiltonians. The one-axis twisting Hamiltonian is one of the most important models studied in generating 
spin squeezing, both theoretically and experimentally. It also describes a nonlinear rotator, and was studied 
in Ref. [204] before its applications in spin squeezing. The proposal of using these two types of twisting 
Hamiltonians to generate spin-squeezed states is directly inspired by using the nonlinear Hamiltonian (10) 
to produce bosonic squeezing. 

In Sec. 3.1, we first present the analytical results of the evolution of the one-axis twisting. Then we 
discuss how to implement this Hamiltonian in a two-component BEC, and by using large-detuned light- 
atom interactions. The experimental progresses are reviewed in Sec. 8.1. Then in Sec. 3.2, we discuss the 
two-axis twisted state. As compared with the one-axis twisting, one can obtain higher degree of squeezing by 
using the two-axis twisting Hamiltonian. However, this type of Hamiltonian is not easy to be implemented 
in experiments, and analytical results are not available for arbitrary system size. 
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3.1. One-axis twisting Hamiltonian 

Here we discuss the generation of spin squeezing by using the one-axis twisting Hamiltonian [1, 204]. 
The one-axis twisting model is very simple, and is one of the most widely studied models in generating 
spin-squeezed state [167, 107, 22, 119, 121, 123, 122, 205, 126, 127, 206, 129, 130, 207, 110, 157, 111, 131, 
132, 18, 22, 135, 208, 209, 210]. It works in analogy to the squeezing operator, Eq. (15), in photon system 
and has been implemented in BEC via atomic collisions [4, 17, 18, 19, 20] and in atomic ensembles [167, 
107, 111, 22, 210, 209]. Moreover, it allows simple derivations of various analytical results. 

3.1.1. One-axis twisted states 

Consider now an ensemble of N spin-1/2 particles with exchange symmetry, and assume that its dynam- 
ical properties can be described by collective operators J a , a = x,y,z. The one-axis twisting Hamiltonian 
reads 

N 

Hoat = xJx = \ ^2 a kx<rix, (128) 
k,i=i 

which is a nonlinear operator with coupling constant \ arL d involves all pairwise interactions, which indicates 
that the spin-squeezed states generated by this Hamiltonian may exhibit pairwise correlations. The most 
commonly used twisting Hamiltonian is along the z-axis, with the initial state being a CSS pointing along 
the £-axis. Since the twisting is along the x-axis, we choose the initial CSS along the z-axis to make our 
analysis consistent with the z-axis twisting version. Here we prefer the x-axis twisting Hamiltonian satisfying 

[H OATl P}=0, (129) 

where P is the parity operator given by Eq. (93). We choose the initial state as \j, —j) = |1)® . Considering 
its dynamic evolution, the spin-squeezed state at time t is formally written as 

|tt(t)> = exp(-i9J*/2)\l)® N , (130) 

where 

= 2 X t (131) 

is the one-axis twisting angle. This state is the one-axis twisted state with even parity, and the MSD is 
along the z direction, thus the results derived in Sec. 2 can be directly used here. 



Table 3: Expectation values of local observables for the one-axis twisted state. 





-cos^- 1 (9/2) 




i(l + cos^) 


(<7l + <7 2 -) 


1(1 -cos*" 2 *) 


(GTl_(72-) 


-1(1- cos N - 2 9)- 1 - sin (6/2) cos*" 2 (9/2) 
8 2 



The expectation values needed for calculating spin squeezing parameters are derived in Appendix D, and 
are summarized in Table 3. By substituting expressions (<T\-\-<T2-) and (<ti_<72_) into Eq. (104), we obtain 



e s = i - c r = i - (jv - i)c, 



where 



C = 



1 - cos*" 2 9) +16 siir 



cos 



2N-4 



1/2 



— II — COS 



N-2 , 



(132) 
(133) 
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It will be clear in Sec. 4 that the quantity C is the concurrence [14], measuring entanglement of two spin 
s = 1/2 particles. 

In the case that N > 1 and |0| < 1, while N |0| » 1 and N \9\ < 1, we can expand Eq. (133) and find 
that the spin-squeezing parameter scales as [1] 



at |0| = O = 12 1/6 {N/2)~ 2/3 . Since is very small and N is large, at O we find 



exp 



In 



-1/3 



(134) 



(135) 



and then, for large enough TV, we have 

£| - iV~ 2/3 . (136) 

Therefore, the projection noise is reduced by a factor of iV~ 2 / 3 . The optimal squeezing angle for the one-axis 
twisting is [1] 

(137) 



6 ~ -arctan(iV~ 1/3 ) , 



which varies with the particle number N; although if N is large enough, S is close to 0. 



0.2 
0.1 

*. 
-0.1 
0.2 




-0.2 



0.2 



0.2 
0.1 

<■ 
-0.1 
0.2 




-0.2 





x 



0.2 



0.2 
0.1 

-0.1 



-0.2 




-0.2 



0.2 



0.2 
0.1 

*. 
-0.1 
-0.2 




-0.2 0.2 

x 



Figure 2: (Color online) Husimi Q function of spin-squeezed states generated by the one-axis twisting model for various times: 
(a) = 0, (b) = 0.1, (c) \t = 0.2 and (d) \t = 0.3. The system size is TV = 60, and at the beginning x* = 0, the state is a 
CSS and thus the Q function is a circle. The optimal squeezing angle rotates with time. 



The squeezing and the dynamic evolution of the state can also be illustrated by calculating the Husimi-Q 
function 

|2 



Q(0 o ,0o) = K0 o ,to|¥(t)>r 



(138) 
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Table 4: Spin-squeezing parameters in terms of the rescaled concurrence C r given by Eqs. (132) and (133) for the one-axis 
twisted state. 



e s 


l-C r 


a 


1-Cr 




1-Cr 



which represents the quasiprobability distribution of \^ (t)), and \9 ,4>o) is the CSS given in Eq. (17). The 
Husimi-Q function is shown in Fig. 2, where the coordinate we used is determined by 

x = Q cos 6*o cos <po , 

y = Q cos 6*0 sin ^ . (139) 

As we can see that, the initial CSS is a circle given by the Husimi function, and during the evolution, the 
Husimi-Q function becomes squeezed and elliptical, while the squeezing angle rotates. 

Now we discuss another two parameters and Parameter £ R is easily obtained, as we know both 
£i and (cr lz ) 2 . To obtain we also need to derive C zz and (<Xi • CT2) as seen from Table 2. For this state, 
(<?i -172) = 1, thus the expression of reduces to 

£l = min{£U 2 } 

= udn{l-C r ,l + (N-l)C zz }. (140) 

By substituting expressions of (<Ji z <J2z) 2 and (ci z ) 2 (Table 3) into definition of the correlation function, we 
obtain, 

C zz = \{l + cos N - 2 9) - cos 2N - 2 {9/2) > 0. (141) 

The proof of the above inequality is given in Appendix C of Ref. [197]. 

As the correlation function C zz and the rescaled concurrence C r are always larger than zero, Eq. (140) 
reduces to 

fl = & = 1 - C r . (142) 

So, for the state given by Eq. (130), the spin-squeezing parameters £g and £,g are equal , and we summarize 
these results in Table 4. 

3.1.2. One-axis twisting with a transverse field 

In Ref. [119], it was found that a one-axis twisting Hamiltonian with a transverse control field 

H = X J 2 X +BJ Z , (143) 

is more effective in generating squeezed states, where B is the strength of the external field. Reference [211] 
proved that the optimally spin-squeezed states that maximize the sensitivity of the Ramsey spectroscopy 
are eigensolutions of this Hamiltonian. This type of Hamiltonian was considered in BEC [123, 122, 129, 130, 
132, 135], and atomic ensembles [111]. The initial state is also a CSS while the dynamic evolution 

cannot be solved analytically, except for N < 3, and numerical results show that the external field leads to 
an improvement of the degrees of squeezing in an extended period of time. 

3.1.3. One-axis twisting in Bose-Einstein condensates 

Below, we discuss how to derive the one-axis twisting Hamiltonian from a two-component BEC, which 
can be regarded as a BEC with atoms in two internal states, or similarly, a BEC in a double-well potential. 
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We consider the case of N atoms in two internal states and \B). The Hamiltonian of the system is given 
by [4, 125] (ft = 1) 

H= f dr [$ if ^* if + 9 -f^k if # if & (f) $k if) 

k=A,B 

dr g AB $ A (r) ^ B (r) i> A (r) 4> B (f) (144) 
where the single-particle Hamiltonian 

h k = ~ + V k (f) (145) 
Ira 

governs atoms in the internal state k, with atom mass m and trapping potential Vfc (r). The bosonic field 
operator i\) k if) annihilates an atom at position r in the internal state k, which obeys 

$k(f ,$t (?)]=*(?- (146) 

The interaction strengths 

fjki = , (147) 

m 

with dki the s-wave scattering length and g AA , 9b b, and g AB are for collisions between atoms in states A, 
B, and interspecies collisions, respectively. 

Now we employ the single-mode approximation for each of the components, 

VM (r) = <\> A if) a, rpB (f) = (j> B (f) b, (148) 

where <p A (t) and <p B if) are assumed to be real, a and b are the bosonic annihilation operators that satisfy 
[a, at] = [6, fet] = 1 and [a, 6] = 0. In the single-mode approximation, the Hamiltonian is rewritten as 

H = (u>a — Uaa) a) a + (w B - U BB ) tfb 

+ U AA (a! a) 2 + U BB (tfb) 2 + 2U AB a^atfb, (149) 

where 

w fc = J dr <p* k (r) h k <p k (r) , 

Uu = ^ [dr \fa (r)\ 2 \<t>t (r)| 2 . (150) 

Note that, <j> k if obeys the coupled Gross-Pitaevskii equation, thus u) k and U k i depend on the time t. Now 
the effective Hamiltonian is rewritten by using angular momentum operators via the following Schwinger 
representation 

J z = \ {a) a - 6+6) , 

J+ = a f &, J_ = ab\ (151) 

where the operator 2J Z measures the population difference between states \A) and \B), and J± describe 
the atomic tunneling between the two internal states. The total atom number operator N = a) a + b^b is a 
conserved quantity here. Using the Schwinger representation (151) the Hamiltonian (149) can be written as 

H = e(t)N + E(t)N 2 + 6 (t) J z + X (t) J 2 , (152) 

where the nonlinear interaction coefficient is 



X(t) = U AA + U BB -2U AB , 
2G 



(153) 



and the other coefficients are 



e (t) = {lu a +ljb- Uaa - U BB ) /2, 
E(t) = {U AA + U BB +2U AB )/4, 

5(t)=UJ A -L0B + (UAA-UBB)(N -1). (154) 

As the number operator N is a conserved quantity, if X (t) 7^ as presented in Sec. 8.1, spin-squeezed states 
are generated as discussed previously. The validity of using the single-mode assumption to calculate spin 
squeezing was verified in Refs. [120, 125], based on both Bogoliubov theory and positive- P simulations. 
Note that, the original Hamiltonian (144) entangles the internal and motional states of the atoms, which 
is a source of decoherence for the spin squeezing. With direct numerical simulation [4, 154], spin squeezing 
produced via the original Hamiltonian (144) is roughly in agreement with the one-axis twisting Hamiltonian 
(152). 

We now consider that a driving microwave field is applied [129], and the Hamiltonian becomes 



ft A (f)i>B(r)n R e- lAt + h.c. , (155) 



where A is the detuning of the field from resonance and £Ir is the effective Rabi frequency assumed to be 
positive. Following the same steps, we effectively obtain the Hamiltonian in terms of the angular-momentum 
operators as 

Hi =5{t)J z + X {t)J 2 z +n [j + e~' l5t + J_e tSt ] , (156) 

where 

n = Sl R fd?<l>* A (T f )(l> B (f). (157) 

Usually, we can assume 6 (t) = 0. As presented in Ref. [119], by using the control field to assist the 
one-axis twisting, spin squeezing can be maintained for an extended period of time. In Ref. [131], xM 
is assumed to be independent of time t, and the external field is turned off rapidly at a time £m, so that 
(t) = 0_r 0(£m — t), where 6 (t) is the usual step function, and the maximal-squeezing time £m is obtained 
analytically 

(l58) 

which is valid for large N (> 10 3 ). The time-dependent field f2(t) provides the control for storing the spin 
squeezing. 



3.1. 4- One-axis twisting from large-detuned atom-field interaction 

Next we discuss the derivation of the one-axis twisting Hamiltonian from a collection of two-level atoms 
interacting with a large detuned field [107, 108, 109, 110, 111, 22]. Consider an ideal model, a collection of 
TV two-level atoms laid in a cavity, the Hamiltonian for the whole system is 

H = H +Hj, (159) 

where 

Hq = wqJ z + oj c a?a, 

Hi = g (J+a + J-a)) . (160) 

Hq describes the free dynamics of the atoms and field, and Hj is the interaction term under the rotating 
wave approximation. The spin operators J z ,± describe the atomic system, a, oJ describe the cavity field, 
and g is the atom-field interaction strength. We can see from the Hamiltonian that, atoms interact with the 
cavity field, while no direct interaction exists between atoms. However, under the large detuning condition, 
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i.e. the atom- field detuning A = ujq — lo c is very large as compared to g such that |A| 3> g\N, and we can 
perform the Froich-Nakajima transform to obtain an effective Hamiltonian describing nonlinear atom-atom 
interaction. The Froich-Nakajima transform is performed as, 

H s = e~ s He s , (161) 

where S is of the same order as the interaction term Hi. Expand the above transform to the second order 
of S, we have 

H s = H + [H, S] + - [[H, S] , S] 

= H a + (Hj + [H a , S]) + - [(Hi + [H , S]) , S] + * [Hi, S] , (162) 
and let Hi + [Ho, S] = 0, which gives 

H s = Ho+ l - [Hi, S] , (163) 



S = ^(a s J--aJ + ). (164) 



where 

q — _ 
A 

Inserting Eq. (164) into Eq. (163), we obtain the effective Hamiltonian 

H s = Ho-ri[Jl- (2ot + l) J z ] , (165) 

and the factor rj — ,g 2 /A. Note that, the effective Hamiltonian (165) contains a one-axis twisting term J|, 
and a dispersive interaction term proportional to cvaJ z , which is employed in the QND measurement and 
shall be discussed in Sec. 8. 

If photon loss is taken into account, the effective Hamiltonian is of the same form as Eq. (165), while 
the interaction strength is modified to be [107] 

r,= 9 9 A — , (166) 
' A 2 + 7 2 /4 V ; 

where 7 is the decay rate of the cavity field and satisfies 7 <C A. The derivation of Eq. (166) can refer to 
Ref. [107], where the atoms are laid in a cavity that is highly detuned from the atomic transition frequency, 
while the atomic dissipation induced by spontaneous emission was considered to be negligibly small as 
compared to the time scale of interaction. Spin squeezing of atoms in cavity was also studied in [108, 109, 
110]. In experiments, the effective one-axis twisting Hamiltonian was demonstrated for squeezing individual 
high spin (F = 4 [22], and F = 3 [89, 111]) Cs atoms that interacts with off-resonant light field. 

References [207, 210, 209] proposed an interesting method to generate effective one-axis twisting Hamil- 
tonian via cavity feedback. The large detuned atom-field interaction in cavity induces an effective QND-type 
Hamiltonian, 

H = aa)aJ z , (167) 

while the backaction of the cavity light causes the photon number operator a) a to be linearly proportional 
to J z , thus the QND Hamiltonian becomes a one-axis twisting Hamiltonian. In a recent experiment [209], an 
effective one-axis twisting Hamiltonian was generated by cavity feedback, and achieved 5.6(6) improvement 
in signal-to-noise ratio for \F — 1, mp = 0) -H- \F = 2, mp = 0) hyperfine clock transition in 87 Rb atoms. 



3.2. Two- axis twisted states 

Although spin squeezing can be produced by the one-axis twisting model effectively, the optimal squeezing 
angle depends on the system size and evolution time. This problem is solved if the twisting is performed 
simultaneously clockwise and counterclockwise about two orthogonal axes in the plane normal to the MSD. 
The initial state is also \j, —j), and the twisting is about two axes in the 9 = ir/2, <fi — ±7r/4 directions. The 
relevant two spin operators are written as 
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J f , f = cos J x + sin J y = — (j, + J, 

1 

71 



(168) 



J f ,-f = cos (|) J * - sin (|) J y = ^ (Jx - Jy) • (169) 
The two-axis twisting Hamiltonian is written as [1] 

^TAT = Jl f - J|,_ f = x(^^ + JyJx) = |(4 - ( 17 °) 

which is analogous to the Hamiltonian (10) for producing squeezed light that creates and annihilates photons 
in pairs. By replacing a and a) with J^/^fN and J+/y/~N, respectively, we will obtain the two-axis twisting 
Hamiltonian shown in Eq. (170). Various approaches for implementing this Hamiltonian were studied in 
Refs. [137, 212, 213, 205, 138, 214]. The MSD is also along the z-axis. Unfortunately, the two-axis twisting 
model cannot be solved analytically for arbitrary N, except for N < 3. Below, we list two advantages of the 
two-axis twisting compared with the one-axis twisting case. 



Table 5: Comparison of squeezing parameters when using either one-axis twisting or two-axis twisting Hamiltonian. 





One- axis Twisting 


Two-axis Twisting 


Minimum £| and £^ 


1 

a 7V2/3 


1 

°< N 


Optimal squeezing angle 


S ~ — arctan 
2 


f N -i/s\ 




Unchanged 


Physical implementations 


(i) Bose-Einstein conden- 
sation [17, 18, 19, 20]; 

(ii) Large detuning 
atom-held interac- 
tion [107, 111]. 


Effective atom-atom in- 
teraction via photon ex- 
change [215, 212]. 



(i) The optimal squeezing angle is invariant during the evolution. Since the MSD for the two-axis twisting 
Hamiltonian is along the z-axis, according to Eq. (54), the optimal squeezing angle is determined by two 
quantities, — Jy), and (J x J y + J y J x ). Due to 

[J X Jy + JyJx, H TAT ] = 0, (171) 

(JxJy + JyJx) is invariant during the time evolution. Here, the initial state is the CSS thus 
(J x Jy + JyJx) — 0. Then, the optimal squeezing direction is <p = 0,7r/2 during the evolution. The Husimi 
function is shown in Fig. 3, and it is clear that the optimal squeezing angle is invariant. 

(ii) The degree of squeezing is high. By numerical calculations, the spin squeezing parameters scales 
as [1] 

^cx- (172) 

in the two-axis twisting model. Thus, according to Eq. (65), the phase noise approaches the Heisenberg 
limit. Comparisons between the one-axis twisting and two-axis twisting Hamiltonians are displayed in Table 
5. 
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Figure 3: (Color online) Husimi Q function of spin-squeezed states generated by the two-axis twisting model for various times: 
(a) %t = 0, (b) \t = 0.1, (c) x~t = 0-2, and (d) \t = 0.3. The system size is N = 60. Compared with the one-axis twisting, the 
optimal squeezing angle here is invariant during the time evolution. 



4. Spin squeezing, negative pairwise correlations, and entanglement 

In this and the following sections, we mainly concentrate on spin squeezing, quantum correlations, and 
entanglement. Here, we first consider spin squeezing and pairwise correlations in a system with exchange 
symmetry, and consider the squeezing parameters £| and Note that for a CSS, which has no pairwise 
correlation, the variance of Jn ± is evenly distributed on individual spin components in the n^-direction, 
thus if a state has negative pairwise correlations in the n^-direction, the variance J^ ± could be reduced as 
compared with CSSs. Below, we show that, as discussed in Ref. [13], spin squeezing with respect to £| < 1 
implies negative pairwise correlation in the nj_ direction. Moveover, the minimum pairwise correlation is 
associated with parameter and ^ < 1 is equivalent to the existence of negative pairwise correlation. 

4-1- Spin squeezing and pairwise correlations 

As discussed in Refs. [1, 13], the spin-squeezing parameter £| is well defined in the j = N/2 subspace, 
where the states are of exchange symmetry. Thus, the expectation values and variances of the angular 
momentum operators J a can be expressed in terms of the local pairwise correlations, which helps us to obtain 
the relations between spin squeezing parameters and the pairwise correlations. The pairwise correlation 
function is defined as 

GmjH = (o l f L a ] n] — {(Tm) (&jn) — (oW2fi) — (o~m) , (173) 

with i, j being particle indices, and the second equality holds due to exchange symmetry. One can further 
define the minimum pairwise correlation as 

G m = min Gm,2n, (174) 

n 

where the minimization is over an arbitrary direction n. 
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4-1.1. Spin-squeezing parameter £| and correlation G ft x ft± 

First, we consider the spin-squeezing parameter £|. From the relation J^j = 1/2 J^-, am, one finds the 
variance of Jf L1 which is given by 



( A ^) 2 = zE« 



t ^ \\0~inO~jn) — {o- l a){o-jn)) 

ij 

I [N (1 - (a is ) 2 ) + 7V(7V - 1)G w ,m] 



(175) 



where we used the symmetry property in deriving the last equality. From this equation, the correlation 
function can be written as 



G\n,2rt — 



N(AJ n ) 2 + (J.^ 



N 2 (N-1) N-l 
If the correlation function is along the direction n±, Eq. (176) reduces to 



Gift \ 



(N — 1)N 



(176) 



(177) 



where we have used the fact (Jft±) = 0. From the relation between the spin squeezing parameter £| and the 
expectation value (J~ ± ) , 

mm(4 ± ) = ^e s , (178) 



one obtain the minimum pairwise correlation [13] 



£-1 



mmG n± ,f l± = ^ — -, 

Til iV — 1 



(179) 



where we have used Eq. (177). This implies that min^ ± Gfi J _ft J _ < is equivalent to £| < 1. Therefore, a 
spin-squeezed state (£§ < 1) has negative pairwise correlation in the n± direction (Gn ± ,n ± < 0). 

4-1-2. Spin- squeezing parameter (fe and correlation Gm t 2n 

Now, we investigate the relation between the correlation function Gm,2R and the squeezing parameter 
It is more convenient to rewrite the correlation function as 

Gin,2R = (n T aia^n) — (ri T ai) (a^ n) = n T Gn, (180) 



where the normalized direction n = (n x , n y ,n z ) , and the pairwise correlation matrix G is given by 



(181) 



The matrix elements of G are 



Gk,i = (o-ik&2i) - (oife) (<7 2 i) , k,l = x,y, z. 



(182) 



From Eq. (176), we know that the correlation function Gm,2n is a linear function of the quantity 
N{AJ a ) 2 + (Ja} 2 , which can be written as 



N (AJ n ) 2 + (J a ) 2 = 



N 



(183) 
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The matrix 

n fii^h^y + ( j )( j)T (184) 

is just the matrix T given by Eq. (85). Therefore, we have 

N (AJ a ) 2 + (J n ) 2 = n T Tn. (185) 
Substituting Eqs. (180) and (185) into Eq. (176) leads to the following relation 

G = n*(n- i) " (N^Ty (186) 

where I is a 3 x 3 identity matrix, and thus G and T can be diagonalized simultaneously. This indicates 
that, Gfifi can be expressed in terms of £| or £jj by choosing a specific direction n. 

Now we look for the relation between the parameter ^ 2 E and the minimum pairwise correlation. From 
Eq. (180), the minimum pairwise correlation 

G rn = mm Gn,n = min (nGn T ) = g min , (187) 

n n 

where <? m i n is the minimum eigenvalue of the pairwise correlation matrix G. From Eq. (186), we find 



2 



ffmin - ^"'^ _^ • (188) 

In the symmetric case (j = N/2), as we have assumed here, £| = 4A m ; n /A^ 2 , and thus 

G - = §3i- ( 189 ) 

If the minimal pairwise correlation is in the plane normal to the MSD, the above relation will degenerate to 
Eq. (179). In another form, we write 

1% = 1 + (N - 1) G m . (190) 

This exact result indicates that the spin squeezing defined by the parameter £|, is equivalent to the negative 
pairwise correlation, i.e., G m < implies < 1 and vice versa. 



4-2. Spin squeezing and pairwise entanglement 

One of the most useful applications of spin squeezing is to detect entanglement for many-qubit sys- 
tem. To determine whether a state is entangled, we just need to measure the collective operators, which 
are particle populations in many cases. Moreover, in many experiments, such as BEC, particles are not 
accessed individually, and the spin-squeezing parameter is easier to obtain than the concurrence and the 
entanglement entropy. Different kinds of spin-squeezing inequalities may be used to detect various types of 
entanglement [26, 27, 28, 6]. In Refs. [13, 216, 15, 217, 218], the relationships between negative pairwise 
correlation, concurrence, and £| for symmetric states and even (odd)-parity states were found. References 
[219, 220] showed that for a two-qubit Dicke system, parameter £| is better than to measure entangle- 
ment. A multipartite entanglement criterion for spin-squeezing parameter was given in Ref. [4]. Inspired by 
this work, some other generalized spin squeezing inequalities were proposed. In Refs. [26, 27], by employing 
a positive partial transpose method, they found generalized spin squeezing inequalities as criteria for two- 
and three-qubit entanglement. In Refs. [221, 28], optimal spin squeezing inequalities were proposed. 

At the time when squeezing parameters £| and £p were proposed, Refs. [1, 2, 3] noticed the potential 
relationship between spin squeezing and entanglement. Since £g = 1 for CSS, Ref. [1] expected that for 
an appropriate correlated state, £| < 1. As shown above, the one-axis twisting Hamiltonian can produce 



32 



squeezed states, and this Hamiltonian involves pairwise interactions. This indicates that spin squeezing is 
associated with pairwise entanglement. 

To detect two-qubit entanglement, it was proven [26, 27, 222] that if the inequality 



(&) + 



N (N - 2) 



> 



(■4> + <4,>- 



N 



+ (N-l) 2 (J Hl ) 2 



(191) 



is violated then the state is two-qubit entangled. For symmetric states, the above inequality is simplified to 

(192) 



t 4 (J H ) 2 > 4(AJ S ) 2 



N 2 



N 



and a symmetric state is two-qubit entangled if and only if it violates the above inequality. 

Entanglement of two-qubit systems is characterized by the concurrence [14]. The concurrence C, quan- 
tifying the entanglement of a pair of qubits, is defined as [14] 

C = max(0,Ai - A 2 - A 3 - A 4 ) , (193) 

where the quantities A;'s are the square roots of the eigenvalues, in descending order, of the matrix 

Ql2 = Pl2(<riy ® V2y)p*2(°'ly ® a 2y), (194) 

where p\ 2 is the two-qubit density matrix, and p* 2 is the complex conjugate of pi 2 . From its definition 
(193), the concurrence C > 0, and the existence of two-qubit entanglement is equivalent to C > 0. 
For example, we calculate the concurrence of a pure state 



\ip) = a\00) + 6J01) + c|10) + 

Then density matrix of \ip) is 

and its conjugate p* is a pure state. Since <j\ v ® a 2y is a unitary operation, 

p = a\ y ® a 2y p* X 2 a \y ® < J 2y, 



(195) 
(196) 

(197) 



is still a pure state. Therefore, the rank of pp is less than or equal to one, and it has at most one nonzero 
eigenvalue A, and the concurrence 



C = = ^Tr (pp) = | (i>\a ly ® a 2y | 
= 2\ad- bc\ . 



(198) 



If a = d = l/y/2 and b = c — 0, then is the Bell state and is entangled since C = 1. 

Following the previous discussion, we proceed to give a quantitative relation between the squeezing 
parameter and the concurrence for states with even (odd) parity and nonzero mean spin, i.e., ( J z ) 7^ 0. The 
two-spin reduced density matrix for a parity state with exchange symmetry can be written in a block-diagonal 
form [216] 

(199) 



P12 



v + u 

U V- 



w y 



in the basis {|00), |11), |01), |10)}, where 

N 2 - 2N + 4(J 2 2 ) ± 4(J Z )(N - 1) 



v± 



4N(N - 1) 



N 2 - 4( J 2 ) 

w = u = 

4N(N- 1) 

J(J 2 X + J 2 )-2N 

y 4N(N-1) ' 



(Jl) 



N(N-l)' 



(200) 



33 



which can be written in terms of local expectations as 

v± = ~(l± 2(a lz ) + (a lz a 2z )) , (201) 

w = \ (1 - Wu<r 2 z)) , (202) 

u = {crx-a 2 _), (203) 

y=(a 1+ a 2 -). (204) 

The concurrence is then given by [14] 

C = 2max{0, \u\ - w , y - ^/v^vl} . (205) 

From the above expressions we know that, the spin-squeezing parameters and concurrence are determined 
by the expectation value {(J\ z ), correlations (o"i + o"2_), {a\-<r 2 -), and (p\ z <J 2z ). 

Since L J 2 ) = N/2 {N/2 + 1), from Eq. (200), one can straightforwardly verify that 

w = y. (206) 

Thus the concurrence given by Eq. (205) becomes 

C = 2max{0, \u\ -y,y- ^/vyvl} . (207) 

Since the density matrix should be positive, we find 

y/v^vl > \u\ . (208) 

From this inequality, one may find that if \u\ — y > 0, then y — y/v + v- < and if y — y?v+V- > , \u\ — y < 0. 
In other words, two quantities \u\ — y and y — ^Jv+V- cannot be simultaneously larger than zero. 



Table 6: Spin squeezing parameters £| and £|, as well as the concurrence C for parity states. Symbols are defined in the text 





Pairwise entangled (C > 0) 


Unentangled 


Concurrence 


C = 2{\u\-y)>Q 


C = 2{y- y /v+v-)>U 


C = 


P 2 


e s = l-(N-l)C<l 


e s >i 


e s >i 


F~ 2 


1% = 1 - (N - 1) C < 1 


i% = L-2(N-L)(y + ^v + vJ)G<L 


&>1 



From Table 2 and Eqs. (203) and (204), we can write the spin-squeezing parameter £| in terms of the 
reduced matrix elements u and y as 

e| = l-2(iV-l)(|u|-2/). (209) 

Again from Table 2, £^ contains the quantity 

f = 1 + (JV- l)C„. (210) 

So, to write l~% in terms of the matrix elements, we consider the correlation function C zz . From Eqs. (201), 
(202), and (206), we obtain 

y 2 -v+v_ = -hz zz . (211) 
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This is a key step. Thus, from Eq. (210), we have 

e = 1-4(JV -l){y + ^vZ){y- y/v^JT). (212) 
From the above equation and Table 2, we obtain 

1% = min {CI, 1 - 4(7V -l)(y + ^Q(v - V^^)} ■ ( 21 3) 

The relations between spin squeezing and concurrence are displayed in Table 6. From it we can see 
that, for a symmetric state, £jg < 1 is qualitatively equivalent to C > 0, implying that spin squeezing 
according to £|, is equivalent to pairwise entanglement [223]. Although £| < 1 indicates C > 0, when 
C = 2{y — y/v + V- ) > 0, we find £| > 1. Therefore, a spin-squeezed state (£| < 1) is pairwise entangled, 
while a pairwise entangled state may not be spin-squeezed according to the squeezing parameter £g [15]. 

Below, we give a simple example to illustrate the above results. Consider a simple superposition of Dickc 
states 

\ip D ) = cos 9\j : m) +e lv sin 9\j,m + 2), n = -j,...,j-2 (214) 

with the angle 9 £ [0, 7r) and the relative phase ip £ [0, 2ir), and j = N/2. This state is of even parity, and 
the MSD is along the z direction. The relevant spin-expectation values can be obtained as 

(J z ) = to + 2 sin 2 6>, 

(J 2 ) = to 2 + 4 (m + 1) sin 2 0, 

(^) = ie**'ain2fl v ^, (215) 

where /U TO = (j + to + 1) (j + m + 2) (j — to) (j — to — 1). With the above results we find 

e 1 ^ sin 29 , 



2N(N - 1) 

i r at 



JV-14 -/V 



^ [m 2 + 4(to + 1) sin 2 0] j , 



, \]{N 2 - 2N + 4 ( J|))2 - 16(JV - l) 2 (J z ) 2 

^ V+V - = AN(N 1) ' (216) 



and thus the spin-squeezing parameters are obtained as 

72 



e s = 1 - - | |sin26»| ^ - 2 [to 2 + 4 (m + 1) sin 2 0]- — \, (217) 
q 2 = 4 [m 2 + 4(to + 1) sin 2 0] - 4< ^ r T ^ [to + 2 sin 2 0l 2 . (218) 
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Figure 4: (Color online) Spin-squeezing parameters £| and concurrence C as a function of 8, with N = 3 and m 
The one-to-one correspondence between f|, and the concurrence C is very clear. 



N/2. 



In Fig. 4, we plot these two spin-squeezing parameters and concurrence versus 9 within one period. The 
system size N = 3 and m = —3/2. We observe that for 9 6 (0,7r/3) U (27r/3, n), £|, = £| < 1, therefore 
the state is spin-squeezed in the x-y plane. Moreover, as C > 0, the state is pairwise entangled. For 
9 6 (tt/3, 27t/3), it is obvious that the state is also pairwise entangled, since C > 0, while spin squeezing 
occurs in the z-axis since l~% < 1 and £§. > 1. These results clearly show that £|, < 1 is equivalent to C > 0. 

^.5. 5pm squeezing and many-body entanglement 

To characterize and detect multipartite entanglement is one of the most challenging open problems in 
quantum information theory [9]. The simplest multipartite state is the three-qubit pure state, for which 
there exists a good measure of tripartite entanglement based on the concurrence [224]. One can also use 
the state preparation fidelity F for a iV-qubit state p in order to investigate multipartite entanglement. The 
state p can be either pure or mixed. The fidelity F is defined as [225] 



where |^ G hz) = l/\/2(|00...0) + |H...l)) is the iV-particle Greenberger-Horne-Zeilinger (GHZ) state. The 
sufficient condition for TV-particle entanglement is given by [225] 



We have the freedom to choose other GHZ states such as |^ghz) = l/\/2(|00...01) ± |11...10)), etc. By 
local unitary operations we can transfer these states to the original GHZ state and these operations do not 
change the entanglement. Detailed discussions on this sufficient condition can be found in Refs. [225, 226]. 
They also discussed how to use many-body Bell inequalities to detect multipartite entanglement. 

Spin squeezing inequalities can act as entanglement criteria (see also Ref. [6]). To detect three-qubit 
entanglement, it has been proven that it is necessary to measure the third-order moments of the collective 
angular momenta [26, 27]. It was found that, for a state p and three orthogonal directions Hi, H2, and H3, 
if the following inequality [27] 



F{p) - (^GHZ|/0|*GHZ), 



(219) 



F(p) > 1/2. 



(220) 



N(N — 1)(5N — 2) 
24 



<0, 



(221) 
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is satisfied, the state p possesses a genuine GHZ-type entanglement. If one of the following inequalities 



iV-2/ n iV 2 -4iV + 8 /T . AW-2)(137V-4) n . . 

(2(Jl) + 2(JH 2 > - (J4)J (J S3 > + ^ i < 0, (222) 



2 

^■4) + (Jn 2 JuM - ^(4 3 ) + ^JrtJ + ^~ 2) < (). (22:$ i 



arc fulfilled, then the state p possesses a genuine 3-qubit entanglement. 

4-4- Spin-squeezing inequalities for higher spin-j systems 

The entanglement criteria presented above are all suitable for spin- 1/2 systems. However, in most 
realistic experiments, atoms can have larger spins. In general, we cannot straightforwardly generalize the 
above criteria to spin-j cases. To introduce the spin-squeezing inequality for higher-spin systems, we first 
consider the squeezing parameter £^ [227]. Unlike £|, it is difficult to determine the minimum value of 
since its denominator is not a constant. To find the minimum value of A, which corresponds to the extreme 
spin squeezing, one can use the Lagrange multiplier method, to find the state that minimizes [227] 

f = p(J z ) + (AJ X ) 2 . (224) 

Here the z-direction is assumed to be the MSD and AJ X is the minimum variance in the x-y plane. As 
discussed in Ref. [227], for integer spins the state that minimizes AJ X for a given ( J z ) has vanishing (J x ) and 
(J y ), thus / = p ( J z ) + ( J 2 ). To minimizes /, one could find the ground state of H = pJ z + J^., which is just 
a transverse-field one-axis twisting Hamiltonian. For half-integer spins, the problem becomes more difficult, 
since the state that minimizes (J 2 ) does not minimize AJ X for a given (J z ), and then one cannot formulate 
the problem as the diagonalization of an operator containing a Lagrange multiplier term. A Monte Carlo 
variational calculation that minimizes / is presented in Ref. [227] . 

Now, consider states of N spin-j particles. The collective spin operator is J = J^. J^, with Ji the spin 
operator for the i-th particle. A separable state of N spin-j particles could also be represented as 

Pscp = YtWh* ®Pk ] ®-® Pk N \ ( 22 5) 
k 

where p\ is the density matrix of the i-th spin-j particle. For separable states, the variance of J x obeys 
the inequality [227] 

l — , (226) 



N r i(k) N ( 

(AJ x ) 2 >Y /Pk J2[(^J x ) 2 \ i >£**Z>' F i( 

k i—1 k i=l \ 



where the first inequality comes from the concavity of the variance, the function Fj (x) is the minimum 
variance of J x divided by j for a given x, x — (J z )[ k ^ /j, and Fj (x) is a convex function. Therefore, by 
considering Jensen's inequality, we have 



{AJ x f > ^ Pk N jFj ( £ %-) > N jFj ( J>£ iJ 

k 



i=l J / \ k i=l J 



and if the above inequality is violated, the spin-j system is entangled. 

Another entanglement criterion is based on the following inequality [228], 



(AJ X ) 2 + (AJ V ) 2 + (AJ Z ) 2 > Nj, (228) 
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which holds for separable states of a multipartite N spin-j system. Based on this inequality, Ref. [228] gave 
a spin-squeezing parameter 

,a _ (AJ x f + (AJ y f + (AJ z f 

Ssinglet — j ; (ZZD) 

where J — Nj. The subscript 'singlet' means that this parameter can detect entanglement in singlet states. 
It has been proven in Ref. [228] that, states satisfying 

kinglet < 1 (230) 

are entangled. This inequality can be used to detect entanglement in the vicinity of many-qubit singlet 
states. These are pure states that are invariant under a simultaneous unitary rotation on all qubits. For 
example, for two qubits, the only state is the two-qubit singlet state, which is invariant under rotation. 



4-5. Two-mode spin-squeezed states 

Spin squeezing discussed above is also called one-mode spin squeezing, which describes the fluctuations 
of collective spin operators J a (a = x, y, z) of a total system. The two-mode spin squeezing [21, 229, 230] 
is analogous to the definition of the two-mode squeezing of continuous observables [144], and is proved to 
be a criterion of inseparability between spin variables of two separated atomic samples. 

A bipartite system is separable if and only if its state can be written as 

P = Y i PiPi ) <»P { i\ (231) 

i 

where pi denotes the probability. In continuous case, it has been proven that [231, 232] the following 
inequality 

A (V 1 ) + + A (p« - p^Y < 2 (232) 

is a sufficient condition for entanglement between subsystems 1 and 2. The continuous operators and 
pW belong to subsystems i, satisfying Up- k ', p^l = iS^i (k, I — 1, 2). 
For spin systems, the observables of interest are spin operators 

J( ± )=JW±J( 2 ). (233) 

The usual criterion for two-mode spin squeezing is [233] 

(AJi+)) 2 +(Aj(-)) 2 <(4 +) ). (234) 

It has been shown that [21, 230, 234], two-mode spin squeezing implies entanglement between spin compo- 
nents of the two subsystems. Furthermore, for pure states of two spin systems of equal dimension, two-mode 
spin squeezing after application of local unitaries is a necessary condition for entanglement, except for a 
set of bipartite pure states of measure zero [233]. Account for the Heisenberg-Weyl algebra discussed in 
Sec. 2.9, Eq. (234) can reduce to Eq. (232). In Refs. [21, 230, 234], two spatially separated atomic ensem- 
bles, each containing about 10 12 Cs atoms, were entangled for 0.5 ms via interacting with a polarized field. 
The entangled state generated in this experiment is similar to a two-mode squeezing but not the maximally 
entangled state. Besides, the two-mode spin-squeezed states can also be generated via QND measurement 
with feedback [235]. Due to the considerable long lifetime of entanglement, two-mode spin squeezing was 
proposed to be a valuable resource of quantum information, and could be used to perform atomic quantum 
state teleportation and swapping [229, 230]. 
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5. Spin squeezing, Fisher information, and quantum metrology 

As discussed previously, the squeezing parameter £^ characterizes the sensitivity of a state with respect 
to SU(2) rotations, and has been studied in quantum metrology. In this section, we first introduce the 
QFI [236, 172, 237, 238], which determines the precision of the parameter estimation. Then, we discuss the 
relation between QFI and spin squeezing [239, 240, 241, 242]. Finally, we discuss the applications of spin 
squeezing and QFI in entanglement enhanced quantum metrology, where the Ramsey and Mach-Zehnder 
interferometers are discussed. 



5.1. Quantum Fisher information 

The Fisher information measures the amount of information of a parameter that we can extract from 
a probability distribution. It determines how precise we can attain when estimate a parameter, and with 
a larger Fisher information we can estimate the parameter with higher precision. Firstly, we begin with a 
brief discussion about the Fisher information in probability theory, and for more information please refer 
to Chapter 13 in Ref. [238]. In the regime of the parameter estimation theory, the central problem is to 
estimate the parameter A in a probability distribution p(x\X), where x is the random variable, and below 
we only consider the single parameter case. 

The minimum variance of our estimation is determined by the Fisher information. Consider a general 
distribution p (x|A), to estimate the parameter A, we construct an estimator E (x) which is a map from the 
experimental data x to the parameter A. The expectation value of the estimator is 



(E(x)) = I dxp(x\X)E(x). (235) 
We consider the case (E (x)) — A, i.e., so-called unbiased estimation, thus we have 

dxp (x\X) [E (x) - A] = 0. (236) 
Differentiating both sides with respect to A gives 

dxp (x\X) L (x, A) [E (x) — A] = 1, (237) 



where 

L(x,X)= d ^^. (238) 
Now, square both sides of Eq. (237) and use the Cauchy-Schwarz inequality: 

l(/,5)| 2 <(/, f)(g,g), (239) 



we obtain the Cramer-Rao inequality: 



(AE(x)) 2 >^, (240) 



where 



is the variance of E (x), and 



(AE {x)Y = I dxp (x\X) [E (x) - A] 2 (241) 



1\ = I dxp(x\X)L(x,X) 2 (242) 



is the Fisher information with respect to A. N is the number of independent experiments (here N = 1), 
i.e., by repeating the experiment N times, the precision of E (x) is improved by 1/N. In practice, we need 
to maximize the precision of A, while this is fundamentally limited by the Fisher information I\, which can 
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be regarded as the amount of information of A that we can get from p(x\X), i.e., with more information, we 
can make more precise estimation. 

Now we turn to discuss the QFI [236, 237, 238], which is the extension of the classical Fisher information 
in the quantum regime. Consider an n x n density matrix p (A), by performing the positive operator valued 
measure (POVM) {Pi}, the parameter A resides in the outcomes probability as 

Pi (A) = Tr (p (A) Pi) . (243) 

According to Eq. (242), we have 



E (Rem[Ly(A)P i ]>r 
^ Pi (A) 



where L\ is the so-called symmetric logarithmic derivative determined by the following equation 

^=*[p(A)£ A + L A p(A)]. (245) 

The operator L\ is a quantum analogy of L (x, A) in Eq. (238). Thus, a set of POVM yields a corresponding 
probability distribution, which gives the Fisher information T\ of the parameter A. It was proven that [243] 

lA<FA=Tr[p(A)L 2 A ], (246) 

where F\ is the QFI, which is measurement-independent. The explicit expression of the QFI is given by 

Fx = ±^ + ± 2 Jf^\ { ^)\\ {k <n), (247) 
/'/ . / Pi i Pj 

»=1 ■> 

where pi and \ifi) are the ith nonzero eigenvalue and eigenvector of p(X). Therefore, they are all functions 
of the parameter A. The probability distribution that gives the maximum Fisher information, i.e. the QFI, 
is produced by the optimal POVM built of the eigenprojectors of L\. However, in some cases L\ is not a 
proper physical observable, and the corresponding POVM cannot be realized in experiments, thus the QFI 
is not always attainable. 

The role of the QFI in parameter estimation is given by the quantum Cramer-Rao bound [236, 237], 

A*)* > (AA)^ = — J^, (248) 

where N m is the number of independent experiments, A is the so-called unbiased estimator of the parameter 
A, i.e. (A) = A. Indeed, A is a map from the experimental data to the parameter space. The Cramer-Rao 
bound gives the ultimate limit for the precision of A that can be achieved. In a sense, parameter estimation 
is equivalent to distinguishing neighboring states along the path in parameter space. The QFI has a more 
intuitive geometric explanation, and it is the geometric metric of the state p (A) in parameter space, since 

F x (dX) 2 = 4(ds B ) 2 , (249) 

where (dsn) 2 is the Bures distance [244]. Thus, the Fisher information is equivalent to the so-called fidelity 
susceptibility [245], that has been extensively studied for characterizing QPT. 

5.2. Spin squeezing and quantum Fisher information 

Recently, it has been found that QFI gives a more stringent criterion for entanglement than the spin- 
squeezing parameter £^ [239, 241]. However, for identical particles, it was shown that [162, 163], neither spin 
squeezing nor quantum Fisher information is an entanglement witness. Below, we deal with distinguishable 
particles. Consider an ensemble of N spin-1/2 particles represented by a density matrix p. The QFI with 
respect to 6 is given by 

F[p(e),J H ]=Tr[p(e)Ll], (250) 
40 



where 

p (9) = exp (-iOJn) p exp (i9J H ) , (251) 
and Jfi is the generator of the rotation along the direction n. For pure states, the QFI becomes [239] 

F \p{6), J H ] = 4(AJ fl ) 2 . (252) 

For mixed states 

F [p (9) , J R ] = 2(pi ~ Pj) \{<Pi\J n \<Pi) I 2 ■ (253) 
ijtj Pi+Pj 

Comparing Eq. (247) with Eq. (253), there are no derivatives of pi, since the transformation in Eq. (251) does 
not change the eigenvalues of p. Actually, as explained in the previous section, F [p (9) , Jft] characterizes 
the geometric properties of p with respect to rotation, and is the rotational sensitivity of the state, thus we 
need not to perform a true rotation. 

According to Eq. (248), the lower bound of the uncertainty of 9 is given by 



where we set N m = 1. If we measure the angular momentum operator J, the lower bound of A9 becomes 

ml s = %, (255) 



where £^ i s the spin-squeezing parameter. Since the Cramer-Rao bound [236, 237] gives the ultimate limit 
of the precision of 6, we must have 

(A9) 2 QCB < (A9) 2 SS , (256) 



and thus 



which was proved in Ref. [239]. They also proved that if 

X 2 < 1, (258) 

the state is entangled. It is known that £fj < 1 also indicates entanglement. Thus, x 2 < 1 is more stringent 
than £jj < 1 in detecting multipartite entanglement, regarding the cases when \ 2 < 1 an d £k > 1. 

The inequality (257) can be readily proven when p (9) is a pure state, since in this case the QFI is just 
as Eq. (252). By using the Heisenberg uncertainty relation, we have 



F 



p{9),J n , (AJ n± f >\(J n )\ 2 , (259) 



where the directions n±, n'j_, and n are orthogonal to each other. Then, according to the definitions of £jj 
and x 2 j we can obtain the inequality (257). Both the spin-squeezing parameters and the QFI could be used 
as criteria for entanglement, which is a resource for high-precision measurements. Note that, generally, it is 
difficult for experiments to achieve the precision given theoretically by the QFI due to practical difficulties 
of the measurements. 



41 



Ramsey interferometer 




tt/2 pulse 





Free evolution: T 



tt/2 pulse 

■mm — 



Mach-Zehnder interferometer 

a 



Mirror 




lout 



2out 



Figure 5: (Color online) Schematic diagrams of the Ramsey (above) and Mach-Zehnder (below) interferometers. For the Ramsey 
interferometer, the blue arrows represent the states of the spin. The uncertainty of the spin components are represented by 
the elliptical multi-color disks. The red circle in (c) represent the uncertainty of a CSS. For the Mach-Zehnder interferometer, 
two beams injected from both sides of the first beam spliter, BS1. The parameter <f> is the phase difference gained between the 
two beam spliters BS1 and BS2. After the second beam spliter, BS2, photon numbers are detected, and the number difference 
between the two detectors is related to the phase difference 0. The two beam spliters correspond to the two iz/2 pulses in the 
Ramsey process (a) and (c). 

5.3. Spin squeezing, Fisher information and metrology 

Here, we discuss spin squeezing and QFI in quantum metrology [172, 246, 247, 248, 249]. In general, 
we cannot access parameters of a state or a Hamiltonian directly, since they are not physical observables. 
But in many cases, the parameter is related to an observable, like in the Ramsey spectroscopy, where the 
phase 4> gained in the free evolution is related to the population differences, and the parameters can also be 
obtained via Landau-Zener-Stckelberg interferometry [31]. This parameter estimation scheme is illustrated 
below, 

Prepare a state p, 
Quantum evolution p^ — U (4>) pU^ (</>) , 

Measurement (6) = Tr 6p^ , (260) 

where the information of the parameter <f> is contained in the measurement result (0}<£, and in some cases 
the observable O cannot be accessed directly, and this was studied in Ref. [250]. The fluctuation of the 
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observable O is unavoidable, and the variance of <j> is 

(AO), 



A0 



d{0)^/dct> 



(261) 



Since is periodic, its mean and variance are not calculated as non-periodic variables, and this should be 
declared. To calculate the mean of 4>, we first obtain the mean of O, that is (O)^. Then the mean of </> is 
the inverse function of (O)^. For more details, please refer to Ref. [251]. Since in many cases the evolution 
processes are fixed (or hard to change), we can reduce the variance of <fi by choosing an appropriate initial 
state p and operator O. When we measure J z , the precision is given by squeezing parameter £ R as shown 
in Eq. (65). 

We first give a brief summary of Ramsey processes, shown in Fig. 5. Consider an ensemble of N two-level 
particles interacting with an applied magnetic field B. The Hamiltonian of the two- level particles 

H = -p,-B, (262) 

where \i = (IqJ is the magnetic moment and B — Bori z + B\ 1 with a static magnetic field 

B = -Hojq/pq, pq < 0, (263) 

and a time-dependent field 

B\ = B\ [n x cos(wi) + n y s'm(ujt)] . (264) 

The Hamiltonian now becomes 

H = hu Q J z + ^-^ ( J+e" 4 "* + h.c.) , (265) 

where Q R = \/iqBi \ jh is the Rabi frequency. It is convenient to use a frame of reference rotating around Bq 
with frequency uj, 

H R = h(uj Q - lS) J z + MI R J X1 (266) 

where the second term in Eq. (266) acts as a pulse. When this pulse is applied, the above Hamiltonian is 
approximated as 

h r ~ nn R j x , (267) 

since Q R ^> \luq — u\. 

As shown in Fig. (5), the Ramsey interferometry consists of two 7r/2-pulses of length t^/ 2 — tt/ (2il R ) 
and a free evolution of length t. The first pulse plays as a tt/2 rotation around the x-axis. Thereafter, during 
the free-evolution period, when B\ = 0, the state vector precesses about the z-axis and acquires a phase 
<p = (cj — ujq) t. Assuming t 3> the time of the entire process is 

t f = 2t w/2 +t ~ t. (268) 

After the free period, followed by a second Ramsey pulse, the spin direction is rotated around the x-axis by 
7r/2. The the initial state evolves to 

\iP(t))=UW(0)), 

where the unitary operator is 

U = exp (-i-T^Jx) ex P {i4>Jz) exp (~^^) 

= exp (—iipjy) exp (— in J x ) . (269) 

Then we could measure the number of atoms in the excited energy level |0) to estimate 4>. This is equivalent 
to measuring J z : 
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(J z ) t = (Jx) t=0 sm(t) - (J z ) t=0 cos(/), (270) 

where the subscript t denotes the time, and in the right-hand side, the average (-) f=0 is carried out using 
the initial state. From the evolution operator U , we can further find the variance 

(AJ z f t = cos 2 (AJ z 2 ) t=0 + sin 2 cp (AjJ) t=Q - sin(20)Cov(J x , J z ) t=0 . (271) 

Then ujq is related to the measurement data of ( J z ) via 

1 „ ( (Jz) t 



uq = lj + - arccos I — J . (272) 

In practice, we repeat the above procedure N m times, with a total experimental time T — N m t. Then, we 
can estimate 4> with the uncertainty 

(A0) 2 = {AJz ^ s , (273) 

N m \d{J z ) t /dij>\ 2 ' 

which is obtained by using the propagation of the fluctuation (261). Below, we set N m = 1. 

From Eqs. (270) (271) and (273), the phase uncertainty is related to the initial states and the phase. 
Below, we discuss the uncertainty of <fi for different initial states. If the initial state is prepared as a CSS 

|V-(0)>css = |i,-i) = |l>^, (274) 

i.e., all the atoms are prepared in their ground states and the spin vector points in the — z direction, as 
shown in Fig. 13. For this initial state shown in Eq. (274), we can find 



which is the shot-noise limit, and we have already known in Sec. 2.5 since = 1 for a CSS. We can also 
obtain Fj, = 1 by using Eq. (247). 

Next, we consider a spin-squeezed state 



1 



|V>(0)>SSS = ^ 



li,o) x --^di.+i), + \j,-i) w ) 



(276) 



where the subscript x denotes the state is the eigenstate of J x , i.e., J x \jim}x — m \ji m )x- This state was 
studied in Ref. [212], the expectation values are 

(^)* = o = (^) t= o = Cov(J :c ,J z ) = ) 



and 



This state is spin squeezed, since 



(A J z f t=0 = + 1) - - A . (278) 



^WTv (279) 
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fitf + i) 


~ 4' 
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~ 4' 



and the phase uncertainty is optimized at cj) ~ 7r/2, 



which attains the Heisenberg limit. 

Now, we take the initial state to be a maximally entangled state 



1 
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|V>(0))GHZ = ~f= (\j,.])y + \j, -j)y) , (281) 



which is the A-body GHZ state, or the NOON state in quantum optics. We emphasize that, for this state 
the phase cannot be estimated by measuring the population difference J Zl since 

Wh = Wh = °. ( 282 ) 

that means we cannot get any information about <fi, and the spin squeezing parameter £p is divergent. On 
the other hand, after the Ramsey process 

|V>W>))ghz = W0))ghz = (|j, i>y + e-'^jj, -j>„) . (283) 

By using Eq. (247) we find its QFI Fs = A 2 , which implies that the precision of cf> could attain the Heisenberg 
limit. To achieve this precision, we should measure the parity operator [33, 37, 252] 

N 

P=H<r lz . (284) 
i=i 

The expectation values and variances of P under the state \ip (</>))ghz are 

(P) t = cos N<f>, (P 2 ) t = l, (285) 

thus the optimal phase uncertainty 



( A< ^)ghz = ^2 (286) 

is attained when N(j> ~ tt/2, which is the Heisenberg-limit uncertainty. Although the Heisenberg-limit 
precision is attained, we should note that, the parity operator P shown in Eq. (284) is not easy to measure 
in experiment as compared with the angular momentum operators, especially when A is large. Additionally, 
it is also difficult to prepare an A-body GHZ state with nowadays techniques. In summary, both spin 
squeezing and QFI are related to the precision in phase estimation. Metrology based on spin squeezing is 
comparatively easy to implement in experiments, while the ultimate precision is determined by the QFI. 

In the end we give a brief review about the Mach-Zchnder interferometer, which can be regarded as 
an optical version Ramsey interferometer. It consists of two beam spliters (BS) and two mirrors, and 
the schematic diagram is shown in Fig. 5. Light beams passing through the Mach-Zehnder interferometer 
undergo transformations as 

^ "Wbs ( ? in ) > (287) 



1 



l2out / \ «2i 

where dj (i = 1, 2) are the annihilation operators of the ith path, and 



denotes the transform of the BS, R = sin (a/2) and T = cos (a/2) are the transmission and reflection 
rates, respectively. The image number i arises from the half-wave loss. Below, we consider the 50-50 BS, 
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that is a — n/2, which acts as the two ir/2 pulses in the Ramsey process. The difference between the two 
path lengths gives rise to a phase difference <j), and can be represented by acting the transform 




(289) 



on modes a\ and a-i. It is convenient to manipulate these transformations in the Schwinger representation 
(151) as 

Ja,out = (oIout>4out) Y ( alTnt ^ = UJa ' inU ^ ^ = X ' V > ^ ' ^ 290 ^ 

where the unitary transform U is just as Eq. (269). At the output-port, the information of the phase <j> is 
obtained via detecting the photon number difference. Since the Mach-Zehnder interferometry is an optical 
instrument, it is not easy to prepare well-defined number of photons, unlike the Ramsey interferometer, 
where the atom number is conserved. 



6. Spin squeezing, quantum phase transitions, and quantum chaos 

In this section, we discuss the applications of spin squeezing in identifying QPTs [253] and quantum 
chaos [254]. For this task, spin squeezing mainly has three advantages: (i) It is comparatively easy to be 
measured, (ii) It is an entanglement witness, (iii) It characterizes the sensitivity of a state with respect to 
SU(2) rotations, thus is promising to detect quantum chaos. 

The QPTs occur at absolute zero temperature, and is driven purely by quantum fluctuations. Conven- 
tionally, QPTs were first studied by Landau's order parameter theory in the framework of statistics and 
condensed matter physics. When a QPT occurs, the ground state can change drastically, and the correla- 
tion length diverges. Therefore, considering these intrinsic properties of QPT, researchers investigated it 
by using concepts borrowed from quantum information [11, 12], such as quantum entanglement [9, 6] and 
fidelity [255].. As discussed previously, spin squeezing is closely related to entanglement, and it is easier to 
measure experimentally. Therefore, it is desirable to study spin squeezing in QPTs. 

Then we discuss using spin squeezing as a signature of quantum chaos. In classical regime, one of the 
most distinct feature of chaos is the extreme sensitivity of the trajectories with respects to perturbation, 
however, in the quantum world due to the unitarity of quantum evolutions, the overlap (or fidelity) between 
two initially separated states is invariant during the evolution, and thus there is no well-accepted definition 
of quantum chaos. To solve this problem, various signatures of quantum chaos have been identified [254, 
256, 257]. Entanglement and spin squeezing [1, 3], which are pure quantum effects, have also been identified 
as signatures of quantum chaos. Recently, entanglement, measured by the linear entropy, as a signature 
of quantum chaos has been demonstrated experimentally in an atomic ensemble [258]. The experimental 
results and the theoretical predictions coincide very well [258]. 

6.1. Spin squeezing and quantum phase transitions in the Lipkin-Meshkov-Glick model 

Below, we discuss the spin squeezing for the ground state of the Lipkin-Meshkov-Glick model [259], which 
occurs a typical second-order QPT [222, 260]. It has been widely studied statistical mechanics of quantum 
spin system [261], Bose-Einstein condensations [262], and superconducting circuits [263]. It is an exactly 
solvable [264, 265] many-body interacting quantum system, as well as one of the simplest to show a quantum 
transition in the strong coupling regime. It is convenient to discuss spin squeezing in this model, since it 
consists of collective spin operators. 

The Hamiltonian of the Lipkin-Meshkov-Glick model reads 

H = -±(j*+~,J*)-hJ„ (291) 

where J a = { a — x,y,z) are the collective spin operators; o~ l a are the Pauli matrices; N is the 

total spin number; 7 € [0, 1] is the anisotropic parameter; A and h are the spin-spin interaction strengths 
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Figure 6: (Color online) Squeezing parameter §^ versus the applied magnetic field h (a), and the 'phase diagram' of spin 
squeezing in the h — f plane (b), for the ground state of the Lipkin-Meshkov-Glick model. The system size used here is N = 2 7 . 
The maximal squeezing occurs at the critical point. The two black dots denote £g = 1, at these points the ground states are 
coherent spin states. In (b), we show the squeezed and non-squeezed regions that are separated by ho = ^y, £g = 1 in the 
thermodynamic limit. 



and the effective external field, respectively. Here we set A = f, where the Lipkin-Meshkov-Glick model 
describes ferromagnetism. The QPT of this model originates from the competition between the spin-spin 
interaction and the external field. With mean-field approach [266]; we can see that, when h > 1, all spins 
tends to be polarized in the field direction ((ct*) = l). However, when h < 1, it is two-fold degenerate with 
(cr*) = h. Therefore, a spontaneous symmetry breaking occurs at h = 1, which is a second-order QPT 
point between the so-called symmetric (h > 1) phase and symmetry broken (h < 1) phase. However, by 
considering the quantum effects, the exact ground state is not degenerate in the symmetry broken phase 
(7^1). Since the Hamiltonian is of spin-flip symmetry, i.e., H, J\iLi o\ = 0, we have (J x ) = (J y ) = 0, 

(JxJz) — (JyJz) — 0, and the MSD is along the z direction. In addition, [if, J 2 ] = 0, and the ground state 
lies in the j = N/2 symmetric sector. 

In the isotropic case, 7=1, the ground state is simply the Dicke state, as discussed in Sec. 2, and 
£l r > 1 for all Dicke state, indicating that there is no spin squeezing with respect to £| and £ R . For the 
Dicke state \j, m) , we have 

£% = ^ < 1, (292) 
3 

thus the squeezing occurs along the z-axis direction with respect to 

In the anisotropic case, 7 7^ 1, there is a second-order QPT. Since the system is of exchange symmetry, 
and the ground state has a fixed parity, £| and are closely related to the concurrence, as demonstrated 
in Refs. [222, 260, 223]. In the thermodynamic limit, we could calculate the squeezing parameter by using 
the Holstein-Primakoff transformation. The results are [240] 

2 _ / y /(h-l)/(h- 7 ), for h>l, , . 

^ - \ V(l-W(l-7), h < 1. [2J3 > 

As shown in Fig. 6, there is no spin squeezing when h is smaller than ho = N ] ^ 1 y/j. At the point h = ho, 
the ground state is a CSS, for which the squeezing parameters are equal to one. 

To demonstrate the QPT, we calculate the values of and in the vicinity of the critical point. 
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Figure 7: Stroboscopic phase-space dynamics of the classical kicked top for k = 3. Two hundred stroboscopic trajectories are 
plotted, each for a duration of 233 kicks. This figure is from Ref. [193]. 



Reference [260] found that at the critical point 

£§ - i/7v-o.33±o.oi j 
dxt 2 s ~ - l/N-°- 33±Q 03 . (294) 

The above results are immediate consequences of the concurrence C, as shown in Ref. [260]. If £| < 1, we 
find £| = 1 — (N — 1)C due to the exchange symmetry. Generalized spin squeezing inequalities can be used 
to detect the ground state entanglement for lattice models [6]. However, as these inequalities are not given 
in terms of parameters like £|, there are no such scaling laws, and we do not focus on them here. 

Besides spin squeezing and concurrence, other measures such as negativity [267], geometric entangle- 
ment [268], and entropy [269] have also been studied. Due to its symmetry, the spin-squeezing parameter 
£| is closely related to the concurrence in this model. At the critical point, the derivative of the parameter 
£| with respect to the external driving field strength tends to diverge. 



6.2. Spin squeezing and quantum chaos in the Quantum kicked-top model 

In this section and the following, we discuss the behaviors of spin squeezing in two typical quantum 
chaotic systems: the quantum kicked-top model and the Dicke model. To facilitate the correspondence 
between the classical and quantum regimes, the initial states are chosen as CSSs. To connect the quantum 
dynamics to classical chaos, we choose the initial states located in both chaotic and regular regions. The 
dynamics of the squeezing parameters distinguish well between regular and chaotic regions, and we conclude 
that spin squeezing could be a good signature for quantum chaos. 

The quantum kicked-top model can be realized in cold atomic ensembles [38, 111], and it is familiar to 
us since it has the form of a one-axis twisting Hamiltonian with a periodic pulse, 

oo 

H=^-J 2 z+P J y Y, Ht-nr), (295) 

n— — oo 

where the second term is the periodic driven pulse (with discrete kicks). The dynamics of this model is 
described by the Floquent operator 

F = exp ( Jj ) exp (-ipJ y ) , (296) 
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where we choose p = it/2 and r = 1 for convenience. The evolution of an initial state \tp (0)) is given by 
\tp (n)) = F n \ip(0)). The spin-squeezing dynamics can be obtained from this state. The chaotic behavior 
of this model has been recently demonstrated experimentally via linear entropy [258], by using the tensor 
part of an effective Hamiltonian for the dispersive atom-field interaction, and just in the same system, spin- 
squeezed states have been created [38, 111]. Therefore, the theoretical proposal of using spin squeezing as a 
signature of quantum chaos could be realized using current techniques. Below, we present numerical results 
for the dynamics of the spin-squeezing parameters. 
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Dynamical evolution of the spin squeezing parameter £| for the initial CSS with 9q 



Figure 

The other parameters are fixed here at re 
is from Ref. [193]. 



2.25 for different values of <f>Q. 
3 and j = 25. The inset displays the initial dynamics when 4>q = 0.63. This figure 



To demonstrate the chaotic properties for this model, one could represent the model Hamiltonian un- 
der the CSS in large- j limit. The classical limit of the quantum kicked top is obtained by expressing 
(X, Y, Z) = ((J x ) , (J y ) , (J z )) /j, and factorizing all the second moments to products of first moments, 
such as (J x Jy) /j 2 = XY. Then the classical equations of motion can be obtained from the Heisenberg oper- 
ator equations of motion with the factorization rule. The classical normalized angular momentum variables 
(X, Y, Z) can be parameterized in spherical coordinates as (X, Y, Z) = (sin 9 cos <fi, sin 9 sin (j), cos 9), where 
and <j) are the polar and azimuthal angles, respectively. Thus, the map is essentially two dimensional. 

The stroboscopic plot with parameter k — 3 is displayed in Fig. 7. The parameter k is chosen such as to 
yield a mixture of regular and chaotic areas. Elliptic fixed points surrounded by the chaotic sea are evident. 
Two such elliptic fixed points have coordinates (9, <f>) — (2.25, —2.5) and (9, <fi) — (2.25, 0.63). As we will 
see, this phase-space structure of the classical kicked top determines spin squeezing in the quantum kicked 
top. If the initial state |0o>^o) is chosen at one of these two points, the classical dynamics is regular. 

Now, we return to the quantum Hamiltonian. The initial state is also a CSS, while the evolution is via 
the Floquent operator. The numerical results for the squeezing parameter £| are shown in Fig. 8. The 
polar angle 9 = 2.25 is chosen for all the initial states. The four plots show very clear variations of the 
dynamics of £| for the fixed point (4>o = 0.63), regular region (<pQ = 0.5), and the chaotic region (<po = 0.0 
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Figure 9: Poincar section for the spin degrees of freedom (section with 52 = 0.0 and p2 > 0.0) in the resonant case (u = dJo = 1), 
energy i? = 8.5. This figure is from Ref. [272]. 

and </>o = —1.0). 

Spin squeezing is drastically suppressed in the chaotic region. As shown in Fig. 8, spin squeezing occurs 
frequently for the initial state located at the elliptic fixed point </>o = 0.63, whereas spin squeezing vanishes 
after a very short time (n = 2) for the initial state is in the 'deep' chaotic region (4>q = —1.0). When the 
initial state is in the regular region ((f> = 0.5), spin squeezing exists over a relatively long time, and only 
when steps n > 318, the spin squeezing vanishes. When the initial state is centered in the 'shallow' chaotic 
sea (f>o = 0, spin squeezing disappears after n = 4 and never revives. The strong spin-squeezing oscillations 
mainly originate from the periodic kicks. Without kicks, the spin squeezing exhibits only periodic regular 
oscillations. Thus, we see that spin squeezing is very sensitive to classical chaos, and classical chaos strongly 
suppresses spin squeezing. 

6.3. Spin squeezing and quantum chaos in the Dicke model 

In this section, we consider the Dicke model, which characterizes the interaction between optical field 
and atoms. The Dicke model describes an ensemble of two-level particles interacting with a single-mode 
cavity bosonic field, and the corresponding Hamiltonian reads (H = 1) 

R R f 

H = uj q J z + ua ] a + —= (,J + a + J_a f ) + -== ( J + a f + J_a) , (297) 

where the collective pseudo-spin operators are J a = X)i=i CT ia/2, (a = x,y,z) with <7, a being the Pauli 
operator for the i-th particle, and R, R' are the coupling strengths. The denominator y/2] arises from 
that the dipole coupling strength is proportional to l/yF, where V is the volume of the cavity. When 
we consider the ground state, the system undergoes a second-order QPT between a normal phase and a 
super-radiant phase [270, 271]. This Hamiltonian is integrable under the rotating- wave approximation, i.e. 
R' = 0, which is used to demonstrate the transfer of squeezing between the field and atoms, while in this 
case no chaos exists. In the non-rotating-wave-approximation regime, R' 7^ 0, its classical dynamics exhibits 
chaos. Therefore, this system is a desirable model for studying the interface between quantum and classical 
chaos. 

As discussed previously, to study the quantum chaos, we write the Hamiltonian in the classical limit 
(in the coherent-state representation). As this system consists of the atomic and optical parts, the initial 
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quantum states are chosen as follows: 



|V(0)> = \v)®\a), 



(298) 



where I77) and \a) are the coherent spin and bosonic states, respectively. Here the variables 77 and a can 
be written as functions of the classical variables in the corresponding phase spaces, (qi,Pi) for the atomic 
degree of freedom, and ((72, Pa) for the bosonic field r\ — (pi + iqi) /\/4j — (jp\ + q\), a = ~^(P2 + ^92), 
where q\ , p\ , qi , pi describe the phase space of the system under consideration, and the indices 1 and 2 
denote the atomic and field subsystem, respectively. The classical Hamiltonian corresponding to Eq. (297) 
can be obtained under the CSS representation. 





Figure 10: (Color online) Dynamical evolution of the spin-squeezing parameter £jj for the initial states with q\ = and 
different p\. (a) The regular region corresponding to Fig. 9: pi = 2.0 (solid line) and pi = —3.5 (dashed line); (b) Chaotic 
region corresponding to Fig. 9: p± = (solid line) and pi = —1.0 (dashed line). This figure is from Ref. [273]. 

The classical dynamics associated with the Dicke Hamiltonian was explored in Ref. [274]. It was shown 
that the integrable cases are recovered when either R or R' is zero, and the most chaotic dynamics is 
associated to the condition R = R', In Fig. 9, we show the Poincar section for the classical counterpart of 
the spin degrees of freedom, defined by the section qi — in the four-dimensional phase space so that every 
time a trajectory pierces this section with p 2 > 0, the corresponding point (qi, pi) is plotted. Here, the 
total energy is fixed at E — 8.5, j — 4.5, u = luq = 1. The limit of atomic phase space is indicated by a 
border of radius equal to \f£j. In Fig. 9, we choose the coupling parameters R = 0.5 and R' = 0.2, which 
yield a mixture of regular and chaotic areas of significant sizes. Many fixed points and near-periodic orbits 
surrounded by the chaotic sea are evident. 

The dynamic behavior of the spin-squeezing parameter f| is computed by choosing the initial states in 
the chaotic and regular regions with respect to the Poincar section. The results are shown in Fig. 10. We 
choose two fixed points at q\ — 0, p\ — 2.0 and q± = 0, p\ = —3.5, which are in the regular regions, and two 
additional points q\ — 0, p\ — and q\ = 0, p% = —1.0, well inside the chaotic sea. 

For the regular regions [see Fig. 10 (a)], spin squeezing vanishes after a relatively longer time (t w 9.0), 
whereas for the chaotic region [see Fig. 10 (b)], the spin squeezing vanishes after a short time (t « 1.8). 
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In summary, spin squeezing is sensitive to and suppressed by chaos. In contrast to this, chaos enhances 
the bipartite entanglement quantified by the linear entropy in the Dicke model [272]. So, the underlying 
chaos affects differently the two typical pure quantum mechanical phenomena, spin squeezing and bipartite 
entanglement. These results are very similar to the case of quantum kicked top model. Numerical results 
for the quantum kicked top and Dicke models suggest that spin squeezing is a good signature of chaos. 
Furthermore, by considering the advantage of the spin-squeezing parameter in measurement, it may have 
good applications in experiments. 



7. Spin squeezing under decoherence 

In this section, we take into account decoherence [275, 276, 277, 278], which is induced by unavoidable 
technical noise, and interactions with the environment. We shall discuss the effects of decoherence in three 
aspects. 

(i) Decoherence in the generation of spin squeezing. In Sec. 7.1 we shall show that, if the decoherence 
effect is sufficiently weak, although the ability of squeezing generation is weakened, considerable large amount 
of squeezing can still be obtained. 

(ii) The robustness and lifetime of spin squeezing in the presence of decoherence. In Sec. 7.2, we show 
that spin-squeezed states are more robust to particle-loss [279] as compared to the GHZ state. In Sec. 7.3, 
the lifetime of spin squeezing is studied in three types of decoherence channels, which are quite general to 
describe decoherence processes, and spin squeezing sudden death is observed. 

(iii) Decoherence in the quantum metrology process. In Sec. 7.4, we discuss how decoherence affects the 
ability of spin-squeezed states to the improvement of phase resolution in Ramsey processes. Although the 
degree of spin squeezing is weakened by decoherence, sub-shot noise phase resolution can still be achieved [13], 
in contrast with the maximal entangled state, with which we could only attain the shot-noise level precision 
in the presence of decoherence [280] . 

7.1. Decoherence and spin squeezing generation 

In Sec. 3, we have discussed generating spin-squeezed states via ideal one-axis twisting Hamiltonian, 
where decoherence effects are not taken into account, and the best squeezing scales as £jj oc 1/iV 2 / 3 for 
the one-axis twisting and £^ oc 1/N 2 for the two-axis twisting. It is natural to expect that, decoherence is 
detrimental to spin squeezing generation, while if the decoherence effect is sufficiently weak, considerable 
large amount of squeezing can still be obtained. Below, we mainly consider the particle loss induced deco- 
herence, which is an unavoidable source of decoherence in cold atom systems, such as BEC, due to collisions 
of condensed atoms with the background gas. 

We first consider the one-axis twisting Hamiltonian. In Ref. [4], to estimate the effect of particle losses 
in a two- mode BEC, they have performed a Monte Carlo simulation of the evolution of squeezing from the 
one-axis twisting Hamiltonian H = xJ 2 . The particle loss is phenomenologically taken into account by 
introducing a loss rate T. Their simulations indicated that even under the conditions that approximately 
10% atoms are lost, squeezing of £^ ~ 0-01 ma y be obtained. In Rcfs. [153, 154], the authors presented a 
more detailed analysis of spin squeezing generation in the presence of particle loss in a two-mode BEC. The 
time evolution of the squeezing parameter is derived as 
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where £fj (t) is the squeezing without particle loss, and 

r^E 1 ^* and r ?= mrW . ( 30 °) 

rn 

where r( m ) is the m-body loss rate. They then found that, for one-body losses, the best obtainable spin 
squeezing scales as 

oc 7V~ 4 / 15 for N — y oc, (301) 
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for two-body loss, in the limit N — > oo the best squeezing is independent of N, and for three-body losses, 
the best squeezing 

f R oc N 4/15 for N -> oo. (302) 

Thus for one-body losses, we can obtain arbitrary small squeezing parameter by increasing the particle 
number, while for three-body losses case, there is a finite optimum number of particles for squeezing. For 
two-axis twisting model, the decoherence effects induced by particle loss was studied in Ref. [212], where 
they established an Raman-scattering based approach to implement the two-axis twisting Hamiltonian as 
Eq. (170). They found that, to achieve the Heisenberg limit squeezing, it is required that \ > T, where 
X is the twisting strength in Eq. (170), and T is the one-body loss rate. To achieve any squeezing, it 
requires N\ > I\ where N is the particle number. Quantum-limited metrology in the presence of collisional 
dephasing was studied in [281, 242], although these works are associated with nonlinear estimation scheme, 
their analysis can also be extended to the discussion of spin squeezing. Otherwise, a detailed analysis of 
decoherence in quantum light-atom interfaces, such as QND-type interactions, was studied in [282]. 

7.2. Spin-squeezed states under particle loss 

In the above subsection, we discussed how decoherence affects the production of spin squeezing. Below, 
we consider the effects of particle loss on the generated spin-squeezed states, which was studied in Refs. [279, 
126, 153]. For convenience, the system is now assumed to be symmetric under particle exchange. Based on 
this assumption, the spin squeezing after particle loss is equivalent to the spin squeezing of the subsystem 
composed of the remaining particles. The state after particle loss is obtained by tracing over the degrees of 
freedom of the lost particles, and this procedure is what we used when deriving the reduced density matrix 
of the state for the subsystem. Therefore, we immediately arrive at the conclusion that concurrence is not 
affected by particle loss. 

Consider an iV-body system with exchange symmetry, represented by a density matrix . The reduced 
states also have exchange symmetry. After particle loss, the reduced state for the remaining N r particles 
becomes 

P Nr = Ti(N-N r ){p N }, (303) 
where we trace over the (N — N r ) lost particles. As the state is of exchange symmetry, we have 

(o~lct) = (o~ia) n = (®ia) N r ' 
(o~l a C2a) = {o'io.O j a ) N — {o'iaO ja) Nr j (304) 

where the subscripts N and N r denote the state before and after particle loss, respectively. The indices i, j 
denote any two remaining particles. Note that, according to Eq. (304), the MSD is invariant under particle 
loss. 

From the expression J 7 i = | 53i=i a ini one can verify that 

N N 2 — N 

(4)n = T + 4 fogo-aa), (305) 
where we have used the symmetry of the state. Conversely, one finds 

(ffi»cr M ) = ^ ( - l) . (306) 

For the state after loss, similarly, we obtain 



N- 1 V N 



<^> = ^ - 1 ) ■ (307) 



Thus, from Eqs. (306) and (307), we obtain the following relation 

H4)N r N r -14(4) N N — N r 



N r N - 1 N N-l 
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(308) 



which is applicable to an arbitrary direction n. Then, we have a similar expression for the direction n± 
(perpendicular to the invariant mean-spin direction), 



H4jNr _ N r -lH4 ± ) N N-N r 

N r N-l N + ^V^T' (309) 

The two quantities 4:(J^ ± )N r /N r and 4(J|^)tv/A~ are just linearly related to each other, and N r — 1 > 0. 
So, if we choose ri± be the direction where the minimal variance is obtained for the state before loss, the 
direction n± also corresponds to the minimal variance after loss. Finally, from the definition of the squeezing 
parameter £| (Table 1), we obtain [279] 

which gives the relation of the spin-squeezing parameters before and after loss. 

Now, we discuss the spin-squeezing parameter which is simply related to the parameter £| via 

N 2 1 
= = W^s,N- (311) 



A similar expression holds for the squeezing parameters after loss, and it is given by 

1 

Multiplying Eq. (310) by (<7i)~ 2 , one immediately obtains [279] 

C 2 _ N r - 1 2 N - N r 1 

6?jv„ - - n -=y£r,n+ n _ 1 {gi)2 - (313) 

The squeezing parameter ^ N only requires the knowledge of £| N , while the parameter £|j N is related to 
both £|j N and (<?i) 2 . 

Note that the GHZ state, which is maximally entangled, is extremely fragile under particle loss. Indeed, 
the entanglement is totally destroyed when only one particle is lost. However, the spin-squeezed states, 
which are only partially correlated (e.g., pairwise correlated) are more robust. We emphasize that the 
robustness is different in the depolarizing channel. It has been found that in the limit N — > oo, GHZ states 
can stand more than 55% local depolarization, while for spin-squeezed states it can stand more than 29% 
local depolarization [283]. 



7.3. Spin squeezing under decoherence channels 

In the previous subsection, the decoherence is induced by particle loss. Here, the decoherence effects 
are described by three types of decoherence channels: the amplitude damping channel (ADC), the phase 
damping channel (PDC), and the depolarizing channel (DPC). In general, decoherence processes can be 
described by these three typical channels. They are prototype models of decoherence relevant to various 
experimental systems. They provide "a revealing caricature of decoherence in realistic physical situations, 
with all inessential mathematical details stripped away" [284]. But yet this "caricature of decoherence" 
leads to theoretical predictions being often in good agreement with experimental data. Examples include 
multiphoton systems, ion traps, atomic ensembles, or a solid-state spin systems such as quantum dots or 
nitrogen-vacancy centers, where qubits are encoded in electron or nuclear spins. Let us first introduce the 
three decoherence channels. 
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7.3.1. Decoherence channels 

Amplitude- damping channel. — The ADC is defined by 



£adc(p) =e oP eI+e iP e\ 



(314) 



where E = y/s\0)(0\ + |1)(1| and E\ — 1) (0| are the Kraus operators, p = 1 — s, s = cxp(— 7i/2), and 
7 is the damping rate. In the Bloch representation, the ADC squeezes the Bloch sphere into an ellipsoid 
and shifts it toward the north pole. The radius in the x-y plane is reduced by a factor y/s, while in the 
^-direction it is rescaled by a factor s. The ADC is a prototype model of a dissipative interaction between 
a qubit and its environment. For example, the ADC model can be applied to describe the spontaneous 
emission of a photon by a two-level system into an environment of photon or phonon modes at zero (or very 
low) temperature in (usually) the weak Born-Markov approximation. 

Phase-damping channel. — The PDC is a prototype model of dephasing or pure decoherence, i.e., loss of 
coherence of a two-level state without any loss of the system's energy. The PDC is described by the map 



and obviously the three Kraus operators are given by Eq = \/sI, E\ = v / p|0)(0| and E2 = v / p|l)(l|, where 
I is the identity operator. For the PDC, there is no energy change and a loss of decoherence occurs with 
probability p. As a result of the action of the PDC, the Bloch sphere is compressed by a factor (1 — 2p) in 
the x-y plane. It is evident that the action of the PDC is nondissipative. This means that, in the standard 
computational basis |0) and |1), the diagonal elements of the density matrix p remain unchanged, while 
the off-diagonal elements are suppressed. Moreover, the qubit states |0) and |1) are also unchanged under 
the action of the PDC, although any superposition of them (i.e., any point in the Bloch sphere, except the 
poles) becomes entangled with the environment. The PDC can be interpreted as elastic scattering between 
a (two-level) system and a reservoir. The PDC is also a suitable model to describe T 2 relaxation in spin 
resonance. 

Depolarizing channel. — The definition of the DPC is given by the map 



We see that for the DPC, the spin is unchanged with probability ,s = 1— p, or is depolarized to the maximally- 
mixed state 1/2 with probability p. It is seen that due to the action of the DPC, the radius of the Bloch 
sphere is reduced by a factor s, but its shape remains unchanged. 



£pdc(p) = s/> + p(Ax>|0)<0|+pii|1><1|) 



(315) 



£dpc(p) = sp + p-. 



(316) 
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Table 7: Analytical results for the time-evolutions of all relevant expectations, correlations, spin-squeezing parameters, and 
concurrence, as well as the critical values p c of the decoherence strength p [197]. This is done for the three decoherence channels 
considered here. For the concurrence C, we give the expression for C' r , which is related to the rescaled concurrence C r via 

C r = max(0, C' r ). The parameters xq = 1 + 2 (<T Z ) + (ci z <T2 z )q, ao = (N — 1) (l — (criztTzz) o) • The critical values p^' , p^?\ 

and pi 3 ' correspond to the concurrence, squeezing parameter and respectively. See Ref. [197] 
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7.3.2. Spin squeezing parameters under decoherence 

Now we begin to study spin squeezing under the above three different decoherence channels. The spin- 
squeezed states are prepared by a one-axis twisting Hamiltonian, and the spin-squeezing parameters £|, £|, 
and £fj can be evaluated by using Eqs. (104) and (111), since the system has exchange symmetry, and the 
decoherence channels act independently on each particle. Therefore, all the spin-squeezing parameters and 
the concurrence are determined by some correlation functions and expectations, and they can be calculated 
via the Heisenberg approach shown below. 

The spin-squeezing parameters and concurrence for the initial state are shown in Table 7. The state 
under decoherence is given by 

£{p)= E (®£i*PiM®£i4i) > ( 31? ) 
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where Kn t denotes the Kraus operator for the i-th particle. The spin-squeezing parameters and concurrence 
consist of some correlation functions and expectations, and we need just to calculate the evolution of these 
quantities by employing Heisenberg picture without writing the state explicitly. The evolution of operator 
A can be put forward by employing the Heisenberg picture with 

(A) = Tx[A£(p)]=Tr[£HA)p\, (318) 

where 

£HA)= £ {®tiKl)A(®tiK^) (319) 
is relatively easy to calculate for our cases. 
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Figure 11: (Color online) Spin-squeezing parameters (red curve with squares), f|, (top green curve with circles), and the 
concurrence C r (solid curve) versus the decoherence strength p = 1 — exp(— jt) for the amplitude damping channel, where 7 is 
the damping rate. Here, 80 is the initial twist angle given by Eq. (323). In all figures, we consider an ensemble of N = 12 spins. 
Note that for small initial twist angle 80 (e.g., 60 = D.l7r), the two squeezing parameters and the concurrence all concur. For 
larger values of do, then Q%, fjj,, and C become quite different, and all vanish for sufficiently large values of the decoherence 
strength. This figure is from Ref. [197]. 

In order to characterize spin squeezing more conveniently, we now define the following squeezing param- 
eters: 

C, 2 = max(0,l-^), ke{S,R,E}. (320) 

Spin squeezing appears when > 0, and there is no squeezing when vanishes. Thus, the definition of 
the first parameter Qg has a clear meaning: namely, it is the strength of the negative correlations as seen 
from Eq. (179). More explicitly, for the initial state, we have £l = 1 — [N — 1)Cq as shown in Eq. (132), so 
d is just the rescaled concurrence 

C r (0) = (N - 1)C . (321) 
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To study the decoherence of spin squeezing, we choose a state which is initially squeezed. One typical 
class of such spin-squeezed states is the one-axis twisting collective spin state [1] , 

|tf(0 o )) = e - ie ° J * /2 \l}® N , (322) 
which could be prepared by the one-axis twisting Hamiltonian H = xJj, where 

6 = 2 X t (323) 

is the one-axis twist angle and x is the coupling constant. 

The numerical results for the squeezing parameters and concurrence are illustrated in Fig. 11 for different 
initial values of 6. For the smaller value of 8 , e.g., #0 = tv/10, we see that there is no entanglement sudden 
death (ESD) and spin squeezing sudden death (SSSD). The spin squeezing and the pairwise entanglement 
are completely robust against decoherence. Intuitively, the larger is the squeezing, the larger is the vanishing 
time for the squeezing. However, here, in contrast to this, no matter how small the squeezing parameters 
and concurrence are, they vanish only in the asymptotic limit. This results from the complex correlations 
in the initial state and the special characteristics of the ADC. For larger values of 9q, as the decoherence 
strength p increases, the spin squeezing decreases until it suddenly vanishes, so the phenomenon of SSSD 
occurs. There exists a critical value p c , after which there is no spin squeezing. The vanishing time of £|, 
is always larger than those of £jj and the concurrence. We note that depending on the initial state, the 
concurrence can vanish before or after This means that in our model, the parameter ^ 2 E < 1 implies the 
existence of pairwise entanglement, while £jj does not. 

We can calculate all the relevant correlation functions, squeezing parameters, concurrence and the critical 
values of the decoherence strength p, which are given in Table 7. Initially, the state is spin-squeezed, i.e., 
£l(0) < 1 or C r (0) > 0. As seen from the Table, one can find that £| < 1, except in the asymptotic limit 
of p — 1. Thus, we conclude that according to £|, the initially spin-squeezed state is always squeezed for 
p ^ 1, irrespective of both the decoherence strength and decoherence models. In other words, there exists 
no SSSD if we quantify spin squeezing by the first parameter 




Figure 12: (Color online) Critical values of the decoherence strength p^ (blue solid curve), p^ (red curve with squares), p£ 3 ' 
(top green curve with circles), and the squeezing parameter £g (black dashed curve) versus the initial twist angle 6>o given by 
Eq. (323) for the amplitude-damping channel, PDC. Here, p c is related to the vanishing time t v via p c = 1 — exp(— ytv)- At 
vanishing times, SSSD occurs. This figure is from Ref. [197]. 
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The SSSD for the three decoherence channels with respect to different initial twisting angles was also 
studied [197], here, we only take the PDC for example. In Fig. 12, we plot the decoherence strength p c 
versus the twist angle 9q of the initial state for the PDC. For this decoherence channel, the critical values 
Pc's first decrease, until they reach zero. Also, it is symmetric with respect to 9q = it, which is in contrast 
to the ADC. There are also intersections between the concurrence and parameter and the critical value 

Pc 3 '' is always larger than pi 1 ' and p c 2 ^ . 

As summarized in Ref. [197], the common features of the three decoherence channels are: (i) The critical 

(3) 

value p c is always larger than or equal to the other two, namely, the spin-squeezing correlations according 
to £jj are more robust; (ii) There always exist two intersections between the concurrence and the parameter 
for 6 from to 2tt, irrespective of the decoherence channels; (iii) When there is no squeezing (central 
area of Fig. 12), all vanishing times are zero. Table 7 conveniently lists all the analytical results obtained in 
this section. 



7.4- Effects of decoherence on spin squeezing in Ramsey processes 

In this section, we present the effects of decoherence on spin-squeezed states in Ramsey intcrfcrometry, 
which was discussed in Sec. 5.3. Now we take into account decoherence, and each spin-1/2 particle is assumed 
to evolve independently via the master equation 



(324) 



where 7 is the decoherence rate. Under the above master equation, which describes a dephasing process, 
nondiagonal elements of p decay as exp (— ft), while the diagonal terms are invariant. Spin squeezing 
weakened by dissipation effects was studied in Refs. [62, 285]. By using the results obtain in Table 7, we 
now derive 



(Jz) t =-cos4> (Jz) t =o exp(-Tt), 



; <t> ( J x) t=0 + cos 2 </> (AJ Z ) 2 =0 - sin(2</>)Cov( J x , J z ) t=0 exp{-2ft) 



N 

-[I- exp(-2 7 t)] 



(325) 



where the exponential terms come from the free evolution, when the system is rotating around the z-axis. 
Thus, from the above equations, the decoherence effects on the spin-squeezing parameter £^ are clear: the 
length of the spin becomes shorter, while the fluctuation, which maybe squeezed at the beginning, becomes 
larger. Therefore, the degree of squeezing becomes smaller. 

Since <f> is time dependent, it is more convenient to estimate u)q. By using Eqs. (273) and (325), we 
obtain 



Aw = 



m 2 ) 



cos' 



•(Aj z r 



[exp(2 7 i) - 1] 



-1 1/2 



tT sin 2 (j>(J z 



(326) 



where we omit the subscript t = 0, and assume Cov ( J x , J z ) t=0 = 0, which is usually satisfied when the 
initial states have parity, e.g. the spin-squeezed states generated via the twisting Hamiltonian. 

Now, we search for the point (<^ pt:^opt) where the minimum of Auj is attained. The minimization of 
Awq is equivalent to the minimization of (Awq) 2 , which is given by 



(Awo) 2 - { ''' } - (A f + (AA)i .. 

tT(J z ) 2 tT sin 2 4>(J Z 



f [exp(2 7 i) 



(327) 



It is evident that (Aw ) 2 decreases monotonically with the increase of sin</>, thus 



,J opt 



{k odd) . 



(328) 
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Figure 13: (Color online) Schematic diagram of the state evolution for Ramsey spectroscopy. Here the initial state is spin- 
squeezed in the rr-direction, with spin length J. At each "snapshot", the initial and final states are joined by a white "path" or 
trajectory on the Bloch sphere. The phase <f> gained in the free-evolution period is determined by measuring the population 
differences, between particles in states |0) and |1), which is equivalent to measuring J z . The precision of is thus determined 
by the spin length and the fluctuation in J z . In the absence of decoherence (a), a squeezed state (fuzzy elliptic rainbow) 
performs better (i.e., has smaller fluctuations in J z ) than a coherent spin state (blue circle), shown in a3. When decoherence 
is considered, during the free-evolution period, the spin length is shorter and the fluctuations become larger. Thus, in the final 
population measurement, b3, the precision of <j> becomes lower. 



GO 



And at this point, (Aloq) becomes 



(Aw ) (0 = (popt) 



(J*) + £ [exp(2 7 i) - 1)] 



tT (J, 



(329) 



It is interesting to see that this minimal value is independent of the variance (AJ Z ) , irrespective of deco- 
herence. Now, we minimize (Aujq) with respect to t. The value of t opt is determined by the conditions 



0_ 
Of 

dt 2 



(Aw ) {<t> = 0opt) 

(Aw ) 2 {4> = 4>o V t) 



= 0, 

> o, 



and thus t opt is the solution of 



where 



(2 7 <-l)exp(2 7 + l = e i 2 



(330) 
(331) 



4 a 2 ) 



N 



(332) 



which is equal to £| if (>/ 2 ) is the minimal spin variance in the n± direction. When t = t opt , one finds 



dt 2 



(Acj ) (</> = <^o P t) 



iV 7 2 exp(2 7 i) 
Tt(J z ) 2 : 



(333) 



which is obviously larger than zero. Thus, at the optimal point, we obtain the minimal value. 
From Eqs. (329) and (331), we obtain the minimal Awq as 



min (Awq) = 



'2 7 exp(2 7 t opt ) 



TN Vz 



(334) 



where T] z — 4 (J z ) 2 /N 2 , which is always less than or equal to one. For the CSS \j,—j), we have rj z = 1, 
^ = 1 and t opt = 1/ (2 7 ). Thus according to Eq. (334), we have 



min (Awq) 



2 7 exp (1) 



A relative improvement in the precision over the above minimal uncertainty may be given by 



P(tapt,Vz) = 1 - 

Then in the idealized case, t opt — > and ry z - 
improvement is 



min (Awo) 
min (A wq' 



= 1 - 



CSS 



/cxp (2 7 t opt ) 
exp(l)% 



(335) 



(336) 



1, the theoretical absolute decoherence limit in the precision 

(337) 



1 — exp 



0.39. 



which is independent of N. Here, we emphasize that, the limit t — > implies £| — ► according to Eq. (331). 

In Ref. [280], a similar issue was studied with the maximally entangled state, the GHZ state, which is not 
spin-squeezed but can be used to achieve the Heisenberg limit [33, 37, 252]. Compared with spin-squeezed 
states, for Markovian dephasing processes as discussed in this section, the best precision given by the GHZ 
state is the shot-noise limit, the same as a CSS. Recently, the lifetime of spin-squeezed states under the 
influence of dephasing noise was measured in Ref. [44], where the spin-squeezed state was generated in a 
collection of atoms in a cavity via feedback mechanism. For experimental conditions in [44], the lifetime is 
r ~ 600 )j,s for an initial spin-squeezed state with ~ 4 dB. Since the CSS is also affected by the same 

dephasing noise, at the time, r = 600 /is, even though the initially squeezed states becomes unsqueezed, 
£p = 1, it still improves the signal-to-noise ratio by ~ 3dB over that of an initially CSS. 
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8. Theoretical proposals and experimental realizations of squeezing generation 

Above we have discussed many theoretical works about spin squeezing and its applications. Below, 
we review the generation of squeezing in several physical systems, with both theoretical proposals and 
experimental realizations. In Sec. 8.1, we review experimental progresses of generating spin squeezing in 
the BEC, for which the theoretical proposals have been discussed in Sec. 3.1.3. Then in Sec. 8.2, we discuss 
the transfer of squeezing from squeezed lights to atoms. Afterwards, spin squeezing produced by QND 
measurement is reviewed in Sec. 8.3. 



8.1. Generating spin squeezing in Bose-Einstein condensations 

In this section, we discuss how to generate spin-squeezed states in BEC. There are two main advantages 
for producing spin-squeezed states in this system: the considerable long coherence times, and the strong 
atom-atom interactions, which induce nonlinearity and squeezing. The simplest and most widely studied 
scheme to generate spin-squeezed states in a two-component BEC is to utilize particle collisions. Under the 
single-mode approximation, the BEC Hamiltonian can be mapped to the one-axis twisting Hamiltonian in 
the Schwinger representation [4, 17, 119, 120, 121, 122, 123, 124, 125, 126, 127, 129, 130, 131, 132, 133, 18, 
134, 135, 136], which was realized in experiments [17, 18, 19, 20]. The validity of the two-mode approximation 
has been confirmed by using perturbation-based Bogoliubov theory [125], and also by employing the positive- 
P method [120], both gave a more detailed analysis than using a single-mode approximation. References 
[137, 138] presented an effective two-axis twisting Hamiltonian resulting from Raman scattering. Spin 
squeezing in a spinor-1 BEC was studied in Refs. [145, 144, 146, 147, 142, 148, 143]. There, spin-squeezed 
states could be generated via spin-exchange interactions [145, 144]. For a spinor-1 BEC, which should be 
described by the SU(3) algebra, spin-squeezed states can be created in an SU(2) subspace [146, 147, 143], and 
even simultaneously in all the three SU(2) subspaces [147]. Decoherence in the generation of spin squeezing 
was studied in the cases of particle collisions [151, 141], and particle losses [126, 153, 154]. Below, we present 
some recent experiments about producing spin-squeezed states. 

8.1.1. Experimental realizations 

In the BEC system, we use the concept of spin squeezing since physical observables associated with the 
TV two-level atoms can be described by a fictitious spin J = N /2 in the Schwinger representation, shown in 
Eq. (151). Below, we show some very recent experimental results [19, 20]. To illustrate these results, we 
explain the so-called number-squeezing parameter, which is defined as 

& - ^ (338) 

It is called number-squeezing parameter because J z measures the number difference between the two com- 
ponents, it is related to the metrological squeezing parameter as 

& = (339) 

where 

C-3$l, (340) 

is the contrast of the Ramsey interferometer. 

The intrinsic physics of the generation of squeezed states in BEC is the particle collisions, as illustrated in 
Sec. 3, where we employed a simple two-component BEC, for which an effective one-axis twisting Hamiltonian 
is derived. In experiments, the BEC can be loaded in a optical lattice, and is described by a Bose-Hubbard 
model, of which the atom number statistics undergoes a drastic change from a superfluid phase to a Mott 
insulator phase [149], where the atom number fluctuation is strongly suppressed. The atomic number 
squeezing has been implemented in Ref. [17]. In this experiment, the condensates were loaded in an optical 
lattice, and the system can be described by a Bose-Hubbard Hamiltonian, of which the atom number 



62 



statistics undergoes a drastic change from a superfiuid phase to a Mott insulator phase [149] . For N atoms 
loaded in an M-site lattice, The Bose-Hubbard Hamiltonian reads 



H = 7 53 aj&j +J2 e * 
{i,j) 1 



9_l 
2 



53 hi (hi - 1) , 



(341) 



where a, and dj are bosonic operators of atoms on the ith site, hi = ajd; is the corresponding number 
operator, and q denotes the external potential term. The interaction strength g is associated with the 
s-wave scattering length, and is shown in Eq. (147), and 



(342) 
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The tunneling amplitude 7 between the i and j sites is 

h 2 



7 = / d r 4> (f ~ r i) 



2m 



V 2 + U (r) 



<f>(r-rj). 



(343) 



In the strong tunneling regime gP/'y <C 1, i.e., the superfiuid phase, the interaction term is negligible and 
the single atomic state will spread the lattice. The many-body ground state is given as [286] 



1^)9/3=0 oc ^53 & lj 



N 



|o>, 



of which the on-site number fluctuation obeys Poissonian distribution, i.e., 

An,- = (hi). 



(344) 



(345) 



Conversely, if the ratio g/3/j overcomes a critical point, the ultracold gas enters the Mott insulator phase, 
where the atomic states are localized and the many-body ground state wavefunction becomes 



M n 

|-0mi) 7 =o oc Yl (af) |0), 



(346) 



i=l 



which is a product of local Fock states for each lattice site, with rii the on-site atom number, and the 
corresponding number fluctuations vanish. 

In Ref. [17], atom number squeezing was observed in the Mott insulator phase. They loaded about 
10 s 87 Rb atoms into a time-orbiting potential (TOP) trap and produced about 10 4 condensates of F = 2, 
rrif — 2 atoms. After an adiabatic relaxation, the condensate is then loaded in a one-dimensional vertically 
oriented optical lattice, which consists of ~12 weakly linked mesoscopic wells. The squeezing of the particle 
number was indirectly observed through detecting the phase-sensitivity interference of atoms released from 
the lattice. The phase variance is related to the number variance a n at each lattice as 



1 



(347) 



due to the phase number uncertainty relation. The increase of the phase variance indicates the reduction 
of the number variance. As discussed in Sec. 3.2.3, the spin squeezing can be produced in BEC due to 
nonlinear particle collisions, while in this experiment, besides particle collisions in the same well, there 
exists a tunneling between neighboring wells. The physics of the formation of number squeezing can be 
illustrated by considering a simple two- well version of the Bose-Hubbard Hamiltonian (341), written as 



H = 7 (a[«R + 4 ffl L) + ^ ( a L flL ) + ( a R a R.) 



(348) 
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where the subscripts L and R represent the left and right well, respectively. In the Schwinger representation 
(151) the above Hamiltonian becomes 

H = jJ x +xJ 2 z +^ r , (349) 

where N is the total particle number, and the squeezing strength \ — g/3. To create number-squeezed states, 
they increase the ratio Nx/"f by varying the intensity of the lattice laser and the initial condensate density. 
The atomic interference is realized by releasing the atoms to fall under the gravity. The contrast of the 
resulting interference pattern is used as a phase probe to distinguish between the limits of coherent and 
Fock states. With the adiabatic increase of the lattice potential depth, they found the increase of the phase 
variance, characterized by the squeezing factor S, defined by 

4 = Sa%, (350) 

where <t? corresponds to non-squeezed case. In the case 

NgP/ 1 = 10 5 , (351) 

where N refers to the particle number in the central well, they found that the squeezing factor 

S dB = 10 log 10 -4- ~ 25 dB, (352) 

corresponding to a number variance a n ~ 1 in the central well. 

In Ref. [149], by using a probe that is sensitive only to the presence of atom pairs at a given lattice, 
they observed a continuous suppression of number fluctuations when the ultracold sample evolves from the 
superfluid regime to the Mott insulator regime. In this experiment, the specific two-particle probe are spin- 
changing collisions, which convert spin / = 1 atoms in the m = Zccman sublevel to pairs with one atom 
in m = +1 and the other in the m = — 1. The advantage of the spin-collisional based probes is that they are 
nondestructive, and they can be resonantly controlled using the differential shift between Zeeman sublevels 
induced by an off-resonant microwave held. The experiments confirmed that the atom number exhibits near- 
Poissonian fluctuations for shallow lattices, and strongly suppressed fluctuations for deep lattice. Besides, 
their results indicate that number squeezing is robust with respect to experimental manipulations. 

Experiments of Ref. [160] presented extended coherence times by a factor of 2 over those expected with 
coherent state BEC interferometry. The coherent time, which is limited due to the mean-field interaction 
induced decoherence, is probed by using the decay of Bloch oscillations. The theoretical treatment for a 
BEC in an optical lattice begins with the Bose-Hubbard Hamiltonian (341). The tunneling is weakened 
by applying a large energy gradient across the array, which benefits from the advantage that the on-site 
mean-field energy is unaffected. They have demonstrated that, for number squeezed states prepared at 
potential depth U = 22.5Er, where Er = h 2 k 2 /2m and k = 2tt/X, with A = 852 nm, the coherence time 
r c = 19.3 ± 3.5 ms, which represents an increase of a factor of 2.1 over the expected decay time of an array 
of coherent states in the same lattice potential. It is interesting to note that here squeezing extends the 
coherence time: typically, the enhanced fragility of squeezed states to loss mechanisms results in reduced 
coherence times. Then in Ref. [159], phase coherence of two BEC confined in a double- well potential on an 
atomic chip was observed. It holds for times up to ^200 ms after splitting, a factor of 10 longer than the 
phase diffusion time expected for a coherent state for the experimental conditions. 

In the previous experiments [17, 286, 149, 159, 160], the suppression of atom number fluctuations in a 
BEC was observed indirectly. The direct determination of squeezing parameter requires access to the on-site 
atom occupations. This was demonstrated in Ref. [18] by imaging the condensate with a resolution of 1 fim 
(full- width at half- maximum) , which is well below the lattice spacing of 5.7 /im, and fulfils the criterion of 
local measurement. The condensate was loaded in an optical lattice of two to six wells, and the occupation 
number per site ranges from 100 to 1,100 atoms. The wells of the lattice are fully resolved, and thus the 
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Figure 14: (Color online) Comparison of linear and nonlinear interfcrometry [19]. a, Classical Mach-Zehnder interferometer. 
The two-mode field is mapped to a Bloch-sphere under the Schwinger representation (151). Here the input state is equivalent 
to a CSS. b, By using a nonlinear beam splitter, which acts as one-axis twisting Hamiltonian and generates spin-squeezed state, 
the precision of interfcrometry was improved, c, Six independent BECs of 87 Rb were prepared in a one-dimensional optical 
lattice [19]. State-selective time-delayed imaging causes the atomic clouds to have different shapes. For this trap geometry, the 
dynamics is described by a one-axis twisting Hamiltonian. This figure is from Ref. [19]. 



atom number in each lattice site can be determined by direct integration of the atomic density obtained 
by absorption imaging. The observed spin-squeezed states allow a precision gain of 2.3 dB (two wells) and 
3.8 dB (six wells) over the standard quantum limit for interferometric measurements. 

Recently, some notable improvements have been reported [19, 20]. Although squeezed states have been 
generated previously, Ref. [19] first performed a direct experimental demonstration of interferometric phase 
precision beyond the standard quantum limit in a novel nonlinear Ramsey interferometer. The experimental 
scheme, shown in Fig. 14, involved six independent BECs of 87 Rb (total number 2,300) prepared in a one- 
dimensional lattice, enables parallel experiments performed on six wells independently, which results in 
increased statistics for a given measurement time. The hyperfine states \a) = \F = l,mp = 1) and 
1 6) = \F = 2,mp = —1) form a two-level system. The single traps are almost spherical and have dipolc 
frequencies of w trap = 2tt x 425 Hz. For this trap geometry, the single mode approximation is well justified, 
and the dynamics is governed by a one-axis twisting Hamiltonian 

H/h = Au J z + X J 2 Z + VJtp, (353) 

where the third term describes spin rotations around an axis 

J</> — Jx cos 4> + Jy sin <j>, (354) 

due to the coupling of \a) and \b) by using two-photon combined microwave and radio- frequency pulses 
with a Rabi frequency of Q and phase <j). In the experiment, the Rabi frequency Vt can be switched rapidly 
between and 2tt x 600 Hz, allowing for fast diabatic coupling of the states. 

Then we consider the nonlinear interaction strength x- When the two states |a) and \b) are complete 
overlap spatially, according to the discussion in Sec. 3.2.3 we have 

X oc a aa + a bb - 2a ab , (355) 

where the s-wave scattering lengths in this experiment [19] satisfy 

a aa ■ abb ■ a ab = 100 : 97.7 : 95, (356) 
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Figure 15: (Color online) Characterization of the quantum state within a nonlinear interferometer [19]. a, In the left, the 
twisting effects of the one-axis twisting Hamiltonian. In the right, the blue data have been corrected for photon shot noise and 
the red data additionally take the technical noise into account [19]. The best number-squeezing factor was £L = — 8.2^'2 dB. 
b, Noise tomography of the output state of the nonlinear beam splitter. The dotted box indicates the region detailed in a. 
The largest fluctuations measured have a number-squeezing factor of £^ max = +10.3^q'| dB. This figure is from Ref. [19]. 



which results in very small nonlinear interaction strength, \ ~ 0. To overcome this problem, they used 
the Feshbach resonance to control scattering lengths. By choosing the magnetic field B = 9.10 G, a a b can 
be decreased by a narrow Feshbach resonance and they achieved x — 27r x 0.063 Hz. The magnetic field 
is preserved constant throughout the whole interferometric sequence, while during the coupling pulses the 
interaction strength \ is negligible since the system is in the Rabi regime, that is, -C 1. 

The experiment begins with a CSS that polarized to the —z direction according to the chosen coordinate 
system. After a ir/2 pulse the state turns to be polarized in the x direction, thus (J z ) = (J y ) = 0. Next, 
by means of Feshbach resonance, the state evolves under the one-axis twisting Hamiltonian, which causes a 
shearing effect, shown in Fig. 2, and finally the state has a squeezing angle ao with respect to the z direction. 
As shown in Fig. 14, to make the state useful for squeezing metrology, it is first rotated around the center of 
the uncertainty ellipse by a = ao + 7r/2, and then evolves freely for r = 2 /is, during which a phase difference 
ip is gained between the two states \a) and \b). Followed by a it/ 2 pulse, this phase difference is estimated 
by a population measurement. This standard Ramsey process is illustrated in Fig. 14. The experimental 
results for number squeezing are shown in Fig. 15. The minimal fluctuations are found for rotation angle 
a = 16°, where the corresponding number squeezing factor is 

^ = 4(A^ = _ 69± o.8 dBj (35?) 

for which the photon shot noise due to the imaging process is removed. By additionally taking into account 
the technical noise due to coupling-pulse imperfections, the squeezing factor can be improved to 

& = -8.2t°JdB, (358) 

which is close to the atom-loss-limited theoretical optimum for this system [154]. Furthermore, they pointed 
out that the best spin squeezing factor = cos ip = —8.2 dB is obtainable, which implies the entangle- 
ment of 170 atoms [227]. 

In another recent work [20], spin squeezing was generated in an atom chip by controlling the elastic 
collision interactions through the wavefunction overlap of the two states with a state-dependent microwave 
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Figure 16: (Color online) Spin squeezing and entanglement through controlled interactions on an atom chip [20]. a, Bloch 
sphere representations of the dynamic evolution of the BEC internal states. The red disks represent spin noise. The one-axis 
twisting Hamiltonian x^l reduces the spin noise at an angle 6 m i n , and the subsequent pulse rotates the state around —x by 
an angle 9 [20]. b, Control of the nonlinearity \ on the atomic chip. \ depends on the difference of intra- and inter-state 
atomic interactions, c, Experimental sequence and motion of the two BEC components [20]. d, Measured Ramsey fringes in 
the normalized population difference. This figure is from Ref. [20]. 
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Figure 17: (Color online) Spin noise tomography and reconstructed Wigner function of the spin-squeezed BEC [20]. a, Observed 
spin noise for the spin-squeezed state (filled circles) and for a CSS (open circles). In the squeezed state, a spin noise reduction 
of — 3.7± 0.4 dB is observed for m i n = 6°, corresponding to gfj = — 2.5 zfc 0.6 dB of metrologically useful squeezing for Ramsey 
spectroscopy [20]. The blue line describes the simulation results which take account of particle loss but not technical noise, 
give rise to a minimum variance of -12.8 dB. The technical noise is included in the red line, which is in good agreement with 
the experimental data, b, Reconstructed Wigner function. The black contour line indicates where the Wigner function has 
fallen to X/i/e of its maximum. For comparison, the circular contour of an ideal coherent spin state is shown. This figure 

is from Ref. [20] . 
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potential. The hyperfine states |0) = \F = 1, mp = — 1) and |1) = \F = 2, mp = 1) of 87 Rb forms a two-level 
system. The system Hamiltonian is the same as in Eq. (353), while here the notation of the spin operator J 
is replaced with S. A key problem in this experiment is also to produce nonvanishing nonlinear interaction 
strength, since the s-wave scattering lengths here satisfy 

a 00 ■ a i : a n = 100.4 : 97.7 : 95.0. (359) 

They realized x — 0.49 Hz though spatially separating the two components of the BEC with a state- 
dependent potential, during which the U ab term in Eq. (153) is zero. The Ramsey process here is similar to 
that in Ref. [19], and is displayed in Fig. 16. At first, a CSS is prepared by a 7r/2 pulse of duration 120 j«s. 
After that, the two components are separated, and the nonlinear interaction is active for a best squeezing 
time to avoid 'oversqueezing'. The two states overlap again after 12.7 ms, and the nonlinear interaction 
stops. As shown in Fig. 16, the transverse spin component 

S e = S z cos9- S y sm9, (360) 

is measured by rotating the state in the y-z plane by 9 before detecting S z . The experiment results are shown 
in Fig. 17, where comparisons between squeezed and coherent states are presented. The normalized variance 
A n Sg — 4ASg/(N) is plotted so that A n Sg = OdB corresponds to the shot noise limit. For squeezed state, 
a minimum of 

A n Sf = -3.7±0.4dB (361) 

is found at m ; n = 6°, with interference contrast C = 2\(S X )\/N = (88 ± 3)%, indicating a squeezing 
parameter of 

= -2.5±0.6dB, (362) 

and reduction of spin noise by — 3.7±0.4 dB was the observed. This could be used to improve interferometric 
measurements by —2.5 ±0.6 dB over the standard quantum limit. The experimental scheme and results [20] 
are summarized in Figs. 16 and 17. It is pointed out that, the technique for tuning of interactions in a BEC 
through wavefunction overlap also works in magnetic traps and for atomic state pairs where no convenient 
Feshbach resonance exists. 

8.2. Squeezing transferred from light to atoms 

Squeezing of photons is not easy to preserve, hence researchers considered how to transfer squeezing from 
photons to atomic ensembles [50, 112, 51, 16, 54], which are more convenient for storing quantum information. 
Many efforts have been devoted to use atomic ensembles to realize quantum memories [48, 287, 288], and 
using photons as information carriers. A practical analysis of generating spin squeezing in an ensemble of 
atoms, by absorption of squeezed light, was studied theoretically in [51], and soon after it was realized in 
experiments [16]. Spin-squeezed states can also be generated via electromagnetically induced transparency 
[64, 66, 67, 68, 69]. In this section, we present the experimental [16] results for generating spin-squeezed 
states when the atomic ensemble absorbs the squeezed vacuum field of the photons. 

8.2.1. Theoretical proposals 

Here discuss the theoretical proposal [51] of creating spin-squeezed states in atomic ensembles via the 
absorption of squeezed vacuum, which is illustrated in Fig. 18. The experiment is realized in Ref. [16]. We 
consider a cloud of Cs atoms [51]. The energy level structure is shown in Fig. 18, where we only consider 
levels |65i/2, F — 4), |6P 3 / 2 , F = 5) and |6P 5 / 2 , F = 6). There are two transition configurations: A V-type 
configuration consisting of levels |0) (|6Si/ 2 , m = 0)), |1) (|6P 3 / 2 , m = —1)) and |2) (|6P 3 / 2 , m = 1)), 
which is used to generate spin-squeezed states; a A-type configuration consisting of levels |1), |2) and 
\^F> 5 /2, m = 0), which is used for polarization measurement. 

Below we mainly consider the y-type structure and only investigate the subspace spanned by levels |1) 
and 1 2). The transition between these two levels is forbidden due to symmetry. However, effective transition 
processes between these two levels are realized via the level |0), and the processes are described by the 
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operators F12 and F21. Therefore, we should first derive the dynamical evolution between levels |1) (or |2)) 
and |0), which involves some assumptions and the Heisenberg-Langevin theory [51]. 
The pseudospin operators are given by 

Fx — ^ (-^12 + ^21) 5 
Fy = 2i (yFl2 ~ F21 ^ ' 

F z = X - (F u - F 22 ) , (363) 

where the collective operators Fij = ^2 a ij 1 (h 3 = 1; 2), and fi is the particle index and = Note 
that, since we use the hyperfine states, the angular momentum operators are denoted by F instead of J. 

As shown in Fig. 18, a left-circularly-polarized coherent field and a right-circularly-polarized squeezed 
vacuum field propagate along the z-direction through a cloud of Cs atoms. Due to the selection rule, the 
coherent field is coupled to levels |0) and |1), and the squeezed vacuum is coupled to |0) and |2). The atomic 
decay rate 7 is assumed to be equal for the two upper states |1) and |2). To describe the interaction between 
field and atoms, Ref. [51] employs the continuous field operators a(r,t) and atomic operators <7y (r, t). The 
continuous atomic operators are given by [51, 289] 



aij m = p~Lv 5> p ^Pt ( z - zM )] 4- ( 364 ) 



where p is the atomic density, and the sum is performed over atoms enclosed in the volume AV around the 
position r. The volume AV is only a small fraction of the entire atom cloud interacting with the field. The 
frequency between levels |z) and \ j) is Uij = (Ei — Ej)/h. 

1/3 

Now we consider an important assumption A <C L ~ (AV) ' <C A' < I, where A is the wavelength 
of the field, A' is the length-scale of the field changes due to absorption and dispersion, and I is the length 
of the atomic cloud along the field propagation direction. This assumption ensures a slow-varying field 
equation [51, 289] 

19 d 



c dt dz 



a(z,t) = gp dx dy u(x, y) a i (r, t) , (365) 



which separates the transverse component of the field a (r, t). Here, u (x, y) is the transverse mode function 
satisfying the normalization condition J dxjdy u 2 (x,y) = 1. The interaction strength g is determined by 
the atomic dipole moment. 

The Heisenberg equation of motion for the continuous atomic dipole operator is given by [289] 

d -a i (r,t) = - ( ^7 + iu) i J ct i (f,t) + g u (x,y)a(z,t) [a n (r,t) - a 00 (r,t)} 



+ V7 Wn (r, *) - aoo (r, t)] b in (r, t) , (366) 

where 7 is the atom decay rate, 6 ln (r, t) is the continuous operator of the atomic noise. For different atoms, 
the noise operators are independent, and thus obey the commutation relations 

p [b in (r, t) , 6 in t {fxM)] =5(f-f f 1 )8(t-h). (367) 

Now we assume the field to be so weak that the number of excited atoms is negligibly small, and replace 
coo (f, t) — an (r, t) = 1 in Eq. (366). Then we obtain 



7^01 (r, t) = - Q7 + iu) i \ 001 (f, t) - g u(x, y) a(z, t) - ^7 b m (f, t) , (368) 
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Figure 18: (Color online) Schematic diagram for a polarization measurement [51]. Three input beams propagate through a 
Cs atom cloud along the ^-direction, a 7r/4-polarized probe beam which is nearly resonant with 6P3/2 to 6-D5/2 transition, a 
left-polarized coherent beam that induces transition [51] between levels |0) (6S1/21 fn = 0) and |1) (6P3/21 m = — 1)> an d a 
right-polarized squeezed vacuum which coupled with levels |0) and |2) (6P 3 /2i m = !)■ The coherent field is much stronger 
than the squeezed vacuum. The atom cloud is optically thick for the squeezed vacuum and coherent field, and we assume that 
these two fields are completely absorbed, while the probe beam is only partially absorbed. Due to the population differences 
between levels |1) and |2), the refraction index for left- and right-polarization components of the probe beam are different, thus 
the polarization plane of the probe beam is rotated along the z-direction. After a polarization beam splitter (PBS), the probe 
beam is split into x- and y-polarized beams and is detected by two photon detectors. As the probe beam is it /A polarized, for 
the case of no atom clouds, the photocurrent difference is zero. Therefore, when an atom could exist, the photocurrent different 
is measured to gain information about the population fluctuations of levels |1) and |2). As the absorption is complete, a 50% 
degree of spin squeezing can be obtained. 
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The solutions of Eqs. (368) and (365) can be obtained [51] by using the Fourier transform of ctqi, 6 ln and 



ft. 



ctoi (f, w) 



3 ti (x, y) a («, w) + 6 m (r, w) 



7/2 -« (w - Woi) 
a(*,w) = a(0,w)e ifc ^ z - 

with the new dispersion relation 



77 g p e *(")(*-*i) 
7/2 -l(U>- U>oi) 



ig 2 p 



c 7/2 - i (u - W01) ' 
The redefined continuous noise operators become 



6 m (z,w) = / / « (as, 6™ (r, w) , 



(369) 
(370) 

(371) 
(372) 



with commutation relation p [6 m (z,u>) , b lnj; (zi,u>i)] — S (z — zi) 5 (w — u>i). For a specific field, the expec- 
tation values of Eq. (370) can be calculated, and inserted into Eq. (369) to obtain the dynamics of the local 
dipole operators. Finally, by using Eq. (364) the collective atomic dipole of the sample can be obtained [51]. 
The above results can be used to compute the collective operators 



(*0 

ij 



p / draij (r) , = 1,2) 



(373) 



Assume now that the coherent field cii and squeezed vacuum a 2 have an identical transverse mode function 
u (x, y), and also that the coherent amplitude is much stronger than the squeezed vacuum fluctuations. Thus 
the operator products can be linearized as 



0-12 (r, t) = (a 10 (r, t)) a m (f, t) , 
a 12 (r*, A = oj - ojx) = (cr 10 (r, t = 0)) <t 02 (f, u) , 



(374) 



where (<7io (r,t)) = (a\o (r,t = 0)) exp(— iu>\t). Now we assume u>\ = loqi = W02, and consider the practical 
condition 

2 _ 7A 



9P» 



(375) 



which implies complete absorption, and that the light traveling time is shorter than the atomic decay time. 
Finally, the collective operator F12 can be obtained [51] 



F 12 (A) ~ — — [4 n («i + A) + 4 n (wi + A)] 

7 — til 

= ne-** [4 n (w! + A) + 4 n (wi + A)] , 



(376) 



where a' 2 n (t) = a\ n = 0, t), and «, are the argument and phase of a* / (7 — iA), respectively. The spin 
squeezing is obtained from 

(377) 
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Figure 19: (Color online) Schematic diagram of the experiment in Ref. [16]. The squeezed vacuum is y-polarized, while the 
pumping beam is x-polarized. The probe is also x-polarized, and after the A/2 wave plate, it turns to w/4 polarized. This is 
from Ref. [16]. 



For the vacuum input field the variance is _ = 1/4, and thus (i*^ v ) = \ (F z ) = \ (F) for coherent 
spin state [51]. For squeezed vacuum, as discussed in Sec. 2, 

4X 2 + = ex P(- 2 \V\), 4Xi = ex P (2 | ^7 1 ) , (379) 

where rj = \r)\ exp(i2</>) is the bosonic squeezing parameter. Therefore, if the input field is a broadband 
vacuum, in the limit X\ = 0, then (F%) = (F z ) /4 and the spin-squeezing parameter becomes (see Table 1) 

= \ (380) 

Thus for an optically thick atomic ensemble, it is possible to achieve 50% spin squeezing when the squeezed 
vacuum is completely absorbed. 

8.2.2. Experimental realizations 

Soon after the above theoretical proposal [51], an experiment was carried out in Ref. [16], with the 
experimental setup shown in Fig. 19. They observed a spin-squeezed ensemble of 10 7 cold Cs atoms, 
which is produced via transferring a state of free propagating squeezed light to the atomic ensemble. This 
experiment also demonstrated storage of quantum information of light in atoms. In this experiment, they 
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Figure 20: (Color online) Schematic diagram of the electric field for the mixed output beam after PBS1 in Fig. 19 [54], The 
fluctuations in the composite field of the coherent pumping beam and the squeezed vacuum depend on their phase difference. 
If the phase difference 6 = 0, the so-called out-of-phase case, the field will have reduced fluctuations in the polarization axis 
direction, while if 6 = 7r, the so-called in-phase case, the field will have reduced fluctuations. This figure is from Ref. [54]. 



prepared about 10 Cs atoms in a magneto-optical trap, and then excited the 6P collective spin with 
the coherent pump polarized along the x-direction to establish the quantum limit of the spin noise. The 
coherent pump is of 50 mW with a beam diameter of 4.0 mm, resonant with the 852 nm transition between 
I6S1/2, F — 4) and I6D5/2, F — 6), provides weak excitation for the 6P collective spin of about 10 7 atoms. 
Then the atom states are detected by a probe beam, which is linearly polarized parallel to the coherent 
component of the pump and near resonant with the |6P 3 / 2 , F = 5) — > \6D 5 / 2 , F = 6) (917 nm) transition. 
The properties of atoms arc related to the change in the probe differential photocurrent Si. In the absence of 
atoms, the differential photocurrent is zero on average, and fluctuates due to shot noise. In the experiment, 
the fluctuation of the change of the differential photocurrent [52] 

(Si) 2 (A) = - [1 - exp {-a A )} + s a exp (-2a A ) (sj^j ^ (381) 

is measured, where (SJ) 2 is the atomic noise contribution per atom depending on the geometry of the 
experiment, s is the probe saturation parameter, and 

(7/2) 2 o ( 7 /2) 2 

Qa = ao 9 = -r^V 9 (662) 

A 2 + ( 7 /2) 2 A A 2 + ( 7 /2) 2 

is the probe optical depth at the detuning A [16, 76], with TV the number of atoms interacting with the 
probe, a the atomic cross section for unpolarized light, A the cross section of the probe beam, and 7 the 
atomic decay rate. 

The probe beam diameter is 3.7 mm and the power is about 250 mW for the off-resonant spin squeezing 
measurements s — 0.4. The squeezing transfer is performed via injecting squeezed vacuum which is linearly 
polarized along the y-direction. After the first polarized beam splitter (PBS1), the coherent pumping beam 
and the squeezed vacuum are mixed, and thus the output field will have fluctuations in either the polarization 
axis direction or in the ellipticity (depending on the phase as shown in Fig. 20) reduced below the standard 
quantum limit [16]. 

In the case of squeezed light, the fluctuations in the intensity difference between the right- and left- 
polarized components of the pump field is 

AS 2 Z = - (i?L cos 2 9 + R+ sin 2 0) , (383) 

where S z is the Stokes operator (E.3), 9 is the phase difference between the squeezed vacuum and the 
coherent pumping field, and R± are determined by the efficiency of the squeezing source, for coherent states 
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Figure 21: Squeezed spin noise of atoms [16]. Dashed line at zero: spin noise for uncorrelated atoms. Dots (left axis): 
antisqueezed spin noise. Triangles (right axis): squeezed spin noise. Solid line: square Doppler broadened dispersion function. 
Dotted line: expected spin squeezing spectrum. The detuning uncertainty is 0.5 MHz. Note that the scales for spin squeezing 
and spin antisqueezing are different. This figure is from Ref. [16]. 



R± = 1. Therefore, if R- < 1 < R+ and in the case of 9 = 0, the output field has ellipticity fluctuations, and 
the fluctuations of the intensity difference between the right- and left-polarized components of the output 
field are reduced below the shot-noise limit [16]. Intuitively, this reduces the fluctuations in the population 
difference between the +m and —m Zeeman levels. 

With the squeezed vacuum in phase, the average squeezing of the pump light available at the trap site 
is (— 1.8 ± 0.2) dB. The squeezing parameter is defined as [16] 



(«)cdh 



g= v " ? ' , (384) 



where Si 2 is the change in the probe differential photocurrent noise caused by the atoms. The subscripts 
"sq" and "coh" denote the spin-squeezed state and CSS, respectively. The best squeezing shown in Fig. 21 
is —3.0% at ±6 MHz. The average value for all of the 183 individual measurements [16] at ±6 MHz is 
£ = — (1.4 ± 0.4) %. The actual degree of squeezing should be greater, since the relative contribution of the 
readout efficiency is not known from the overall efficiency r\. 



8.3. Generating spin squeezing via quantum nondemolition measurements 

Quantum nondemolition measurements have been studied to generate non-classical photon states [170], 
and also spin-squeezed states [74, 76, 53, 75, 77, 290, 78, 80, 82, 40, 42, 92, 23, 96, 25, 94, 98]. At first, we 
briefly introduce the criteria for QND measurements. Then we present a typical QND-type Hamiltonian, 
implemented via dispersive interactions between light field and atoms. Utilizing this Hamiltonian, spin 
squeezing was generated in experiments [76, 23, 25]. In this QND-type Hamiltonian, the light field plays the 
role of a probe, by readout of the light field, the atomic spin state could be squeezed, but is conditioned to 
the specific readout of the light field, and thus is called conditional spin-squeezed state. By using feedback 
techniques [53, 290, 235, 78, 73, 99], unconditional spin-squeezed states can be generated. 

8.3.1. Theoretical background 

We first briefly discuss the basic requirements for a QND measurement (see, e.g. Ref. [170]). Consider 
two systems A and P, where A is the signal, and P acts as a probe coupled to the signal. By detecting the 
probe (e.g., measuring a physical quantity P a £ {Pi} that belongs to P), we could gain information about a 
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signal operator A a € {Ai} without perturbing its subsequent evolution. The two sets {Ai} and {Pi} contain 
conjugate operators of the probe and signal systems, respectively. For example, for a free particle, the two 
conjugate operators could be the position operator x and momentum operator p. The total Hamiltonian for 
the signal-probe system reads 

H = H A + Hp + Hi. (385) 

The conditions for a QND-type measurement of a signal operator A a via detecting a probe operator Pp are 
listed below: 

(i) For a QND-type measurement of A a , the interaction term Hi must contain the operator A a . 

(ii) Backaction-evasion in the measurement of the operator A a requires 

[Hi, A a ] = 0. (386) 

(iii) The dynamics of Pp should be associated with A a , so that the readout of Pp contains information 
about A a ; therefore 

[Hi, Pp] ± 0. (387) 

(iv) Since the successive measurements of A a should be predictable, then 

[H A , A a ] = 0. (388) 

From the criteria given above, one of the simplest QND-type Hamiltonians is 

Hi = A a P a , (389) 

and the probe operator to be detected is Pp with f3 =^ a. This QND-type Hamiltonian can be implemented 
through far-off resonant dispersive interactions between light and atomic ensembles, as shown in Fig. 22(b). 
Two probe light beams, ai and a 2 , are coupled to different hypcrfine levels, with detuning Ai and A 2 , 
respectively. In the large detuning case, the effective Hamiltonian is [23] 

5?A, _ t _ T _ to ( a T , 1 5 

4' 



H = hJ2 2/4 T \ 2 a \ ^ J ^ = m ( SzJz + jnN ) , (390) 

i=l I 1 



where the 74 's are the spontaneous decay rates for the excited levels, and 

J ll =Y,\^^\ ( 391 ) 

measures the population of atoms in the level and n, N are the photon and atom numbers, respectively. 
The Stokes operator S z (E.3) measures the population difference between the left- and right-polarized states. 
In experiment [23], the detuning Ai and A 2 are chosen in such a way that 

rfAx _ ,1A 2 _ n (392) 



7l 2 /4 + A2 Tl/4 + Al 2- 

This Hamiltonian acts as a Faraday rotation to the photon polarization, and as a fictitious magnetic field 
to atoms. The Heisenberg equations for the operators S and J are 
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Figure 22: (Color online) (a) Schematic diagram for the QND measurement. The input probe beam is x-polarized, and the 
atoms are initially polarized along the ic-axis. The two circular polarized components of the probe beam interact dispersively 
with atoms, with detuning Ai and A2, as shown in (b), f2i and O2 are atom-field couplings. The dispersive interactions 
induce a rotation of the polarization plane, the Faraday rotation. The A/2 wave plate together with the following PBS and 
two detectors play the role of QND measurement of J z by detecting S y , which is readout by the photocurrent difference, (c) 
Joint probability distributions for a system with N = n = 30 and \ = 0.1. (d) When the readout of S y gives m = —7, the 
conditional squeezed state is plotted and compared with the coherent spin state. 
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where x = ^ t. Therefore, a QND measurement on J z is performed through the readout of S x (t) [or S y (t)]. 

Now, we consider a specific case, where the initial states for the light field and atoms are both CSSs 
along the x-direction, 

\f (0)) = \S, S) X \J, J) x , (394) 

where S = n/2, J = N/2, and the subscript x means that the states are represented in the S x and J x basis. 
A QND measurement on J z is performed by detecting S y (or S x ), the readout of S y (t) becomes 



(S y (t)) = ~{sm( X J z )) 



(Sy) 



nx 



(Jz 



(395) 



where the average (•) is carried out under the initial states since we are using in the Heisenberg picture. 
Thus we see that (J z ) and (S y ) are related, the mean value of (S y (t)) is approximately zero. From (395), 
the fluctuation of S y (t) is also related to J z [74]. 

Below, we demonstrate the QND-measurement-induced spin squeezing from the point of view of state 
evolution, shown in Fig. 22(c)-(d). As presented in Fig. 1, for a coherent spin state \J,J) X , if we measure 
J z , the probability of the readout of states \J,M) Z obeys a binomial distribution 



Pj (M ) = P[\J, M) z | J, J) x = Z (J, M\J, J) 



1 

2N 



N 

N/2 + M 



where P (jJ,M) z \\J, J)x) is a conditional distribution. Similarly, in the basis of \S, 
basis, we have 



P s (m) = p[\S,m) y \\S,S) x ) = | y (S,m\S,S) x 
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(397) 



The evolution of the state becomes 

lv(*)) 



exp (-ixS z J z ) \S, S) X \J : J) x 

exp {-i X MS z ) y/Pj{M)\S, S) X \J, M) 2 
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At time t, we measure S y . For a single trial measurement, the readout will fall in an eigenvalue m of 
S y , and the photon state collapses into the corresponding eigenstate \S, m). After this readout, the atomic 
state is also changed according to the measurement result of S y . Therefore, the focus is on the conditional 



distribution P ( | J, M) z \ \ S, 



which gives information of the atomic state after a single trial measurement 



of S y . To analyze this problem, we consider the system size to be very large, in which case the binomial 
distribution becomes a Gaussian normal distribution, 

1 
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With the Bayes theorem 



P (| J, M) z \m) P s (m) = P (\S,m) y \M) Pj (M) , 



(399) 



(400) 



the conditional probability distribution for P (\S, m) y \M} is easy to obtain. If the atomic system is in the 
state | J, M) z , according to Eq. (398) the photon state is rotated around the z-axis with an angle 4> = %M. 
If the rotation angle is very small, then the photon state after the rotation is still a Gaussian distribution 
with mean value at sin (<fi) n/2 ~ <t>n/2, and 



P(\S,m) y \M) = ^_ exp 

y/TTn/2 
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(m — (f>n/2) 2 
n/2 



(401) 



Therefore, the conditional probability distribution of the atomic state becomes 



P(\J,M) z \m) 
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where 
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and the state of the atomic system after detecting S y is squeezed, while the expectation value of J z is shifted 



Note that the spin length is generally shortened due to photon scattering and other noises during the 
QND measurement, so the above Eq. (403) holds under the condition that the spin length is not affected 
too much after the QND measurement. An illustrative description of the above procedure is displayed in 
Fig. 22. The spin state derived above is the so-called conditional squeezed state, which is conditioned to 
the first measurement output of S y . In a "single trial" measurement of S y , the output value is random with 
respect to a specific probability distribution. The conditional squeezed state is generated randomly, so it is 
not repeatable. To obtain an unconditional squeezed state, we could employ a feedback or filter scheme to 
shift (J z ) back to zero, which is demonstrated in Refs. [53, 290, 235, 78, 73, 99]. 

8.3.2. Experimental results 

Below, we first give a brief introduction about some experimental works, then we present a detailed 
discussion about a recent notable improvement. In Ref. [76], QND measurement was performed on a 
collective spin of ~ 10 6 Cs atoms, and sub-shot-noise fluctuations of the collective spin were observed. In 
this experiment, the probe beam was short pulses of laser light that linearly polarized at ±45°. Then, in 
Ref. [77], the authors performed continuous QND measurement on a collection of ~ 10 7 Cs atomic spins with 
an off-resonant laser beam. The continuous QND measurement is an important feature of this work, and 
the probe laser also serves for measurement of the spin rotation signal. They observed squeezed atomic spin 
states with spin noise reduction to 70% below the shot-noise limit, and measured a small spin rotation with 
an accuracy exceeding the shot-noise limit of the phase measurement. Reference [40] proposed a sequence 
of QND measurements and spin rotations, which require interferometric QND measurement of the atomic 
population with a sensitivity at the projection-noise level. To achieve this they demonstrated experimentally 
that a shot-noise- limited fiber optical interferometer at the white- light setting can reach a sensitivity sufficient 
to detect the projection noise under conditions not far from those of the QND measurement. Increasing the 
optical density can significantly improve the QND figure of merit k 2 , and thus will enhance the degree of 
squeezing. In Ref. [92], the authors demonstrated QND measurement of atomic ensembles in a high-finesse 
optical cavity, which can enhance the optical depth of the atomic sample. In this experiment, a cloud of 87 Rb 
atoms was cooled in a MOT, located at the center of a hemispherical cavity. The maximum atom-cavity 
coupling, gj (2tt) = 53 kHz, is realized for the \F — 2, rap = 2) — > \F' = 3, nip = 3) cycling transition. The 
atoms are released from the MOT, further cooled by optical molasses, prepared in an equal superposition of 
the \F = 1, rap — 0) and |F = 2, mp — 0) clock states, and probed by a standing wave of intracavity light. 
They observed 3.8 dB in the variance of J z below the projection noise with spontaneous emission of < 30%. 

Reference [103] demonstrated sub-projection-noise sensitivity of QND spin measurements of 1 x 10 6 
cold 87 Rb atoms in a broadband atomic magnetometer. The high-bandwidth system avoids decoherence 
effects through rapid probing. The spin readout noise is 2.8 dB below the thermal spin noise and 1.6 dB 
below the CSS projection noise. Since the Hamiltonian (390) only appropriate for spin-1/2 or other two- 
level systems as pointed out by studies of Ref. [104], where they considered larger spins and proposed a 
decoupling mechanism based on a two-polarization probing technique to recover the ideal QND behavior. 
The decoupled QND measurement achieves a sensitivity 5.7(6) dB better than the spin projection noise. 



to xm£ R N/2. 
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Figure 23: (Color online). Experimental schematic diagram of Ref. [102]. (a) Experimental setup, (b) Level structure of 
rubidium atom, (c) Experimental sequences. Dashed lines denote the microwave pules while solid lines denote the probe 
pulses. The fr pulse represents a composite 7r pulses. Various operations are performed between two readouts (Mi and M2), 
denoted by X, such as the measurement of the variance of the CSS, the operation of a clock. A-C: Semiclassical probability 
distribution functions of the state. This figure is from Ref. [102]. 

Recently, Ref. [102] demonstrated the generation of spin-squeezed state of 5 x 10 4 trapped 87 Rb atoms 
on an atom-clock transition by resonator-aided QND measurement. For photon number p = 3 x 10 5 , 
they observed noise reduction of 3.0(8) dB below the shot-noise limit. The hyperfme clock states used 
are |1) = |5 2 Si /2 , F — 1, mj? = 0) and |2) = |5 2 Si /2 , F = 2, m F = 0). They used the metrological 
spin-squeezing parameter defined as 

2S iD (A^) 2 

Cm = , . 2 ' ( 405 ) 



s 

where S m is the initial spin of the uncorrelated ensemble. Actually, this parameter is equal to £fg when 
= N/2. The experimental setup is shown in Fig. 23. A strong ensemble-light coupling was achieved by 
using a near-confocal optical resonator, and the experiments are performed on an ensemble containing up to 
N a = 5 x 10 4 laser-cooled 87 Rb atoms optically trapped inside the resonator in a standing wave of 851-nm 
light. The wavelength of the probe light is 2ir/k = 780 nm, and the cavity has a finesse T = 5.6(2) x 10 3 , a 
linewidth n = 2tt x 1.01(3) MHz, and a mode waist w — 56.9(4) //m at the atoms' position, corresponding 
to a maximal single-atom cooperativity r/o — 24J 7 / (7rfc 2 u> 2 ) = 0.203(7). One resonator mode is tuned 
3.57(1) GHz to the blue of the \5 2 S 1/2 ,F = 2) -> \5 2 P 3/2 ,F' = 3) transition in 87 Rb to result in a mode 
frequency shift u> proportional to the population difference = A 2 — Ni between the hyperfme clock states 

\l) = \5 2 S 1/2 ,F = l,m F = 0), 

\2) = \5 2 S 1/2l F = l,m F = 0). (406) 
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Figure 24: (Color online). Normalized spin noise cr 2 = (A5^) Mi / (A5^) cg g (open diamonds), contrast C (open square), and 
metrological squeezing parameter £ m (solid triangles). Solid line denotes the fits. This figure is from Ref. [102]. 



= CSa. as shown 



The experimental processes are shown in Fig. 23(c). At first, the atomic ensemble is prepared in state 
|1), after a tt/2 pulse, the atomic state polarizes in the x-axis, then after measurements Mi, the atomic 
ensemble is conditional spin squeezed as Fig. 23C. To quantify the spin squeezing, they measured AS^ and 

(^S^ , the latter is obtained from the observed contrast C of Rabi oscillations since (§ 

in Fig. 24, from which the metrological squeezing parameter achieves 

C 1 - 3.0(8) dB 

for photon number p = 3x 10 5 . 

The QND measurements discussed above were all performed on alkali atoms, which have complicated 
multilevel structures. In Ref. [25], the authors demonstrated a QND measurement with a collective spin of 
cold ytterbium atoms, 171 Yb, which have the simplest ground state with a nuclear spin of one-half and have 
no electron spin, and thus the system is robust against stray magnetic fields since the magnetic moment of 
nuclear spin is a thousandth of that of electron spin. By using short light pulses with a width of 100 ns, more 
than a hundred time operations are expected to be performed within the coherence time, and experimental 
results showed 1.8^5 dB spin squeezing. 

With the above experimental efforts, recently, two successive notable improvements in QND experiment 
were reported [23, 24], where a complete sub-shot-noise Ramsey process was demonstrated by means of 
QND measurement. The experimental schemes for these two works [23, 24] are similar, thus here, we mainly 
discuss the very recent one [24]. The experimental setup is shown in Fig. 25(a). An ensemble of ~ 10 5 cold 
Cs atoms are located in one of the arms of the Mach-Zehnder interferometer. The atoms are first loaded into 
a standard MOT, and are then transferred into a far-off resonant optical dipole trap (FORT), generated by 
a Versadisk laser with a wavelength of 1032 nm and a power of 2.3 W. After loading the FORT, the MOT 
is turned off and a 1.22 Gauss magnetic field is applied, which defines a quantization axis orthogonal to the 
trapping beam. The ground levels 

||) = \6S 1/2 ,F = 3,m F =0), 

\V = \6S 1/2 ,F = 4,m F =0), (407) 

are referred to the clock levels, and by optical pumping, the atoms are prepared in the clock level | |). 
As described in Ref. [23], atoms remaining in states other than | \) due to imperfect optical pumping are 
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(a) Experimental setup 



(b) Cesium D 2 -line level diagram 




(c) Pulse sequence and noise data (d) Projection noise and spin squeezing 




Figure 25: (Color online) (a) Schematic diagram of the QND measurement experimental setup, a Mach-Zehnder interferome- 
ter [24, 23] . The elongated Cesium atomic cloud is placed in the probe arm. The probe beam consists of two identical linearly 
polarized light beams Pf and Pj_, which acquire phase shifts proportional to the number of atoms in the clock states JV-j- and 
respectively, as also shown in (b). In (c) A, we show an experimental sequence for generating the state \J, J) x and QND 
measurement of S v , and in B and C, noise data of <f>2 and <pi show correlations between <f>2 and <f>i, and reduced noise of 
(<j>2 — C^i)- 1° (A), we display the experimental results of the projection noise and spin squeezing. Blue stars: variance of the 
second measurement A 2 02 ■ Solid blue line: quadratic fit to A 2 (j>2 ■ Green area: atomic projection noise of the CSS. Red 
diamonds: conditional variance A 2 (<j>2 — C<f>i)- Red line: reduced noise as predicted by the fits to the noise data. This data 
was obtained by acquiring 4800 experimental runs. Inset: spin squeezing £|j as a function of the number of pulses combined 
to form the second measurement. The data is fitted with an exponential decay (solid line). Figures (a), (b), and (d) are from 
Ref. [24], and figure (c) is from Ref. [23]. 
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subsequently pushed out of the trap. To prepare the CSS 10 x 1 



To (I i) + I t)) ! a resonant 7r/2 microwave 



^2 



? = 1 

pulse at the clock frequency is applied. Then, successive QND measurements of the population difference, 
N-f — N±, will be performed by measuring the phase shift of the probe light in a balanced homodyne 
configuration. After the QND measurement, all atoms are pumped into the F = 4 level to determine the 
total atom number Na- The sequence, shown in Fig. 25(c), is repeated several thousand times with a cycle 
time of ~ 5 s. 

As shown in Fig. 25(a), two identical linear polarized beams P-j- and Pj., which are generated by two 
extended-cavity diode lasers, enters the interferometer from the same port [Fig. 25(a)], off-resonantly probe 
transitions \F = 3) — > \F' =4) and \F = 4) — > \F' = 5), respectively. Each beam gains a phase shift 

<t>f(l) = kf{i) (408) 

which is proportional to the number of atoms N-f and iVj. corresponding to the clock states | f) an d | I) , 
respectively. The detuning and are tuned to ensure = k\,, such that the Hamiltonian for this 
system is Eq. (390), and thus 

J z = (N t - Ni) /2, (409) 

where iV-j-m are the photon number operators, Na = + N± is the total atom number. The photocurrent 
difference is equivalent to measuring Sy jOU t, which is derived from Eq. (395) as 

71 

Sy,out — Sy t i n + — XJz,in- (410) 

The phase difference between the two arms of the Mach-Zehnder interferometer is defined as 

<f>=^^, (411) 
n 

where n — + is the total photon number. In experiments [23, 24], the photon numbers of the two 
probe beams and are set to be almost equal. If the probe beam is in a coherent state, the variance of 
the phase difference is 

A 2 0=-+ x 2 (AJ,, in ) 2 , (412) 
n 

which is related to both photon shot-noise 1/n, and atomic spin noise (AJ^^) 2 . As the atomic ensemble 
is prepared in a CSS, AJ z in = i/Na/2. From Fig. 25(c), two successive measurements of (f>\ and <j>2 are 
performed, using probe photons n\ and ni- Thus, 

A^ 1(2) = — +X 2 ^, Cov(^> 2 ) = x 2 ^. (413) 
"i,(2) 4 4 

The covariance (413) is calculated for the atomic ensemble, which is the same for the two measurements, 
while the two probe beams have no correlations. As <p2 is conditioned to 4>i and can be predicted with 
a better precision, after the measurement of cf>±, an appropriate estimate of 02 is chosen as C^U anc i ^ ne 
variance of their difference becomes 



A 2 (0 2 -OM = (A 2 0i)C 2 -2 Cov(0 1 ,0 2 )C + A 2 </) 



2 



1 1 9 N A / \ 

n-i 1 k 4 

where the second line is the minimum of A 2 (<j>2 — C4>i) a t 

* A 2 4>i i + « 2 ' 

Here, 

K 2 = ni N x 2 /4:, (416) 
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which is the same result derived in Eq. (403). Now, according to Eq. (414), the spin-squeezing parameter 



& = ih' (417) 

A measurement of k 2 = 1.6 for Na — 1.2 x 10 5 atoms gives rise to a reduction of the projection noise by 
— 4dB compared to the CSS projection noise, displayed in Fig. 25(d). 

Now consider the decoherence during a QND measurement. In experiments [23, 24], the employ of 
two-color QND can suppress the dephasing effects induced by an inhomogeneous AC stark shift due to the 
transverse intensity profile of the probe beam [96] . Otherwise, since the probe interacts dispersively with 
the atoms, decoherence is inevitable induced by spontaneous photon scattering [23, 24], which shortens the 
mean-spin length as 

(J) (1-7?) (J), (418) 

where n = 1 — exp (na) , with n the total number of photons in the probe pulse and the parameter a = 
—2.39 x 10 -8 [24]. The coherence between the clock levels decreases while the phase shifts are almost 
unaffected, so in the absence of the classical noise, we have 

& = «- . (419) 

R (I-77) (1 + k 2 ) 

When using ni/2 = 3 x 10 6 photons per probe, the reduction is r\ = 14%, and experimental data gives 
= 2.7 dB as shown in Fig. 25(d). This is very close to the theoretical prediction from Eq. (419), 
— 2.8 dB. The above results indicate that experiments must carefully choose photon numbers for 

optimal spin squeezing. 



9. Conclusion 

We have reviewed both basic notions and recent progresses on spin squeezing. After a short introduction 
on bosonic squeezing, in Sec. 2, we reviewed various definitions of spin-squeezing parameters and their phys- 
ical significance. In order to understand these parameters better, we derived explicit analytical expressions 
of them for states with parity. At the end of Sec. 2, we showed that, in the limit of a large number of atoms 
and small number of excitations, spin squeezing reduces to bosonic squeezing. 

In Sec. 3, we reviewed theoretical approaches toward producing spin-squeezed states via one-axis twisting 
and two-axis twisting Hamiltonians. The one-axis twisting Hamiltonian has attracted considerable atten- 
tions, and has been implemented in BEC, as reviewed in Sec. 8. For this model, we presented analytical 
results for the spin-squeezing parameters £| and which scale with the particle number as I/A 2 / 3 . 

The two-axis twisting Hamiltonian can produce more squeezing than the one-axis case. Moreover, its 
optimal squeezing angle does not vary with time. However, analytical results are only available for par- 
ticle number N < 3, and numerical results showed that £| and scale as 1/N. This two-axis twisting 
Hamiltonian can be realized in spinor BEC. 

Section 4 showed that spin squeezing has very close relations to both entanglement and negative pairwise 
correlations. A system consisting of spin-1/2 particles is entangled if £^ < 1. For states with parity and 
exchange symmetry, close relations between the spin-squeezing parameters £g, fjg, and the concurrence C 
have been established. In this case, spin squeezing according to £fj is qualitatively equivalent to the existence 
of pairwise entanglement. Moreover, recent results generalized the proposal of using spin squeezing to detect 
entanglement, and also developed a set of spin inequalities. In the end, the two-mode spin squeezing is 
reviewed. It can detect the entanglement between spin variables of two separated atomic samples, and 
was proposed to be a valuable resource of quantum information processing. As explained in this section, 
spin-squeezing-based entanglement witnesses should be useful in measurements and are then expected to be 
of wide applications in experiments. 

In Sec. 5, we discussed the relation between spin squeezing and QFI, as well as their applications in 
quantum metrology. Applications of spin squeezing also deserve more considerations, and in Sec. 6, we 
explained how spin squeezing can be used to detect quantum chaos and QPTs. 



84 



In Sec. 7, we first consider how decoherence affects the generation of spin squeezing. If the decoherence 
effect is sufficiently weak, considerable large amount of squeezing can still be obtained. Then we showed 
that spin-squeezed states are more robust than the GHZ state in the presence of particle loss. The life- 
time of squeezing is discussed for three typical decoherence channels, and the spin-squeezing sudden death 
was derived and analyzed. These channels are prototype models, however, are quite general to describe 
decoherence processes. The results derived by using these channels still make sense when using concrete 
physical models. In addition, the one parameter in the dephasing, dissipation, or the depolarizing channels 
may be a very complex function depending on the realistic physical models of decoherence. Afterwards, we 
show that in the Ramsey spectroscopy, decoherence reduces the ability of spin-squeezed states to help in 
precision-improvement, however, compared with the GHZ states, spin-squeezed states are more robust to 
decoherence. 

In Sec. 8 we reviewed several recent experiments, as well as their theoretical backgrounds. Spin-squeezed 
states were realized in BEC and atomic ensembles. In the past two decades, spin squeezing has been 
studied extensively. It was first proposed in order to overcome the classical shot-noise limit and make high- 
precision atomic interferometers. This encouraged many theoretical and experimental efforts in the control 
and manipulation of atomic systems. Afterwards, it was realized that the notion of spin squeezing can 
provide an insight on entanglement, and this has attracted considerable recent interests. Experiments can 
now produce spin-squeezed states, and even use them to perform sub-shot noise interference. However, in 
practice, using spin-squeezed states in an atomic clock or a gravity interferometer still faces many challenges, 
including how to suppress the decoherence and enlarge the particle number of the spin-squeezed states. 

Finally, we summarize and list several main results, some of which are of broad interest: 

(a) There are several definitions of squeezing parameters (see Table 1), which one is preferable depends 
on various situations. For example, when studying the accuracy of Ramsey spectroscopy, the parameter (fa 
is more suitable. In this review, we mainly focused on two parameters, £| and £fj, which are widely studied 
in the literatures. 

(b) In general, spin-squeezed states are entangled; however this depends on specific squeezing parameters. 
From Sec. 4, we know that, for the parameter states are entangled if £^ < 1; while for the parameter 

states with parity have pairwise entanglement if £jg < 1. Generalized spin inequalities have also been 
developed to detect entanglement. 

(c) Since spin-squeezing parameters only involve expectations and variances of collective spin operators, 
they can be measured relatively easy in many physical systems, e.g., BEC and atomic ensembles, while 
addressing individual atoms is an arduous task. 

(d) Among proposals for generating spin-squeezed states, the twisting Hamiltonian, which describes 
nonlinear two-body interactions, plays an important role. This type of Hamiltonian has been studied and 
implemented in physical systems like BEC. 

(e) In parameter-estimation processes, compared with the GHZ state, spin-squeezed states can also attain 
the Heisenberg limit, but are more robust to decoherence. 

In summary, spin squeezing has provided very useful insight into quantum processes, and will continue 
to considerably impact quantum optics and quantum information science in the future. 
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Appendix A. Principal squeezing parameter 

From Eqs. (7) and (8), we have 

cos 2 9(AX) 2 + sin 2 9(AP) 2 + sin(20) Cov(X, P) 
(AX) 2 + (AP) 2 1 



Cb = mm 

9e[0,27r) 



mm 

0e[O,27r) 



(cos(2<9) [(AX) 2 - (AP) 2 ] + 2sin(26>) Cov(X,P)} 



i j Var(X) + Var(P) - \J [Var(X) - Var(P)] 2 + 4Cov 2 (X, P) 



where the variance 
and the covariance 



Var(X) = (X 2 ) - (X) 2 



Cov(X,P)= l -(XP + PX)-(X)(P). 



From the expressions for X and P in Eq. (2), one finds 

XP + PX a 2 - at 2 



which leads to the covariance 



From the identity 



we have 



which leads to 



Cov(X, P) = - Var(a) - Var(a t ) 
i L 

/ N \ N 

Var I ^ X* J = Var(Xi) + 2 ^ Cov(X l , X 



3J> 



Var(X) = Var(a) + Var(a f ) + 2Cov(a, 0+), 
Var(P) = -Var(o) - Var(a f ) + 2Cov(a, a)), 



Var(X) + Var(P) = 40^(0,0*), 
Var(X) - Var(P) = 2 (Var(a) + Var(a f )) . 
Substituting Eqs. (A. 5), (A. 9), and (A. 10) into Eq. (A.l), we finally obtain 



Cl = 2 | Cov(a, a)) - y / Var(a)Var(at) 



which is just Eq. (9). 



(A.l) 

(A.2) 
(A.3) 

(A.4) 
(A.5) 

(A.6) 

(A.7) 
(A.8) 

(A.9) 
(A.10) 

(A.H) 



Appendix B. Variances of Jfj for the CSSs 

We obtain the variance (AJ n ) 2 = (JH) - (J fl ) 2 for the CSS. The expectation values (j 2 ) are given by 



1 N 

( J l) = 7 (VioVjo 



N 1 

4+4 



JV 



N N 2 -N , 
1 + 4~ 



(B.l) 



8(> 



The last equality is valid since we are considering a product state with exchange symmetry. Note that 

X 
~2 



N 

J a ) = ^r {<?!*), (B.2) 



the variance is found to be 



(AJ Q ) 2 = ^(l-<<7 la ) 2 ). (B.3) 



Now, we consider a more general operator as 

N N 



J ft = J ■ ft = 

from which, its variance is given by 



-ft = l^T<J M} (B.4) 



2 ^ 2 

i=l i=l 



Var(J K ) = (AJ ri f 



j^Vax(a«) 

i—l 
1 N 



4 
i=i 



N 

T 



(l-K.) 2 ). (B.5) 



Here, the second equality follows from Eq. (A. 6) and the fact that the CSS is a product state and contains 
no correlations (the covariance Cov(Xi, Xj) = {[Xi, Xj] + /2) — (Xi)(Xj) — 0). The third equality is from 
the result of- = I. 

From Eqs. (B.4) and (100), we obtain 

N N 
( J n) = -r,Wifi) = y |(<?i)K • n, 

(AJnf = j (l - (a lfl ) 2 ) = f [l - I (^i)l 2 (no • n) 2 ] , (B.6) 
which is valid for any product state with exchange symmetry. The second equality is obtained from the fact 

(<?i) = |(<?i)|n = |(<?i)|(sin0cos(/>, sin0sin<^>, cos0). (B.7) 
From the definition of the CSS, the expectation values {J a ) (a = x, y, z) are obtained as 

N 

{J x ) = — sin6»cos0, 
N 

{J y ) = — sin 9 sin <j>, 
N 

{J z ) = -cos9. (B.8) 

with N — 2j, and 9, <j> are the polar and azimuth angles, respectively. Substituting Eq. (B.8) into Eq. (B.5) 
leads to the explicit form of the variances 

(AJ x f = sin 2 6> cos 2 <j>), 

(AJ y f = ^(1- sin 2 (9 sin 2 0), 

(A J z f = ^ sin 2 9. (B.9) 
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Using Eq. (1), one finds the variance along the MSD, 

N 



Var(J fto ) = -(l -!<<?!) | 2 ) 



For the CSS, only Dicke states are populated, and thus, 

(J 2 X + 4 + J 2 Z ) = j (y 



(B.10) 



(B.ll) 



From the above equation and from the relations between the collective operator J a and the local operators 
(B.4), one finds 



N 2 — N 

As the CSS is a product symmetric state, the above equation reduces to 

Applying this result to Eq. (B.10), one immediately has 

Var(J llo ) = 0. 

For the variances along the direction ri± perpendicular to the MSD, from Eq. (100), we have 

Var(J fl± ) = ^(l-( ( 7 1 ^) 2 ) 



N r 
4 

N r 

1 

N 

T 



l-((^i)-n±r 
l-K^l^no-nx) 5 

= j/2, 



(B.13) 
(B.14) 



(B.15) 



where the third equality follows from Eq. (40). 



Appendix C. The minimum value of the quantity n T Tn 

Here, we prove that the minimum value of n T Tn is a minimum eigenvalue of T . The minimum value of 
n T Tn is obtained by searching the vector n = (rti, ri2, nj) T . Since the matrix T is real and symmetric, we 
obtain 



min (n T Tn) = min \n T T (OTO T ) On 

ft V / n L 

where O is an orthogonal 3x3 matrix which diagonalizes T, and 

r d = diag(A 1 ,A 2 ,A 3 ), 



nin (n' T T d n'^j , 



(CI) 



(C.2) 



where we set Ai > A2 > A3. The rotated direction is n 1 — On — {n'^n^n 1 ^) . Now, the minimum value of 
n T Tn is obtained as 



nin (n T Tnj = min [n' T Y d n'^ = min (Xmf + A 2 n' 2 2 + X 3 nf^j = A 3 



(C.3) 



since n\ +n' 2 f + n'^ = 1. Thus, the minimum value of n Tn is obtained by choosing n' 3 — 1, and n[ = n' 2 = 0, 
therefore, n — T n' is the direction where the minimum is achieved. 



Appendix D. Expectations for the one-axis twisted state 



Here, we will use the Heisenberg picture to derive the relevant expectation values for the one-axis twisted 
state. To determine the spin-squeezing parameters £| and £fj, as seen from Table 2, one needs to know the 
expectation {<tu)i and correlations {01+02-) and (eri_<72-)- We first consider the expectation (cri 2 ). 

The evolution operator can be written as, 



U = cxp(—ixtJ x ) = exp 



V k>l 



(D.l) 



up to a trivial phase, where 9 = 2xt. From this form, the evolution of ju can be obtained as 

U^jizU = ju cos[9j^] + j lv sm[9ji% 

where 

Therefore, the expectations are 



l=k 



(ji»> = -^l / |cos[fljW]|l / > 



(D.2) 



(D.3) 



since (l|ji„|l) = 0. Here, |1'> = |1> 2 
values 



|1)jv- So, one can find the following form for the expectation 



(l|cos [0J B ]|1) = «l|e* J -|l)+c.c.) /2 

= (nf =1 (l|e^|l)+c.c.) /2 
= cos N {9'), 

where 9' = 9/2 and |1) = \1)® N . By using Eqs. (D.3) and (D.4), one gets 



JV-l to' 



(D.4) 



(D.5) 



Since the operator (J\ x &2x commutes with the unitary operator U, we easily obtain 

(Vlx&2x) = 0. 

We now compute the correlations {<Ju a 2z)- From the unitary operator, 

V^hzhzU = [j u cos(0ji 2) ) + j ly sin(0j< 2 ))] [j 2z cos[0(j lx + j x 3y )} + j 2y sm[9( Jlx + j^)} 
= ju cos{9j 2x ) cos{9j x V) - ju sm(9j 2x ) sin(9j x 3) ) 
+ji y sm(9j 2x ) cos(9jW) + jiy cos(9j 2x ) sm(9j x ^) 
x j 2z cos(9j lx ) cos(9j x 3) ) - j 2z sin(9ji x ) sin(9j x 3) ) 
+32 V sm(6j lx ) cos(9j x 3) ) + j 2y cos(9j lx ) sin(9j x 3) ) 



(D.6) 
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Although there are 16 terms after expanding the above equation, only four terms survive when calculating 
(hzhz)- We then have 

(jizhz) = (l|ji 2 j 2z cos 2 (0/2)cos 2 (^ 3 )) -j lz j 2x j 2y sin(9) S in 2 (dj^) 

+ ^kyhxhxky sin 2 (9/2) cos 2 [6 - hyhx32z sm(9) sin 2 (#j( 3 >)|l) 



;i'|cos 2 (0j( 3 ))|i'> 



1 • <-'-:^./? . 

cob w - 2 ((9)] , 



II') 



(D.7) 



where |1') = 1 1) 3 (g) ... <g> \1)n- The second equality in Eq. (D.7) is due to the property j x j y = —j y j x = 
iz/2, and the last equality is from Eq. (D.4). Finally, from the above equation, one finds 



JV-2 , 



(D. 



Due to the relation {<j\ x g% x + o\ y 02 y + oiz<?2z) = 1 for the initial state, the correlation (o\ y U2 y ) is obtained 
from Eqs. (D.6) and (D.8) as 

1 



(p ly a 2 y) = 2 (! 



„JV-2 , 



') 



Substituting Eqs. (D.6) and (D.9) into the following relations 

(rix<?2x + G \y G 2y = 2 ((Tl + (T2- + Ol-0"2+) 

leads to 

(o 1+ o 2 -) = \{l-cos N - 2 9), 

where the relation (0-1+02-) = {&1-&2+) is used due to the exchange symmetry. 
To calculate the correlation (<7i_cr2-), due to the following relations 

0"lx0"2x — V\y&2y = 2 (0-1 + 0-2+ + 01_<7 2 _) i 
i (<rix<72y + Ci y (J2x) = 2 (<7i + <72+ - , 

we need to know the expectations (j\ x j2y)- The evolution of ji x j2y is given by 



U ] si x s 2y U = ji x |j 2y cos 9(j lx + j x 3) ) 



J2z sm 



and the expectation is obtained as 



1 



(jlxhy) = [9{ ]lx +j x z) )\ |1') 



JV-2 rai 



Jixe 



3 e~ i6 -' k:c \l') 



cos- - (^)<l|jixe^ - 

W-2 



-cos 



(0')<l|ii*sm(0jix)|l) 



= ^sin (6»') cos^ 2 (9 r ). 



(D.9) 

(D.10) 
(D.ll) 



(D.12) 
(D.13) 



(D.14) 



(D.15) 



Here, |1') = |l)i O |1) 3 
symmetry, and thus, 



\1)n, where 1 1)2 is absent. Moreover, (j\ y 22x) = (jixj2y) due to the exchange 



{hx32y + hvhx) = 2 sin iff) cos N - 2 {&). 
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(D.16) 



For the state (322), we obtain the following expectations 

(a lx a 2y + a ly a 2x ) = 2 sin (0') cos 7 ^ 2 (0') . (D.17) 
The combination of Eqs. (D.6), (D.9), (D.12), (D.13), and (D.17) leads to the correlation 

(oi-ffa-) = - - (1 - cos^- 2 9)-- sin (0') cos N ~ 2 (0') . (D.18) 
8 2 

Substituting Eqs. (D.ll) and (D.18) to the expression £| in terms of local expectation values (Table 2) leads 
to the desirable expression of the squeezing parameter for the one-axis twisted state. 



Appendix E. Photon polarization and Stokes operators 

Here, we briefly introduce the Stokes operators for photons. The light field has two orthogonal polar- 
ization directions. If the light field propagates along the z-axis, we use the bases given by e x and e y , which 
denote the vertical and horizontal polarizations, respectively. The other two sets of bases can be expressed 
with e' x and e! y . If the phase difference between the e x - and ej,-polarized components is ±7r/2, then the light 
is 45° polarized, with new bases expressed as 

e' x = i (e x + e y ) , e' y = {-e x + e y ) . (E.l) 

If the phase difference between the e x and e y components is ±7r/4, then the light is circularly polarized, 
with right- and left-circular bases 

e+ = - — (e x + ie v ) , e_ = (e x - ie y ) . (E.2) 

The two polarized components can build up an angular momentum operator analogous to the Schwinger 
representation (151). Denote the x- and y-polarized field operators as a x and a y , respectively. Then, the 
Stokes operators are defined as 



cy ( a x a x a y a y) ) 



S z = g ( a + a + - a-a-^j ■ (E.3) 

The three Stokes operators (S x , S y and S z ) shown in (E.3) measure the photon number differences of the x- 
to y-polarization, 7r/4— to — 7r/4-polarization, and right- to left-polarization, respectively. The polarization 
squeezing of light field was studied in Refs. [291, 292, 293, 294, 295, 296]. 
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